It occurred to me that since the population genetics and evolutionary biology fields are obsessed with Kimura’s substitution formula to the point of literal unreason, instead of trying to show them how Kimura made an algebraic mistake and why the formula only applies to one specific case instead of everything, it would be much more useful to demonstrate how, with a few modifications, Kimura’s equation could serve as the foundation of a predictive calculator that is considerably more accurate and useful than the original equation.
Kimura’s Fixation Calculator: Providing Neutral Theory With Predictive Capacity
Neutral theory has stood for fifty-seven years on a simple result: the substitution rate k equals the per-site mutation rate μ. This identity, derived by Kimura in three lines, rests on canceling two quantities that share a letter but not a meaning: the census number of breeding adults N (which supplies mutations) and the variance effective population size Nₑ (which governs drift and fixation). The cancellation in the derivation is valid in the special case of asexual bacteria where N ≈ Nₑ. It does not hold in sexually reproducing species, where Nₑ/N is typically ~0.1 (Frankham 1995).
Rejecting the incorrect application of the derivation and treating the realized substitution rate as the minimum of three serial constraints—input flux, polymorphism throughput, and selection cost—yields Kimura’s Fixation Calculator. The selection-cost term is a simple expression in four independently measurable parameters (maximum reproductive differential s_max ≈ 1, Selective Turnover Coefficient d, genome length L, and effective population size Nₑ). The full calculator recovers k ≈ μ for bacteria while predicting the observed compression of rates across sexual eukaryotes, where the selection term sets a ceiling two to five orders of magnitude below textbook expectations based on the standard derivation.
Validated on fourteen sexual species pairs plus the E. coli LTEE (all calibrations independent of molecular clocks), the calculator provides forward prediction of k from organismal parameters, inverse inference of divergence time or Nₑ from observed substitutions, and joint constraint surfaces. Where the textbook supplies a single number, the calculator returns a mechanistically grounded range consistent with observable biological reality.
You can read the whole paper if you are a serious glutton for punishment or if you want to understand why no less than nine scientific fields will be seeing significant future adjustments. This paper will be one of the new appendices in the second edition of Probability Zero, since there really is no need for the Sakana study and the rejection of the MITTENS paper means that there is no reason to add it at the back as well.

