This is an old revision of the document!
How does Seqprep do error correction?
Quality scores Q_i = -10 * log_10(P(i is wrong)) Given that base Xi is x:
= -10 * log_10(1-P(Y=x|Xi=x))
Phred: 66 → 0 (33) 33 - 64 = -31 in most cases
Systematic base-calling biases in Illumina exist. I.e. read one thing going one way, something else going the other.
Quality of an alignment: Q1 (quality of read 1), Q2 (quality of read 2), Qa (alignment quality: agree and disagree).
If we assume that the alignment is perfect: