“…For example, assuming s a and s A are very small that higher order terms involving and above are negligible, for HA, ID and MA, The three mechanisms differed in h, s a and s A . For FS (see Nabutanyi & Wittmann (2021) supplementary Section A.3), and for FD, where, h t , h, μ a , μ A , s A , s a , s a,t , and s A,t are as defined under ”Selection Mechanisms” subsection. We adopted the numerical scheme in Zhao et al (2013); Xu et al (2019) to solve Equation (1) subject to the initial condition f ( x, 0) = x 0 and the boundary condition that conserves probability at all times. The scheme basically discretizes both time t and allele frequency x over the interval 0 and 1 inclusive.…”