The emergence of resistant strains is a complex phenomenon that requires the study at different scales, even with the need to consider stochasticity. Mathematical models describing this emergence are scarce and with poor predictive capabilities to be used for optimization (to for example find the best treatments to minimize resistance). In this project the aim is to develop minimal models able to predict main behaviors at different time and space scales.