Interior regularity for weak solutions of ultraparabolic equations in divergence form with discontinuous coefficients