A main feature of high-energy scattering in QCD is saturation in the number density of gluons. This phenomenon is described by non-linear evolution equations, JIMWLK and BK, which have been derived at leading logarithmic accuracy. In this paper we generalize this framework to include running coupling corrections to the evolution kernel. We develop a dispersive representation of the dressed gluon propagator in the background of Weiszäcker Williams fields and use it to compute O(β n−1 0 α n s ) corrections to the kernel to all orders in perturbation theory. The resummed kernels present infrared-renormalon ambiguities, which are indicative of the form and importance of non-perturbative power corrections. We investigate numerically the effect of the newly computed perturbative corrections as well as the power corrections on the evolution and find that at present energies they are both significant.