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.