We apply Noether's theorem to observations of main-sequence stars from the Gaia Data Release 2 archive to probe the matter distribution function of the Galaxy. That is, we examine the axial symmetry of stars at vertical heights $z$, $0.2 \le |z| \le 3$ kpc, to probe the quality of the angular momentum $L_z$ as an integral of motion. The failure of this symmetry test would speak to a Milky Way, in both its visible and dark matter, that is not isolated and/or not in steady state. The left-right symmetry-breaking pattern we have observed, north and south, reveals both effects, with a measured deviation from symmetry of typically 0.5%. We show that a prolate form of the gravitational distortion of the Milky Way by the Large Magellanic Cloud, determined from fits to the Orphan stream by Erkal et al., 2019, is compatible with the size and sign of the axial-symmetry-breaking effects we have discovered in our sample of up to 14.4 million main-sequence stars, speaking to a distortion of an emergent, rather than static, nature.

We study a set of solar neighborhood ($d < 3$ kpc) stars from Gaia Data Release 2 to determine azimuthal star count differences, i.e., left and right of the line from the Galactic center through the sun - and compare these differences north and south. In this companion paper to Gardner et al. (2020), we delineate our procedures to remove false asymmetries from sampling effects, incompleteness, and/or interloper populations, as this is crucial to tests of axisymmetry. Particularly, we have taken care to make appropriate selections of magnitude, color, in-plane Galactocentric radius and Galactic $|b|$ and $|z|$. We find that requiring parallax determinations of high precision induces sampling biases, so that we eschew such requirements and exclude, e.g., regions around the lines of sight to the Magellanic clouds, along with their mirror-image lines of sight, to ensure well-matched data sets. After making conservative cuts, we demonstrate the existence of azimuthal asymmetries, and find differences in those, north and south. These asymmetries give key insights into the nature and origins of the perturbations on Galactic matter, allowing us to assess the relative influence of the Magellanic Clouds (LMC & SMC), the Galactic bar, and other masses on the Galactic mass distribution, as described in Gardner et al. (2020). The asymmetry's radial dependence reveals variations that we attribute to the Galactic bar, and it changes sign at a radius of $(0.95 \pm 0.03) R_0$, with $R_0$ the Sun-Galactic-Center (GC) distance, to give us the first direct assessment of the outer Lindblad resonant radius.