User Tools

Site Tools


lecture_notes:05-16-2011

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision
Previous revision
Last revision Both sides next revision
lecture_notes:05-16-2011 [2011/05/16 21:45]
eyliaw created
lecture_notes:05-16-2011 [2011/05/16 22:42]
eyliaw
Line 1: Line 1:
 ====== Seqprep ====== ====== Seqprep ======
 +Phred:
 +  * 66 -> 0 (33)  (Or not... makes sense to not have 0)
 +  * 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.
 +
 +
 How does Seqprep do error correction? How does Seqprep do error correction?
  
Line 7: Line 14:
     = -10 * log_10(1-P(Y=x|Xi=x))     = -10 * log_10(1-P(Y=x|Xi=x))
  
-Phred: +Quality of our base call based on an alignment: ​ Q1 (quality of read 1), Q2 (quality of read 2), Qa (alignment quality).
-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: If we assume that the alignment is perfect:
   * Bases 1 and 2 agree: Q1 + Q2   * Bases 1 and 2 agree: Q1 + Q2
   * Bases 1 and 2 disagree: Q1 - Q2?   * Bases 1 and 2 disagree: Q1 - Q2?
-     * (1 - P1)(2/3( 
lecture_notes/05-16-2011.txt · Last modified: 2011/06/13 00:40 by karplus