The study of random Fourier series, linear combinations of trigonometric functions whose coefficients are independent (in our case Gaussian) random variables with polynomially bounded means and standard deviations, dates back to Norbert Wiener in one of the original constructions of Brownian motion. A geometric generalization -- relevant e.g.\ to Euclidean quantum field theory with an infrared cutoff -- is the study of random Gaussian linear combinations of the eigenfunctions of the Laplace-Beltrami operator on an arbitrary compact Riemannian manifold $(M,g)$, Gaussian noise $Φ$. I will prove that, when our random coefficients are independent Gaussians whose standard deviations obey polynomial asymptotics and whose means obey a corresponding polynomial upper bound, the resultant random $\mathscr{H}^s$-wavefront set $\operatorname{WF}^s(Φ)$ (defined as a subset of the cosphere bundle $\mathbb{S}^*M$) is either almost surely empty or almost surely the entirety of $\mathbb{S}^*M$, depending on $s \in \mathbb{R}$, and we will compute the threshold $s$ and the behavior of the wavefront set at it. Consequently, the random $C^\infty$-wavefront set $\operatorname{WF}(Φ)$ is almost surely the entirety of the cosphere bundle. The method of proof is as follows: using Sazonov's theorem and its converse, it suffices to understand which compositions of microlocal cutoffs and inclusions of $L^2$-based fractional order Sobolev spaces are Hilbert-Schmidt (HS), and the answer follows from general facts about the HS-norms of the elements of the pseudodifferential calculus of Kohn and Nirenberg.