In ref.[1] we presented the JIMWLK Hamiltonian for high energy evolution of QCD amplitudes at the next-to-leading order accuracy in α s . In the present paper we provide details of our original derivation, which was not reported in [1], and provide the Hamiltonian in the form appropriate for action on color singlet as well as color nonsinglet states. The rapidity evolution of the quark dipole generated by this Hamiltonian is computed and compared with the corresponding result of Balitsky and Chirilli [2]. We then establish the equivalence between the NLO JIMWLK Hamiltonian and the NLO version of the Balitsky's hierarchy [3], which includes action on nonsinglet combinations of Wilson lines. Finally, we present complete evolution equation for three-quark Wilson loop operator, thus extending the results of Grabovsky [4].