  In this paper, we prove the Lorentz space $L^{q,p}$-estimates for gradients
of very weak solutions to the linear parabolic equations with
$\mathbf{A}_q$-weights $$u_t-\operatorname{div}(A(x,t)\nabla
u)=\operatorname{div}(F),$$ in a bounded domain $\Omega\times
(0,T)\subset\mathbb{R}^{N+1}$, where $A$ has a small mean oscillation, and
$\Omega$ is a Lipchistz domain with a small Lipschitz constant.
