Abstract We extend the quark mean-field (QMF) model for nuclear matter and study the possible presence of quark matter inside the cores of neutron stars. A sharp first-order hadron-quark phase transition is implemented combining the QMF for the hadronic phase with “constant-speed-of-sound” parameterization for the high-density quark phase. The interplay of the nuclear symmetry energy slope parameter, L, and the dimensionless phase transition parameters (the transition density n trans/n 0, the transition strength Δε/ε trans, and the sound speed squared in quark matter ) are then systematically explored for the hybrid star properties, especially the maximum mass M max and the radius and the tidal deformability of a typical 1.4 M ⊙ star. We show the strong correlation between the symmetry energy slope L and the typical stellar radius R 1.4, similar to that previously found for neutron stars without a phase transition. With the inclusion of phase transition, we obtain robust limits on the maximum mass (M max < 3.6 M ⊙) and the radius of 1.4 M ⊙ stars (R 1.4 ≳ 9.6 km), and we find that a phase transition that is too weak (Δε/ε trans ≲ 0.2) taking place at low densities ≲1.3–1.5 n 0 is strongly disfavored. We also demonstrate that future measurements of the radius and tidal deformability of ∼1.4 M ⊙ stars, as well as the mass measurement of very massive pulsars, can help reveal the presence and amount of quark matter in compact objects.