Two-sided Gaussian bounds for fundamental solutions of non-divergence form parabolic operators with H"older continuous coefficients