  Motivated by tensions between experimental measurements and SM predictions in
$b\to s \ell^+\ell^-$ transitions, we present the first study of non-minimal
flavour-violating Minimal Supersymmetric Standard Model (MSSM) scenarios
contributing to the relevant Wilson coefficients to address the observed
anomalies using SuperIso and MARTY. We calculate the full one-loop analytical
contributions of the general MSSM to Wilson coefficients relevant for flavour
anomalies, together with the anomalous muon magnetic dipole moment $(g-2)_\mu$.
We show that, after imposing theoretical constraints on the flavour-violating
parameters we can find scenarios in agreement with the experimental
measurements that can address at the same time the tensions in flavour
observables and in $(g-2)_\mu$.
