Summary: In this paper, we formulate a new Wolbachia infection model in a two-sex mosquito population with stage structure. Some key factors of Wolbachia infection, including cytoplasmic incompatibility (CI), male killing (MK) effect, maternal transmission, fecundity cost due to fitness effect and different mortality rates for infected individuals, are captured. Dynamical analysis has been carried out, and the basic reproduction number \(R_0\) for Wolbachia infection has been calculated. Our analysis shows that Wolbachia can establish in a mosquito population if \(R_0\) is greater than unity. If \(R_0\) is less than unity, Wolbachia establishment still can be achieved if backward bifurcation occurs. Under this circumstance, the initial values lying in the basin of attraction of the stable Wolbachia-established equilibrium are essential to guarantee Wolbachia establishment. In particular, the method to find the basin of attraction and evaluate the threshold initial values is given. Besides, according to a comparison of different releasing strategies, it is shown that, from the perspective of economy and disease control, keeping the number of infected female mosquitoes to a necessary minimum by relying on higher number of male mosquitoes released is a desirable strategy. Moreover, global and local sensitivity analysis and numerical simulation have been performed to explore the impact of model parameters to the success of population establishment. Our results suggest that low levels of MK effect and fitness costs as well as high levels of CI and maternal inheritance are in favor of Wolbachia establishment. Moreover, not considering MK effect and incomplete CI effect may result in the underestimation of the number of infected mosquitoes needed to be released.