  We present the full $\mathcal{O}(\alpha_s)$ supersymmetric QCD corrections
for stop-anti-stop annihilation into electroweak final states within the
Minimal Supersymmetric Standard Model (MSSM). We also incorporate Coulomb
corrections due to gluon exchange between the incoming stops. Numerical results
for the annihilation cross sections and the predicted neutralino relic density
are presented. We show that the impact of the radiative corrections on the
cosmologically preferred region of the parameter space can become larger than
the current experimental uncertainty, shifting the relic bands within the
considered regions of the parameter space by up to a few tens of GeV.
