  We describe the implementation of infrared subtractions for two-loop QCD
corrections to quark-antiquark annihilation to electroweak final states. The
subtractions are given as form-factor integrands whose integrals are known. The
resulting subtracted amplitudes are amenable to efficient numerical
integration. Our procedure is based on the universality of infrared
singularities and requires a relatively limited set of subtractions, whose
number grows as the number of two-loop diagrams, rather than with the number of
singular regions of integration.
