SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
How does "mpileup -Q" option work? dariober Bioinformatics 1 01-31-2012 01:24 PM
bowtie -m option frymor Bioinformatics 1 11-03-2011 06:58 AM
mpileup with -r option reports different number of reads at a site depending on regio rjp Bioinformatics 3 06-09-2011 10:13 PM
No enrichment option Michael L. Altshuler Illumina/Solexa 2 01-03-2011 06:48 AM
Cufflinks with the -G option gen2prot Bioinformatics 0 12-23-2010 08:42 AM

Reply
 
Thread Tools
Old 04-08-2011, 02:57 PM   #1
giverny
Member
 
Location: Canada

Join Date: Mar 2010
Posts: 10
Default mpileup -B option

Hi there
I used mpileup with the -B option to predict variations. To evaluate our pipeline, the same analysis has been done by an other lab and results have been validated in wet lab.
we share true positive prediction... good. But also we missed A LOT of true positive. I checked without the -B option and they are found. A lot of reads were not considered by using -B ... and it was wrong to do that
I'm wondering if you know which parameters I can play with to be less stringent.
Thanks for your time
giverny is offline   Reply With Quote
Old 04-08-2011, 03:15 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 868
Default

That's my experience too. Real, sanger validated SNPs will turn up with -B is selected, but won't be called otherwise.

Here's the result running

../../../../samtools/1.14/samtools mpileup -f ../1/genome6.fa region.bam > mpileup_noB.txt
598255 C 136 t,ttttttT.tttt.ttTTT.ntnn.t.t.,t.,,.TTTTtT..Tttt,T.tTtttTT.tttttttttTT.t,ttttt,t.T...tTTTtTTtttt.tT.TtTTTT.TTT.Tt,.ttttt,TtTTTt,ttttTTT^Ft !U!!!!!!!g!!>!B!!!%!BE!BDf!h!ag!hhYh!!!!!!h]*!!!f0g!%!!!!!Z!!!!!!!!!!>g!c!$!!!h!f!hhh!!!!!,!!!!!h!!g!%!!!$h!!!h!$BS!!>!!g!!!!!!b#!!!!!!!

../../../../samtools/1.14/samtools mpileup -Bf ../1/genome6.fa region.bam > mpileup_B.txt

598255 C 136 t,ttttttT.tttt.ttTTT.ntnn.t.t.,t.,,.TTTTtT..Tttt,T.tTtttTT.tttttttttTT.t,ttttt,t.T...tTTTtTTtttt.tT.TtTTTT.TTT.Tt,.ttttt,TtTTTt,ttttTTT^Ft ggThghghfghgBfBghh^dBEhBDfdhhagghhYhfgfdhgh]WffffPgf]fhggcZdfghhffhhgBghcg_hhchgfghhhefghgThPhafhfdgh^a`hZhfghhe_BS`gBghghghhhhg]hdhgfgf

the -B has for some reason destroyed the quality scores of the letters in this area, and the vcf is empty in this region. When I make it starting with -B, then that SNP pops up.

It's mixed, but the sanger was plain as day.
swbarnes2 is offline   Reply With Quote
Reply

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 11:03 PM.


Powered by vBulletin® Version 3.8.6
Copyright ©2000 - 2014, Jelsoft Enterprises Ltd.