We extend our previous analysis of the decaysB d,s → µ + µ − in the Minimal Supersymmetric Standard Model (MSSM) to include gluino and neutralino contributions. We provide analytic formulae, valid at large values of tan β, for the scalar and pseudoscalar Wilson coefficients arising from neutral Higgs boson exchange diagrams with gluinos and neutralinos. Together with the remaining contributions (W ± , H ± ,χ ± ), and assuming the Cabibbo-KobayashiMaskawa (CKM) matrix to be the only source of flavour violation, we assess their implications for the branching fractions, since (i) the theoretical errors cancel to a large extent, and (ii) it offers a theoretically clean way of extracting the ratio |V td /V ts | in the Standard Model, which predicts R SM ∼ |V td /V ts | 2 ∼ O(10 −2 ). Exploring three different scenarios of modified minimal flavour violation (MFV), we find that part of the MSSM parameter space can accommodate largeB d,s → µ + µ − branching fractions, while being consistent with various experimental constraints. More importantly, we show that the ratio R can be as large as O(1), while the individual branching fractions may be amenable to detection by ongoing experiments. We conclude that within the MSSM with large tan β the decay rates ofB d → µ + µ − and B s → µ + µ − can be of comparable size even in the case where flavour violation is due solely to the CKM matrix.