In this talk CP violation in the supersymmetric models, and especially in B-decays is discussed. We review our analysis of the supersymmetric contributions to the mixing CP asymmetries of $B\to \phi K_S$ and $B\to \eta^{\prime} K_S$ processes. Both gluino and chargino exchanges are considered in a model independent way by using the mass insertion approximation method. The QCD factorization method is used, and parametrization of this method in terms of Wilson coefficients is presented in both decay modes. Correlations between the CP asymmetries of these processes and the direct CP asymmetry in $b\to s \gamma$ decay are shown.
