I've been staring at Zassenhaus's table, and it looks to me like there is a uniform way to describe all the exceptional near fields. Let $\Gamma$ be a finite subgroup of the quaternions, with entries in number field $K$. Let $\pi$ be a prime of $\mathcal{O}_K$ which does not appear in the denominator of any of the entries of $\Gamma$ and let $p$ be the norm of $\pi$. Assume $p$ is odd. Then we can reduce the entries of $\Gamma$ modulo $\pi$ and get a subgroup of $\mathbb{F}_p\langle i,j \rangle / (i^2=j^2=-1,\ ij=-ji)$. This ring is known to be isomorphic to $\mathrm{Mat}_{2 \times 2}(\mathbb{F}_p)$, so we get an action of $\Gamma$ on $\mathbb{F}_p^2$. In all of the examples that occur, we can make this more explicit: In each case $\pi$ is principal and we can choose a generator $\alpha$ which is a sum $a^2+b^2+c^2+d^2$ of four squares in $K$. So we can realize the vector space $\mathbb{F}_p^2$ explicitly as the quotient of the quaternions by the right ideal generated by $g:=a+bi+cj+dk$.
All of Zassenhaus's examples are of the form $\Gamma \times C$ acting on $\mathbb{F}_p^2$, where $C$ is a subgroup of $\mathbb{F}_p^{\ast}$ acting by scalars and chosen such that $p^2-1 = |\Gamma| \times |C|$. In terms of the original question, we then take $n=p^2$ and $G = (|\Gamma| \times |C|) \ltimes (\mathbb{Z}/p)^2$.
The tetrahedral cases Take $\Gamma$ to be the double cover of the rotational symmetries of the tetrahedron. This can be embedded in the quaternions with rational coefficients, as the group $\{ \pm 1, \pm i, \pm j, \pm k, (\pm 1 \pm i \pm j \pm k)/2 \}$. So $|\Gamma|=24$.
We can take $p=5$, $g=2+i$, $|C|=1$ or $p=11$, $g=3+i+j$, $|C|=5$.
The octahedral cases Take $\Gamma$ to be the double cover of the rotational symmetries of the octahedron. This can be embedded in the quaternions with coefficients in $\mathbb{Q}(\sqrt{2})$. So $|\Gamma| = 48$.
We can take $p=7$, $\pi = (3 + \sqrt{2})$ and $g=(1+\sqrt{2}/2) + i+j/2+k/2$ and $|C|=1$. We can also take $p=23$, $\pi = (5+\sqrt{2})$ and $g=(1+\sqrt{2}/2) + 3i/2+j+k/2$ and $|C|=11$.
The icosehedral cases Take $\Gamma$ to be the double cover of the rotational symmetries of the icosahedron. This can be embedded in the quaternions with coefficients in $\mathbb{Q}(\sqrt{5})$. So $|\Gamma| = 120$.
We can take $p=11$, $\pi = (4+\sqrt{5})$ and $g = (1 + \sqrt{5})/2+(1 + \sqrt{5})/2 i + j$ and $|C|=1$. We can also take $p=29$, $\pi = (7+2 \sqrt{5})$, $g=(1+\sqrt{5}) + i$ and $|C|=7$. Finally, we can take $p=59$, $\pi = (8+\sqrt{5})$, $g=(1 + \sqrt{5})/2+(1 + \sqrt{5})/2 i + 2j+k$ and $|C|=29$.