Let's assume $P$ is a principal bundle, $F^A \in \Omega^2(M,Ad(P))$ the curvature 2-form, $Ad(P)$ the adjoint bundle. $d_A$ the covariant differential. For sections in the associated bundle $E=P \times_{(G, \rho)} V$, $d_A$ is just the covariant differential. In local coordinates it is of the form $d_A \rightarrow \partial_{\mu}+\rho_*(A_{\mu})$. $\phi$ is a section in the associated bundle and in local coordinates takes the form $[s(x),\varphi(x)]$ where $s:U \rightarrow P$ is a section in the principal bundle and $\varphi:U \rightarrow V$. The Yang-Mills-Higgs action is \begin{equation} \mathcal{S}_{Y K}: \mathcal{C}(P) \times \Gamma(E) \rightarrow \mathbb{R}, \quad \mathcal{S}_{Y K}[A, \phi]=\int_{M}\left(-\frac{1}{2}\left\langle F^{A}, F^{A}\right\rangle_{\mathrm{Ad}(P)}+\left\langle d_{A} \phi, d_{A} \phi\right\rangle_{E}-m^{2}\langle\phi, \phi\rangle_{E}\right) d \nu_{g} \end{equation} The variation $A\mapsto A+\omega$ gives the equations of motion \begin{equation} \delta_{A} F^{A}=j \end{equation} \begin{equation} \delta_{A} d_{A} \phi + m^{2} \phi=0 \end{equation} with the codifferential $\delta_A$ and $j \in \Omega^{1}(M , \operatorname{Ad}(P))$ implicitly defined by \begin{equation} \langle j, \omega\rangle_{\mathrm{Ad}(P)}=-2 \operatorname{Re}\left(\left\langle d_{A} \phi, \rho_{*}(\omega) \phi\right\rangle_{E}\right)\quad\text{for all }\omega. \end{equation} In physics, the current is defined by \begin{equation} j_{\nu}^{a}=-i\left(\left(D_{\nu} \varphi_{i}\right)^{\dagger}\left(T_{a}^{r} \varphi\right)_{i}-\left(T_{a}^{r} \varphi\right)_{i}^{\dagger} D_{\nu} \varphi^{j}\right) \end{equation} where $T_a$ is a basis of the Lie algebra and $T_a^r=\rho_*(T_a)$, $D_{\nu}=\partial_{\nu}+A_{\nu}^aT_a^r$. $\varphi_i$ is just the $i$-th component of $\varphi$. The $i$'s come into play due to the definition of physicists that every Lie algebra element is multiplied with $I$.
$\mathbf{Question}$: How exactly can one derive the physical local coordinate expression from the mathematical definition?