In Devitt et al. 2013's introduction to quantum error correction, the authors mention (bottom of page 12) how the stabilizer group for $N$ qubits is abelian.
More specifically, here is the quote:
An $N$-qubit stabilizer state $\lvert\psi\rangle_N$ is then defined by the $N$ generators of an Abelian (all elements commute) subgroup $\mathcal G$ of the $N$-qubit Pauli group, $$\mathcal G=\{K^i\,|\,K^i\lvert\psi\rangle=\lvert\psi\rangle,\,[K^i,K^j]=0,\forall (i,j)\}\subset \mathcal P_N.$$
I am confused by this. Is the stabilizer subgroup $\mathcal G$ defined as an abelian subgroup of elements of $\mathcal P_N$ that stabilizes $\lvert\psi\rangle$, or is it instead the case that the subgroup of elements of $\mathcal P_N$ that stabilize $\lvert\psi\rangle$ is abelian?
If the latter, doesn't this introduce ambiguity in the definition? There could be other elements that stabilize $\lvert\psi\rangle$ but not commute with $\mathcal G$.
If the former, how is this shown? I can see why the action of $K^i K^j$ and $K^j K^i$ is identical on $\lvert\psi\rangle$, but how do you show that $K^i=K^j$?