Solutions to nonlinear integro-differential systems are regular outside a negligible closed subset whose Hausdorff dimension can be explicitly bounded from above. This subset can be characterized using quantitative, universal energy thresholds for nonlocal excess functionals. The analysis is carried out via the use of nonlinear potentials and allows to derive fine properties of solutions under sharp assumptions on data and kernel coefficients.