To allow for joint inference we place a Dirichlet prior on the transition matrix at each SNP, centred on the expected values given by Equation (2), and with precision K, 13 14 where vi, j, t = Kρ and ρ is given by Equation (2).