Résumé : New physics in the neutrino sector might be necessary to address anomalies between different neutrino oscillation experiments. Intriguingly, it also offers a possible solution to the discrepant cosmological measurements of $H_0$ and $\sigma_8$. We show here that delaying the onset of neutrino free-streaming until close to the epoch of matter-radiation equality can naturally accommodate a larger value for the Hubble constant $H_0=72.3 \pm 1.4$ km/s/Mpc and a lower value of the matter fluctuations $\sigma_8=0.786\pm 0.020$, while not degrading the fit to the cosmic microwave background (CMB) damping tail. We achieve this by introducing neutrino self-interactions in the presence of a non-vanishing sum of neutrino masses. This strongly interacting neutrino cosmology prefers $N_{\rm eff} = 4.02 \pm 0.29$, which has interesting implications for particle model-building and neutrino oscillation anomalies. We show that the absence of the neutrino free-streaming phase shift on the CMB can be compensated by shifting the value of other cosmological parameters, hence providing an important caveat to the detections made in the literature. Due to their impact on the evolution of the gravitational potential at early times, self-interacting neutrinos and their subsequent decoupling leave a rich structure on the matter power spectrum. In particular, we point out the existence of a novel localized feature appearing on scales entering the horizon at the onset of neutrino free-streaming. While the interacting neutrino cosmology provides a better global fit to current cosmological data, we find that traditional Bayesian analyses penalize the model as compared to the standard cosmological. Our analysis shows that it is possible to find radically different cosmological models that nonetheless provide excellent fits to the data, hence providing an impetus to thoroughly explore alternate cosmological scenarios.