The problem is actually to verify the divergence theorem by computing both $\iiint_E \text{div } \mathbf{F\space} dV$, which was relatively easy to compute and gives $\frac{256\pi}{3}$. To find $\iint_S \mathbf{F}\cdot d\mathbf{S}$, I parametrized the surface with spherical coordinates: $\mathbf{r}=\langle 4\sin\phi\cos\theta,4\sin\phi\sin\theta,4\cos\phi\rangle$ with $0\leq\phi\leq\pi,0\leq\theta\leq2\pi$. Now, noting that $\iint_S\mathbf{F}\cdot d\mathbf{S}=\iint_S \mathbf{F}\cdot \mathbf{n}\space dS$ where $\mathbf{n}$ is the normal vector to the $S$, and since $S$ is a sphere, we have $\mathbf{n}=\langle\sin\phi\cos\theta,\sin\phi\sin\theta,\cos\phi\rangle$, and $$\mathbf{F\cdot n}=4\sin\phi\cos\phi\cos\theta+4\sin^2\phi\sin\theta\cos\theta+4\sin\phi\cos\phi\cos\theta$$
But the integral of this is 0.
Where am I going wrong?