Thermalization of non-Abelian gauge theories at next-to-leading order