SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BFAST and Variant Calling nexgengirl Bioinformatics 8 01-02-2013 06:03 AM
variant calling kjaja Bioinformatics 1 11-04-2011 08:16 AM
PubMed: Targeted next-generation sequencing: microdroplet PCR approach for variant de Newsbot! Literature Watch 0 08-25-2011 11:10 AM
Samtools variant calling questions Chiel Bioinformatics 2 06-07-2011 10:10 AM
variant calling using samtools -v- bcftools ksc Bioinformatics 2 04-13-2011 07:44 AM

Reply
 
Thread Tools
Old 12-21-2011, 12:16 AM   #1
krawitz
Member
 
Location: Bonn

Join Date: Feb 2010
Posts: 30
Default A heretically simple approach to variant calling

Hi everyone,

we had a look at the distribution of heterozygous allele frequencies in NGS datasets and found that their variance is larger than expected by a bionomial distribution (http://www.ncbi.nlm.nih.gov/pubmed?term=22127862). For every variant caller this means a binomial prior distribution is not the right choice and might lead to false negative calls. We also found that a simple frequency classifier (heterozygous if covered by more the 20 reads and variant allele between 14% and 86%) is more sensitive at comparable specificity for high quality data, compared to default setting of most standard calling tools.

Is anyone aware of a fast tool, that allows to apply such a frequency filter directly on a .bam file?

cheers,

peter
krawitz is offline   Reply With Quote
Old 12-21-2011, 03:58 AM   #2
Hena
Member
 
Location: Finland

Join Date: Nov 2009
Posts: 19
Default

You can ask samtools mpileup to print out the nucleotide pileups for each position in bam file. Parsing that should be fairly simple with a script.
Hena is offline   Reply With Quote
Old 02-04-2012, 05:07 AM   #3
krawitz
Member
 
Location: Bonn

Join Date: Feb 2010
Posts: 30
Default

the samtools mpileup output can be piped into VarScan to apply a coverage and frequency filter:
samtools pileup -f reference.fasta myData.bam | java -jar VarScan.v2.2.jar pileup2snp --min-coverage 20 --min-var-freq 0.14
see: http://varscan.sourceforge.net/using-varscan.html
krawitz is offline   Reply With Quote
Old 04-24-2012, 09:10 AM   #4
NGSfan
Senior Member
 
Location: Austria

Join Date: Apr 2009
Posts: 181
Default

Very nice paper.

Sometimes simpler is better.
NGSfan is offline   Reply With Quote
Old 04-25-2012, 04:01 AM   #5
Rocketknight
Member
 
Location: Ireland

Join Date: Sep 2011
Posts: 86
Default

This is awesome. Moving away from big fancy well-established tools to something like the "14-86%" rule is scary though.
Rocketknight is offline   Reply With Quote
Reply

Tags
allele distributions, variant calling

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:38 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO