  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, $\beta$-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$\sigma$, 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 $\nu_e$ and $\bar\nu_e$ disappearance data, showing that there is
a 2.6-3.3$\sigma$ preference in favor of short-baseline oscillations, which is
driven by an updated analysis of reactor spectral ratio data.
