We derive corrections to atomic energy levels from disformal couplings in Galileon theories. Through Bayesian inference, we constrain the cutoff radii and Galileon scale via these corrections. To connect different atomic systems, we assume the various cutoff radii related by a one-parameter family of solutions. This introduces a new parameter α which is also constrained. In this model, we predict shifts to muonic helium of δE He 3 = 1.97