A solution that is simpler still, requiring no differentiation or fancy identities other than one of the Bianchi identities, is as follows.
We know that the Riemann tensor can measure how much the covariant derivatives commute with each other, e.g
\begin{equation}
[ \nabla _{\mu } ,\nabla _{\nu }] k_{\alpha } =-R^{\beta }{}_{\alpha \mu \nu } k_{\beta }
\end{equation}
We can shuffle the indices in (4) and use the antisymmetry of the Lie bracket to get the three equations
\begin{gather}
[ \nabla _{\nu } ,\nabla _{\mu }] k_{\alpha } =R^{\beta }{}_{\alpha \mu \nu } k_{\beta } \tag{5A}\\
[ \nabla _{\mu } ,\nabla _{\alpha }] k_{\nu } =R^{\beta }{}_{\nu \alpha \mu } k_{\beta } \tag{5B}\\
[ \nabla _{\alpha } ,\nabla _{\nu }] k_{\mu } =R^{\beta }{}_{\mu \nu \alpha } k_{\beta }
\end{gather}
We now look at $\displaystyle ( 5\mathrm{C}) -( 5\mathrm{B}) -( 5\mathrm{A})$:
\begin{equation*}
[ \nabla _{\alpha } ,\nabla _{\nu }] k_{\mu } -[ \nabla _{\mu } ,\nabla _{\alpha }] k_{\nu } -[ \nabla _{\nu } ,\nabla _{\mu }] k_{\alpha } =\left( R^{\beta }{}_{\mu \nu \alpha } -R^{\beta }{}_{\nu \alpha \mu } -R^{\beta }{}_{\alpha \mu \nu }\right) k_{\beta }
\end{equation*}
Applying the first Bianchi identity,
\begin{equation*}
R^{\beta }{}_{\alpha \mu \nu } +R^{\beta }{}_{\mu \nu \alpha } +R^{\beta }{}_{\nu \alpha \mu } =0
\end{equation*}
This simplifies to
\begin{equation*}
[ \nabla _{\alpha } ,\nabla _{\nu }] k_{\mu } -[ \nabla _{\mu } ,\nabla _{\alpha }] k_{\nu } -[ \nabla _{\nu } ,\nabla _{\mu }] k_{\alpha } =2R^{\beta }{}_{\mu \nu \alpha } k_{\beta }
\end{equation*}
Expanding the Lie brackets,
\begin{gather*}
[ \nabla _{\alpha } ,\nabla _{\nu }] k_{\mu } -[ \nabla _{\mu } ,\nabla _{\alpha }] k_{\nu } -[ \nabla _{\nu } ,\nabla _{\mu }] k_{\alpha }\\
=( \nabla _{\alpha } \nabla _{\nu } -\nabla _{\nu } \nabla _{\alpha }) k_{\mu } -( \nabla _{\mu } \nabla _{\alpha } -\nabla _{\alpha } \nabla _{\mu }) k_{\nu } -( \nabla _{\nu } \nabla _{\mu } -\nabla _{\mu } \nabla _{\nu }) k_{\alpha }\\
=\nabla _{\alpha } \nabla _{\nu } k_{\mu } -\nabla _{\nu } \nabla _{\alpha } k_{\mu } +\nabla _{\alpha } \nabla _{\mu } k_{\nu } -\nabla _{\mu } \nabla _{\alpha } k_{\nu } +\nabla _{\mu } \nabla _{\nu } k_{\alpha } -\nabla _{\nu } \nabla _{\mu } k_{\alpha }
\end{gather*}
Using the linearity of $\displaystyle \nabla $,
\begin{gather*}
[ \nabla _{\alpha } ,\nabla _{\nu }] k_{\mu } -[ \nabla _{\mu } ,\nabla _{\alpha }] k_{\nu } -[ \nabla _{\nu } ,\nabla _{\mu }] k_{\alpha }\\
=\nabla _{\alpha }( \nabla _{\nu } k_{\mu } +\nabla _{\mu } k_{\nu }) +\nabla _{\mu }( \nabla _{\nu } k_{\alpha } -\nabla _{\alpha } k_{\nu }) -\nabla _{\nu }( \nabla _{\alpha } k_{\mu } +\nabla _{\mu } k_{\alpha })
\end{gather*}
Using Killing's equation $\displaystyle \nabla _{( \rho } k_{\sigma )} =0$ the first and third terms vanish, leaving
\begin{equation*}
[ \nabla _{\alpha } ,\nabla _{\nu }] k_{\mu } -[ \nabla _{\mu } ,\nabla _{\alpha }] k_{\nu } -[ \nabla _{\nu } ,\nabla _{\mu }] k_{\alpha } =\nabla _{\mu }( \nabla _{\nu } k_{\alpha } -\nabla _{\alpha } k_{\nu }) =2\nabla _{\mu } \nabla _{\nu } k_{\alpha }
\end{equation*}
Hence,
\begin{equation}
\boxed{\nabla _{\mu } \nabla _{\nu } k_{\alpha } =R^{\beta }{}_{\mu \nu \alpha } k_{\beta }}
\end{equation}