  With a interpolation method on the P-$\mu$ plane, a hybrid equation of state
is explored. The quark phase is described by our newly developed
self-consistent two-flavor Nambu$-$Jona-Lasinio model. It retains the
contribution from the vector channel in the Fierz-transformed Lagrangian by
introducing a weighting parameter $\alpha$ [Chin. Phys. C \textbf{43}, 084102
(2019)]. In the hadron phase we use the relativistic mean-field theory. We
study the dependence of hybrid EOS and mass-radius relation on $\alpha$. It is
found that increasing $\alpha$ makes the hybrid EOS softer in the medium
pressure. We can get stellar mass larger than $2M_\odot$. Further, we calculate
the tidal deformability $\tilde\Lambda$ for binary stars and compare with
recent analysis GW170817 [Phys. Rev. X \textbf{9}, 011001 (2019)].
