Summary: This paper is concerned with the existence of wavefront solutions to a system of degenerate quasilinear reaction-diffusion equations of mixed quasi-monotone properties in the form \[\partial u_i /\partial t = \nabla \cdot \left(D_i \left(u_i\right) \nabla u_i\right) + f_i \left(\mathbf{u}\right) - \infty < x < \infty, \quad t > 0\] for \(i = 1, \ldots, n\). The important features of this system are that some of the diffusion coefficients \(D_i \left(u_i\right)\) are density dependent and may vanish at certain value of \(u_i\). and that each function \(f_i\) is quasi-monotone increasing for some components of \(\mathbf{u =} \left(u_1, \ldots, u_n\right)\) and decreasing for other components of \(\mathbf{u}\). Such systems model reaction-diffusion processes with density driven diffusion mechanism. Under certain general conditions we prove the existence of a traveling wave solution that is between a pair of coupled upper and lower solutions. A predator-prey model with nonlinear diffusion is used as an illustration of application. The presence of wavefront solutions flowing toward the coexistence states is established by constructing appropriate upper and lower solutions.