[step:Check the only overlap ambiguities using the Jacobi identity]
Because all left sides of the reduction rules have length $2$, there are no inclusion ambiguities. The only possible overlap ambiguity occurs in a word
\begin{align*}
e_k e_j e_i
\end{align*}
with $k>j>i$, where one may first reduce the left pair $e_k e_j$ or the right pair $e_j e_i$.
Let
\begin{align*}
R_{ab}:=e_a e_b-e_b e_a-[e_a,e_b]\in T(\mathfrak g)
\end{align*}
for $a>b$. The ambiguity for $e_k e_j e_i$ is represented by
\begin{align*}
R_{kj}e_i-e_kR_{ji}.
\end{align*}
Expanding gives
\begin{align*}R_{kj}e_i-e_kR_{ji}=-e_j e_k e_i-[e_k,e_j]e_i+e_k e_i e_j+e_k[e_j,e_i].\end{align*}
Reduce the remaining inverted adjacent pairs in the degree-three terms. First, \begin{align*}-e_j e_k e_i\equiv -e_j e_i e_k-e_j[e_k,e_i].\end{align*} modulo the ideal generated by the relations, and then
\begin{align*}-e_j e_i e_k\equiv -e_i e_j e_k-[e_j,e_i]e_k.\end{align*}
Also, \begin{align*}e_k e_i e_j\equiv e_i e_k e_j+[e_k,e_i]e_j.\end{align*} and
\begin{align*}e_i e_k e_j\equiv e_i e_j e_k+e_i[e_k,e_j].\end{align*}
Substituting these reductions, the degree-three ordered terms cancel, and the ambiguity is congruent to
\begin{align*}
-e_j[e_k,e_i]-[e_j,e_i]e_k-[e_k,e_j]e_i+e_i[e_k,e_j]+[e_k,e_i]e_j+e_k[e_j,e_i].
\end{align*}
Using the relation $uv-vu=[u,v]$ for $u,v\in\mathfrak g$, this expression is congruent to
\begin{align*}
-[e_j,[e_k,e_i]]-[[e_j,e_i],e_k]-[[e_k,e_j],e_i].
\end{align*}
By antisymmetry of the Lie bracket, this equals
\begin{align*}
[e_j,[e_i,e_k]]+[e_k,[e_j,e_i]]+[e_i,[e_k,e_j]].
\end{align*}
This is $0$ by the Jacobi identity in $\mathfrak g$. Hence every overlap ambiguity is resolvable.
The standard reduction criterion for a terminating linear rewriting system now applies: if all overlap and inclusion ambiguities are resolvable, then every element of the quotient by the reduction relations has a unique normal form as a $k$-linear combination of irreducible words. In the present system, the quotient by the reduction relations is exactly $U(\mathfrak g)$, because the reductions are precisely the defining relations $e_j e_i-e_i e_j-[e_j,e_i]$ for $j>i$.
[/step]