- Make chrV vcf and fasta (in helperfiles)
- Read in vcf
- initially two types of mutations m1 are the markers of each haploid genome, m3 are SNPs - at this stage there are 17 alleles - and about 8K SNPs
- pick a subset of m3 SNPs to make causative (NSNPs), these become flavor m2 mutations
- give these mutations causative effects - with effects scaled by 1/ (minorAlleleCount)^scaling_constant - the scaling constant allows one to make rarer alleles have bigger effects
- now mix these 17 founders to create the population with counts given by K0 - known mixing proportions * K0 pop size
- random mate for 12 generations
- grow to size K1 before starting actual evolution experiment