In this paper, we examine fluctuations of polynomial linear statistics for the Anderson model on $\mathbb{Z}^d$ for any potential with finite moments. We prove that if normalized by the square root of the size of the truncated operator, these fluctuations converge to a Gaussian limit. For a vast majority of potentials and polynomials, we show that the variance of the limiting distribution is strictly positive, and we classify in full the rare cases in which this does not happen.