We demonstrate that ultracold interacting bosonic atoms in an optical lattice with large on-site population show sub-Poisson on-site and intersite atom-number fluctuations. The experimental observations agree with numerical predictions of the truncated Wigner approximation. The correlations persist in the presence of multimode atom dynamics and even over large spatially extended samples involving several sites