Next to leading order evolution of SIDIS processes in the forward region

We compute the order $\alpha_s^2$ quark initiated corrections to semi-inclusive deep inelastic scattering extending the approach developed recently for the gluon contributions. With these corrections we complete the order $\alpha_s^2$ QCD description of these processes, verifying explicitly the factorization of collinear singularities. We also obtain the corresponding NLO evolution kernels, relevant for the scale dependence of fracture functions. We compare the non-homogeneous evolution effects driven by these kernels with those obtained at leading order accuracy and discuss their phenomenological implications.