% run a simulation to look at the distribution of numbe of mutants that
% will arise, if each cell upon being plated has a certain probability of
% mutating
mut_rate = 1e-5; % probability that a given cell will mutate
n_cells = 128000;
n_simulations = 1000;
outcomes = zeros(n_simulations,1);
for i=1:n_simulations
outcomes(i) = sum(rand(n_cells,1)