In this paper, we consider the one-dimensional stochastic heat equation driven by a space time white noise. In two different scenarios: {\it (i)} initial condition $u_0=1$ and general nonlinear coefficient $σ$ and {\it (ii)}: initial condition $u_0=δ_0$ and $σ(x)=x$ (Parabolic Anderson Model), we establish rates of convergence for the uniform distance between the density of (renormalized) spatial averages and the standard normal density. These results are based on the combination of Stein method for normal approximations and Malliavin calculus techniques. A key ingredient in Case (i) is a new estimate on the $L^p$-norm of the second Malliavin derivative.