What Selberg proved is that for $E\subset \mathbb R$, we have that the limit as $T\to\infty$ of $ \frac{1}{T}\mu(T\le t\le 2T\,|\,\arg(\zeta(1/2+i t)/\sqrt{1/2\log\log t}\in E) $ where $\mu$ is Lebesgue measure, is equal to the integral over $E$ of a Gaussian Random Variable with mean $0$ and standard deviation $1$: $ \frac{1}{\sqrt{2\pi}}\int_E \exp(-x^2/2)\, dx. $
Edit: Here's an example that may help clarify, using the harmonic conjugate $\log|\zeta(1/2+i t)|$ (which is implemented in Mathematica). The analog of Selberg's theorem is true for this function as well. The plot is of $\log|\zeta(1/2+i t)|/\sqrt{1/2\log\log(t))}$, for $50\le t\le 100$. Note it looks nothing like a Gaussian.

Here's a histogram of $50000$ equally spaced values values taken by this function:

Extreme negative values (near the Riemann zeros) are extremely scarce, as are large positive values. The fit to the bell curve is not yet good, but the Riemann zeta function approaches its asymptotic behavior very slowly.