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.