Thinking about $SO(3)$. Any rotation matrix $R$ can be written $$ R = e^{\theta \hat{n}\cdot J} $$ where $J$ is a vector the three skew-symmetric generators of rotation $J_x$, $J_y$, and $J_z$. In this case $R$ representation a rotation of angle $\theta$ about axis $\hat{n}$. These three $J$ matrices can be written as \begin{align*} J_x =& \begin{pmatrix}0 & 0 & 0 \\ 0 & 0 & -1 \\ 0 & 1 & 0 \end{pmatrix}\\ J_y =& \begin{pmatrix}0 & 0 & 1 \\ 0 & 0 & 0 \\ -1 & 0 & 0 \end{pmatrix}\\ J_z =& \begin{pmatrix}0 & -1 & 0 \\ 1 & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}\\ \end{align*} We have $$ [J_i, J_j] = \epsilon_{ijk} J_k $$ Thinking of the space spanned by the $J_i$ matrices as the Lie algebra $\mathfrak{so}(3)$ of the Lie group $SO(3)$ we can understand that $\epsilon_{ijk}$ are the structure coefficients of the Lie algebra.
However, in this case we ALSO have the property that $$ (J_i)_{kj} = -\epsilon_{ijk} $$ So the matrix elements of the action of the $J_i$ matrices on $\mathbb{R}^3$ are also given by $\epsilon_{ijk}$. We have $$ [J_i, J_j] = -(J_i)_{jk} J_k $$ Is this a coincidence? Or is there some deeper Lie or representation theory reason why this is expected?