Recent results from the CDF Collaboration favor a large CP asymmetry in Bs→K-π+ decay, while the standard model prediction is very small. Moreover, the measurement of its branching ratio is lower than the standard model prediction based on the QCD factorization. We compute the gluino-mediated supersymmetry contributions to Bs→K(*) -π+, K(*)-ρ+ decays in the framework of the mass insertion method, and find that for mg2 mq2≤2, the theoretical predictions including the LR and RL mass insertion contributions are compatible with the measurements of Bs→K -π+ decay and B0-B̄0 mixing within 2σ ranges. Using the constrained LR and RL mass insertion parameter spaces, we explore the supersymmetry mass insertion effects on the branching ratios, the direct CP asymmetries, and the polarization fractions in Bs→K*-π+, K -ρ+, K*-ρ+ decays. We find the constrained LR and RL insertions can provide sizable contributions to the branching ratios of Bs→K* -π+, K(*)-ρ+ as well as the direct CP asymmetry and the longitudinal polarization of Bs→K*-ρ+ decay without conflict with all related data within 2σ ranges. Near future experiments at Fermi Lab and CERN LHC-b can test our predictions and shrink/reveal the mass insertion parameter spaces. © 2010 The American Physical Society.