The Color Glass Condensate (CGC) effective field theory (EFT) at leading order describes well the Deep Inelastic Scattering (DIS) inclusive cross section data at small-x as measured by the HERA experiments [1-3]. Recently the inclusive DIS impact factors have been calculated in Next-to-Leading Order (NLO) accuracy in CGC EFT [4-6], and the soft gluon divergence present at NLO has been factorized succesfully .
In this talk we discuss our recent work  on the first comparisons of the NLO DIS cross sections to HERA data. Fitting the HERA reduced cross section data fixes the the non-perturbative initial condition to the BK evolution. Since the available NLO DIS inclusive cross sections are calculated in the massless quark limit, we construct and fit a dataset of light quark only cross sections using an independent parametrization of HERA total and heavy quark data. We find an excellent description of the HERA data. As the NLO BK is computationally expensive , we compare a number of beyond LO prescriptions that approximate the full NLO BK, including the recent evolution parametrized in target rapidity . These beyond LO evolution equations include important higher order contributions by resumming corrections enhanced by large transverse logarithms.
 T. Lappi, H. Mäntysaari, Phys.Rev. D88 (2013) 114020
 E. Iancu, J.D. Madrigal, A.H. Mueller, G. Soyez, D.N. Triantafyllopoulos, Phys.Lett. B750 (2015) 643-652
 B. Ducloué, E. Iancu, G. Soyez, D.N. Triantafyllopoulos, arXiv:1912.09196 [hep-ph]
 G. Beuf, Phys.Rev. D94 (2016) no.5, 054016
 G. Beuf, Phys.Rev. D96 (2017) no.7, 074033
 H. Hänninen, T. Lappi, R. Paatelainen, arXiv:1711.08207 [hep-ph]
 B. Ducloué, H. Hänninen, T. Lappi, Y. Zhu, Phys.Rev. D96 (2017) no.9, 094017
 G. Beuf, H. Hänninen, T. Lappi, H. Mäntysaari, arXiv:2007.01645 [hep-ph]
 T. Lappi, H. Mäntysaari, Phys.Rev. D91 (2015) no.7, 074016
 B. Ducloué, E. Iancu, A.H. Mueller, G. Soyez, D.N. Triantafyllopoulos, JHEP 1904 (2019) 081