  McMule is a framework for fully differential higher-order QED calculations of
scattering and decay processes involving leptons. It keeps finite lepton
masses, which regularises collinear singularities. Soft singularities are
treated with dimensional regularisation and using FKS$^\ell$ subtraction. We
describe the implementation of the framework in Fortran 95, list the processes
that are currently implemented, and give instructions on how to run the code.
In addition, we present new phenomenological results for muon-electron
scattering and lepton-proton scattering, including the dominant NNLO
corrections. While the applications presented focus on MUonE, MUSE, and P2, the
code can be used for a large number of planned and running experiments.
