where generated from a quasispecies population Q as follows. First, a quasispecies sequence q Q is randomly chosen accordingly to its unknown frequency f q . A read starting position is generated from the uniform distribution and then a read r is produced from quasispecies q with j sequencing errors. The probability of this event is calculated as h q r j l lj j , ( )