Additional file 5 Simulation results with Evolved genotype Recruitment/transition selection Mutation W M >1 W M =1 W M >1 W M >1 W M <1 Protected pair Region Outcome Male Female W F =1 W F >1 W F >1 W F <1 W F >1 polymorphism a→A/f→f -I R f/f - f/f 698 0 0 0 0 0 II RH f -/f - f/f - 0 0 90 0 1 127 III R m/f - f -/f - 0 0 1850 475 1369 6 a→A/f→f + I RH m/m m/f + 0 0 136 57 2 341 II R m/f + f + /f + 0 0 2016 483 1538 6 a→A/m→m -I R m -/f f/f 4715 0 0 0 0 0 a→A/m→m + I RH m/m m/m + 0 0 112 27 4 220 II R m/m + m + /m + 0 0 103 27 79 1 III RH m + /m + m + /f 0 0 49 0 0 81 IV R m


), but that the borders are not in exactly the same positions across the range of k and h. Due to computational constraints we used a slightly lower threshold of < 10 -11 for the population dynamics (see Methods).
Additional file 5 Simulation results with Evolved genotype Recruitment/transition selection Mutation W M >1 W M =1 W M >1 W M >1 W M <1 Protected pair Region Outcome Male Female W F =1 W F >1 W F >1 W F <1 W F >1 polymorphism a→A/f→f -I R f/f - f/f 698 0 0 0 0 0 II RH f -/f - f/f - 0 0 90 0 1 127 III R m/f - f -/f - 0 0 1850 475 1369 6 a→A/f→f + I RH m/m m/f + 0 0 136 57 2 341 II R m/f + f + /f + 0 0 2016 483 1538 6 a→A/m→m -I R m -/f f/f 4715 0 0 0 0 0 a→A/m→m + I RH m/m m/m + 0 0 112 27 4 220 II R m/m + m + /m + 0 0 103 27 79 1 III RH m + /m + m + /f 0 0 49 0 0 81 IV R m