Source:
Chapter 11 of Freedman and Van Proeyen’s supergravity [textbook][1]
An infinitesimal symmetry transformation is determined by
a parameter, call it $\epsilon^A$, and
an operation, call it $\delta_\epsilon$.
The operation $\delta_\epsilon$
depends linearly on the parameter $\epsilon^A$, and
acts on fields, i.e. $\delta_\epsilon(\phi^i)$.
For some global symmetry, $\epsilon^A$ does not depend on the spacetime $x^\mu$.
Another way to say "$\delta_\epsilon$ depends linearly on the parameter $\epsilon^A$,'' is to write $$\delta_\epsilon =\epsilon^A T_A$$
where the $T_A$ are some operations on fields.
Let $\{(t_A)^i{}_j\}$ be the matrix generators of a representation of some Lie algebra.
This Lie algebra (LA) is defined by $[t_A,t_B]=f_{AB}{}^C t_C$.
The action of $T_A$ on the fields is defined with the LA basis elements,
$$ T_A(\phi^i)=-(t_A)^i{}_j \phi^i $$
So then we have
\begin{eqnarray*} \delta_\epsilon(\phi^i) &=& \epsilon^A T_A(\phi^i) \\ &=& -\epsilon^A (t_A)^i{}_j (\phi^j) \end{eqnarray*}
Then the product of two symmetry transformations reads,
\begin{eqnarray*} \delta_{\epsilon_1}\delta_{\epsilon_2}(\phi^i) &=& \epsilon_1{}^A T_A(\epsilon_2{}^B T_B\phi^i) \\ &=& \epsilon_1{}^A T_A(-\epsilon_2{}^B (t_B)^i{}_j \phi^j) \\ &=& -\epsilon_1{}^A \epsilon_2{}^B (t_B)^i{}_j T_A \phi^j \\ &=& -\epsilon_1{}^A \epsilon_2{}^B (t_B)^i{}_j (-(t_A)^j{}_k \phi^k) \\ &=& \epsilon_1{}^A \epsilon_2{}^B (t_B)^i{}_j (t_A)^j{}_k \phi^k \\ \end{eqnarray*}
The authors then go on to state the commutator, which I am concerned with,
$$ \begin{aligned} \ [\delta(\epsilon_1), \delta(\epsilon_2)] \phi^i & = \epsilon_2^B \epsilon_1^A [T_A, T_B] \phi^i \\ & = \epsilon_2^B \epsilon_1^A ([t_B, t_A]\phi)^i = -\epsilon_2^B \epsilon_1^A {f_{AB}}^C (t_C\phi)^i\\ & = \epsilon_2^B \epsilon_1^A {f_{AB}}^C T_C\phi^i \end{aligned} \tag{11.6} $$
Notably,
$\epsilon_1{}^A$ and $\epsilon_2{}^B$ are numbers, and are so are commutative, and
$(t_B)^i{}_j,$ $(t_A)^j{}_k,$ and $\phi^k$ are matrices, and so are associative.
Thus, how is $[\delta_{\epsilon_1},\delta_{\epsilon_2}](\phi^i) = \delta_{\epsilon_1}\delta_{\epsilon_2}(\phi^i) - \delta_{\epsilon_2}\delta_{\epsilon_1}(\phi^i)$
not equal to zero?
As far as I can tell, $\delta_{\epsilon_1}\delta_{\epsilon_2}(\phi^i)=\delta_{\epsilon_2}\delta_{\epsilon_1}(\phi^i).$
Disclaimer:
That formula for the product of symmetry transformations on the fields is my work. The authors have,
$$ \begin{aligned} \delta(\epsilon_1) \delta(\epsilon_2) \phi^i & = \epsilon_1^A T_A \epsilon_2^B [-{(t_B)^i}_j \phi^j] \\ & = \epsilon_1^A \epsilon_2^B {(-t_B)^i}_j T_A \phi^j \\ & = \epsilon_1^A \epsilon_2^B {(-t_B)^i}_j {(-t_A)^j}_k \phi^k . \end{aligned} \tag{11.5} $$
which seems the same as mine but I just want to be cautious.