I have calculated the average over initial and sum over final states of the squared amplitude for both Compton scattering $e^-\gamma \rightarrow e^-\gamma$ (QED) and quark-gluon scattering $qg \rightarrow qg$ (QCD).
Both these quantities are in agreement with the literature.
My understanding is that QED is the Abelian limit of QCD. Hence, I expect that to recover the QED analogue of a QCD scattering amplitude, I should be able to send the number of colours $N_C \rightarrow 1$.
However, I do not get that my $qg \rightarrow qg$ scattering amplitude reduces to $e^-\gamma \rightarrow e^-\gamma$ in the $N_C \rightarrow 1$ limit.
Is my understanding incorrect, or is something else going on here that I may naively be overlooking?