I want to compute the magnetic field due to a homogeneously charged, rotating sphere with radius $R$, angular velocity $\vec{\omega} = \omega\hat{z}$ and total charge $Q$. I want to use Biot-Savart law, so I first computed the current density $\vec{j}(r,\theta,\phi) = \frac{\omega Q}{4 \pi R} \delta(R-r) \cdot (-\sin\theta \sin\phi,\sin\theta\cos\phi,0)^T$. After taking out all the constants with respect to integration and removing the integral over the radial component with the $\delta$ function, what stays in the Biot-Savart integral is
$$\vec{B}(\vec{r}) \propto \iint_{\theta,\phi}^{\pi,2\pi}\frac{\vec{j}(\vec{y}) \times (\vec{r} - \vec{y})}{|\vec{r}-\vec{y}|^3}\sin{\theta}\;\text{d}\phi\;\text{d}\theta$$
where $\vec{r} = |\vec{r}|\cdot(\sin\theta_r\cos\phi_r, \sin\theta_r\sin\phi_r,\cos\theta_r)^T$ is the point where the $B$-field is evaluated, $\vec{y} = R\cdot(\sin\theta\cos\phi, \sin\theta\sin\phi,\cos\theta)$ is a point on the sphere. and the $\vec{j}$ is now missing the $\delta$ part, and some constants.
The problem I have is, that the vector inside the integral is
$$\begin{pmatrix} \sin\theta\cos\phi(|\vec{r}|\cos{\theta_r} - R\cos\theta )\\\sin\theta\sin\phi(|\vec{r}|\cos\theta_r-R\cos\theta)\\-|\vec{r}|\sin\theta\sin\phi\sin\theta_r\sin\phi_r-|\vec{r}|\sin\theta\cos\phi\sin\theta_r\cos\phi_r+R\sin^2\theta\end{pmatrix}.$$ Since all terms involving $\sin\phi$ or $\cos\phi$ will integrate to zero, this implies that I will get a $B$-field that is along the $z$ axis, no matter where it's evaluated. This can't be, because in the next part I'll show that far away the field looks like that of a dipole.
I have been sitting over this for hours, can someone help me clear this up? Where am I going wrong?