Many models of one-dimensional local random growth are expected to lie in the Kardar-Parisi-Zhang (KPZ) universality class. For such a model, the interface profile at advanced time may be viewed in scaled coordinates specified via characteristic KPZ scaling exponents of one-third and two-thirds. When the long time limit of this scaled interface is taken, it is expected—and proved for a few integrable models— that, up to a parabolic shift, the Airy₂ process $\mathcal{A} : \mathbb{R} \to \mathbb{R}$ is obtained. This process may be embedded via the Robinson-Schensted-Knuth correspondence as the uppermost curve in an $\mathbb{N}$-indexed system of random continuous curves, the Airy line ensemble. Among our principal results is the assertion that the Airy₂ process enjoys a very strong similarity to Brownian motion $B$ (of rate two) on unit-order intervals. This result yields bounds on the Airy₂ probabilities of a large class of events from the counterpart bounds on Brownian motion probabilities. The result has the consequence that the Radon-Nikodym derivative of the law of $\mathcal{A}$ on say $[−1, 1]$, with respect to the law of $B$ on the same interval, lies in every $L^p$ space for $p \in (1, \infty).$ In fact, the quantitative comparison of probability bounds we prove also holds for the scaled energy profile with Dirac delta initial condition of the model of Brownian last passage percolation, a model that lies in the KPZ universality class and in which the energy of paths in a random Brownian environment is maximised. Our technique of proof harnesses a probabilistic resampling or Brownian Gibbs property satisfied by the Airy line ensemble after parabolic shift, and this article develops Brownian Gibbs analysis of this ensemble begun in [CH14] and pursued in [Ham19a]. Our Brownian comparison for scaled interface profiles is an element in the ongoing programme of studying KPZ universality via probabilistic and geometric methods of proof, aided by limited but essential use of integrable inputs. Indeed, the comparison result is a useful tool for studying this universality class. We present and prove several applications, concerning for example the structure of near ground states in Brownian last passage percolation, or Brownian structure in scaled interface profiles that arise from evolution from any element in a very general class of initial data.