Mathematical model and simulations: list of parameters.

*We assume that nonsense and frameshift mutations render the integrase non-functional if it occurs within the first 500 bp of the integrase.
**Gene cassette acquisition rates were conservatively set at least two orders of magnitude lower than low conduction frequencies as reported for example in [59].
***Hill- coefficient arbitrarily set to 1 [within the range presented in [49]] and decay rates were adjusted to ensure efficient antibiotic selection. Values in brackets for Vx, θ, π, and λ are the widened parameter ranges applied in model predictions.