The significance of the Gallium Anomaly, from the BEST, GALLEX, and SAGE radioactive source experiments, is quantified using different theoretical calculations of the neutrino detection cross section, and its explanation due to neutrino oscillations is compared with the bounds from the analyses of reactor rate and spectral ratio data, β-decay data, and solar neutrino data. In the 3+1 active-sterile neutrino mixing scheme, the Gallium Anomaly is in strong tension with the individual and combined bounds of these data sets. In the combined scenario with all available data, the parameter goodness of fit is below 0.042%, corresponding to a severe tension of 4–5σ, or stronger. Therefore, we conclude that one should pursue other possible solutions beyond short-baseline oscillations for the Gallium Anomaly. We also present a new global fit of νe and $$ {\overline{\nu}}_e $$
ν
¯
e
disappearance data, showing that there is a 2.6–3.3σ preference in favor of short-baseline oscillations, which is driven by an updated analysis of reactor spectral ratio data.