<p>We present the analytic integration of the iterated single-unresolved subtraction terms in the extension of the CoLoRFulNNLO subtraction scheme to color-singlet production in hadron collisions. We exploit the fact that, in this scheme, subtraction terms are defined through momentum mappings which lead to exact phase space convolutions for real emissions. This allows us to write the integrated subtraction terms as parametric integrals, which can be evaluated using standard tools. Finally, we show that the integrated iterated single-unresolved approximate cross section can be written as a convolution of the Born cross section with an appropriately defined insertion operator.</p>