We would like to note that even though our alignment model does assume independence of sites, we model a local dependency in our evolutionary model by conditioning our emission probabilities on the nucleotide context in the ancestral sequence. An undoubtable improvement would be to model the dependency as continuous throughout the evolutionary process [17] . However, as noted by the authors themselves, the elaborate MCMC method developed in order to do this makes it a computationally intractable option.