We study the Higgs boson decays $h^0 \to c \bar{c}, b \bar{b}, b \bar{s}, \gamma \gamma,
g g$ in the Minimal Supersymmetric Standard Model (MSSM) with general quark flavor violation (QFV),
identifying the $h^0$ as the Higgs boson with a mass of 125 GeV. We compute the widths of the
$h^0$ decays to $c \bar c, b \bar b, b \bar s (s \bar b)$ at full one-loop level.
For the $h^0$ decays to photon photon and gluon gluon we compute the widths at NLO QCD level.
We perform a systematic MSSM parameter scan respecting all the relevant constraints, i.e.
theoretical constraints from vacuum stability conditions and experimental constraints,
such as those from K- and B-meson data and electroweak precision data, as
well as recent limits on Supersymmetric (SUSY) particle masses and the
125 GeV Higgs boson data from LHC experiments. From the parameter scan, we find
that the deviations of these MSSM widths from the Standard Model (SM) values can be
quite sizable. All of these large deviations in the $h^0$ decays are due to large
scharm-stop mixing, large scharm/stop involved trilinear couplings
$T_{U23}, T_{U32}, T_{U33}$, large sstrange-sbottom mixing, and large
sstrange/sbottom involved trilinear couplings $T_{D23}, T_{D32}, T_{D33}$.
International Linear Collider (ILC) can observe such large deviations from
the SM at high signal significance. In case the deviation pattern shown here
is really observed at ILC, then it would strongly suggest the discovery of
QFV SUSY (the MSSM with QFV).