We study the moduli stabilization from the viewpoint of modular flavor symmetries. We systematically analyze stabilized moduli values in possible configurations of flux compactifications, investigating probabilities of moduli values and showing which moduli values are favorable from our moduli stabilization. Then, we examine their implications on modular symmetric flavor models. It is found that distributions of complex structure modulus τ determining the flavor structure are clustered at a fixed point with the residual ℤ3 symmetry in the SL(2, ℤ) fundamental region. Also, they are clustered at other specific points such as intersecting points between |τ|2 = k/2 and Re τ = 0,±1/4,±1/2, although their probabilities are less than the ℤ3 fixed point. In general, CP-breaking vacua in the complex structure modulus are statistically disfavored in the string landscape. Among CP-breaking vacua, the values Re τ = ±1/4 are most favorable in particular when the axio-dilaton S is stabilized at Re S = ±1/4. That shows a strong correlation between CP phases originated from string moduli.