Athurva Gore/LabNotes/2009-8-20

From ZhangLabWiki
Revision as of 22:08, 20 August 2009 by >Ajgore (→‎False Positive SNPs)
Jump to navigation Jump to search

Navigation

False Positive SNPs

  • Currently checking to see if clonal reads could be influencing false positive SNPs.
  • Added a step to variant caller that removes clonal reads (rmdup)
  • Testing on HL003_merged "high false positive" set.
    • Removing clonal reads does not have any effect; downsampling seems to remove them in a similar way
    • No advantage gained.
  • MAQ, however, does something more interesting
    • MAQ only calls 2500 SNPs from HL003 NA12878 set.
    • Most of them are at dbSNP locations (only 300 not)
    • MAQ is being much more conservative somehow
      • Taking into account mapping quality? What about taking into account neighboring consensus quality?