Precise predictions of hard scattering processes at the LHC
usually require the inclusion of NNLO corrections in QCD and NLO in
EW. For some specific cases, this might not be sufficient to match
the experimental accuracy, and higher-order radiative contributions
must be taken into account too. In such situations, effects due the
interplay among strong and electroweak interactions may be
relevant. The Drell-Yan process provides a promiment example of
ultra precision physics at the LHC. We recently completed the first
calculation of the mixed QCD-EW corrections to the neutral current
Drell-Yan process for massive leptons. In this proceeding, we
elaborate on the comparison of the exact two-loop virtual amplitude
with its approximation by a pole expansion and show phenomelogical
results in the region of high invariant masses of the lepton pair
system.