$W^{1,p}$ regularity of solutions to Kolmogorov equation with Gilbarg-Serrin matrix

In ${\mathbb{R}}^d$, $d\ge 3$, consider the divergence and the non-divergence form operators \begin{equation} \tag{$i$} -\Delta - \nabla \cdot (a-I) \cdot \nabla + b \cdot \nabla, \end{equation} \begin{equation} \tag{$ii$} - \Delta - (a-I) \cdot \nabla^2 + b \cdot \nabla, \end{equation} where the second order perturbations are given by the matrix $$a-I=c\mathrm{\mid }x{\mathrm{\mid }}^{-2}x\otimes x,c>-1.$$ The vector field $b:{\mathbb{R}}^d\to {\mathbb{R}}^d$ is form-bounded with the form-bound $\delta >0$ (this includes a sub-critical class $[L^d + L^\infty]^d$, as well as vector fields having critical-order singularities). We characterize quantitative dependence on $c$ and $\delta$ of the $L^q \rightarrow W^{1,qd/(d-2)}$ regularity of the resolvents of the operator realizations of ($i$), ($ii$) in $L^q$, $q\ge 2\vee (d-2)$ as (minus) generators of positivity preserving $L^{\mathrm{\infty }}$ contraction $C_0$ semigroups.