  We consider Feynman integrals with algebraic leading singularities and total
differentials in $\epsilon\,\mathrm{d}\ln$ form. We show for the first time
that it is possible to evaluate integrals with singularities involving
unrationalizable roots in terms of conventional multiple polylogarithms, by
either parametric integration or matching the symbol. As our main application,
we evaluate the two-loop master integrals relevant to the $\alpha \alpha_s$
corrections to Drell-Yan lepton pair production at hadron colliders. We
optimize our functional basis to allow for fast and stable numerical
evaluations in the physical region of phase space.
