We consider the generic divergence form second order parabolic equation with coefficients that are regular in the spatial variables and just measurable in time. We show that the spatial derivatives of its fundamental solution admit upper bounds that agree with the Aronson type estimate and only depend on the ellipticity constants of the equation and the L $infty$ norm of the spatial derivatives of its coefficients. We also study the corresponding stochastic partial differential equations and prove that under natural assumptions on the noise the equation admits a mild solution, given by anticipating stochastic integration.