In this paper, we establish $L_p$ estimates and solvability for time fractional divergence form parabolic equations in the whole space when leading coefficients are merely measurable in one spatial variable and locally have small mean oscillations with respect to the other variables. The corresponding results for equations on a half space are also derived.