This article deals with finite-time quenching for the system of coupled semi-linear heat equations ut=uxx+f(v) and vt=vxx+g(u), for (x,t)∈(0,1)×(0,T), where f and g are given functions. The system has the homogeneous Neumann boundary conditions and the bounded nonnegative initial conditions that are compatible with the boundary conditions. The existence result is established by using the method of upper and lower solutions. We obtain sufficient conditions for finite time quenching of solutions. The quenching set is also provided. From the quenching set, it implies that the quenching solution has asymmetric profile. We prove the blow-up of time-derivatives when quenching occurs. We also find the criteria to identify simultaneous and non-simultaneous quenching of solutions. For non-simultaneous quenching, the corresponding quenching rate of solutions is given.