Ylaine/2009-7-9: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Ylaine
No edit summary
>Ylaine
Line 2: Line 2:
* Wrote 'extractSNPlines.pl' to get pileup lines by line number (instead of grep-ing for position)
* Wrote 'extractSNPlines.pl' to get pileup lines by line number (instead of grep-ing for position)
  ./extractSNPlines.pl notwo.pileup snoOrigLocs.txt > snp.pileup
  ./extractSNPlines.pl notwo.pileup snoOrigLocs.txt > snp.pileup
* Discovered that pileup does not output SNP calls. Modify 'heterozygousSNPfilter.pl' to extract heterozygous SNPs 'cns.final.snp', put call into hash keyed by position. Then iterate through lines of snp.pileup and count occurrences of each called nucleotide.  
* Discovered that pileup does not output SNP calls. Modify 'heterozygousSNPfilter.pl' to extract heterozygous SNPs 'cns.final.snp', put call into hash keyed by position. Then iterate through lines of snp.pileup and count occurrences of each called nucleotide.
===Bayesian SNP Calling===
* We want to determine if a SNP called as heterozygous is actually homozygous (both alleles same as reference)
<math>H_{0}: heterozygous</math>


==Calculating False Positives==
==Calculating False Positives==

Revision as of 20:15, 9 July 2009

SNP Filtering

  • Wrote 'extractSNPlines.pl' to get pileup lines by line number (instead of grep-ing for position)
./extractSNPlines.pl notwo.pileup snoOrigLocs.txt > snp.pileup
  • Discovered that pileup does not output SNP calls. Modify 'heterozygousSNPfilter.pl' to extract heterozygous SNPs 'cns.final.snp', put call into hash keyed by position. Then iterate through lines of snp.pileup and count occurrences of each called nucleotide.

Bayesian SNP Calling

  • We want to determine if a SNP called as heterozygous is actually homozygous (both alleles same as reference)

<math>H_{0}: heterozygous</math>

Calculating False Positives

  • dbSNP is a possibly biased database, while the 1000genome SNPs are unbiased (based on sequencing data, not just looking at known SNP locations).
  • 1000genome SNP list may be smaller than list of true SNPs.
  • If there are no false positives, then the proportion of our SNPs in common with dbSNPs should be independent of presence in dbSNP.

Yesterday's SNP Frequency Measurements

  • Possible over-representation of heterozygous SNPs: much more common among presumed false positives than direct hits.
' In dbSNP Not in dbSNP
Homozygous: 55.33% 20.72%
Heterozygous: 44.67% 79.28%
Transition: 76.48% 67.17%
Transversion: 23.52% 32.83%
Total: 7377 999