SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools mpileup not filtering reads that fail mapping/base quality filters tracecakes Bioinformatics 0 04-14-2015 08:11 PM
Impact on quality of SNP calls using samtools mpileup myi Bioinformatics 3 03-03-2014 08:18 AM
samtools mpileup and bcftools not reporting dbSNP ids dmb Bioinformatics 2 08-13-2012 11:18 PM
GATK base quality recalibration suppose to keep old and new quality scores? Heisman Bioinformatics 2 10-21-2011 08:40 AM
samtools SNP Quality Scores in Depth alejo69 Bioinformatics 1 09-08-2011 09:34 AM

Reply
 
Thread Tools
Old 06-04-2015, 01:40 PM   #1
fwuffy
Junior Member
 
Location: Michigan, USA

Join Date: Jun 2015
Posts: 3
Question Why is samtools mpileup reporting all zero quality scores?

Hi everyone.

I'm calling some variants using samtools from a BWA-aligned and sorted BAM.

For some reason, the samtools mpileup is reporting all zero quality scores, but I know the base and read quality scores in the BAM are good (viewed in IGV)



samtools mpileup -uvB -t DP -f ref.fa -r chrX:48,902,600-48,902,700 mapped_sorted.bam


gives me


CHROM POS ID REF ALT QUAL
chrX 48902600 . C <X> 0

...
one for every base in the specified range

here's a variant position:
chrX 48902688 . A C,<X> 0

Base phred quality scores as reported in IGV range from 20-40 and read mapping quality is all 60.

Anyone run into this before?

Thanks!
fwuffy is offline   Reply With Quote
Old 06-04-2015, 02:59 PM   #2
lindenb
Senior Member
 
Location: France

Join Date: Apr 2010
Posts: 141
Default

cross posted : https://www.biostars.org/p/145292/
lindenb is offline   Reply With Quote
Old 06-04-2015, 05:59 PM   #3
fwuffy
Junior Member
 
Location: Michigan, USA

Join Date: Jun 2015
Posts: 3
Default

Quote:
Originally Posted by lindenb View Post
Hope that's ok.
fwuffy is offline   Reply With Quote
Old 06-04-2015, 06:09 PM   #4
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,695
Default

Quote:
Originally Posted by fwuffy View Post
Hope that's ok.
It's just helpful to note it when you post (at each location), so that people can check to see if it's already been answered sufficiently before writing a potentially quite long answer. Furthermore, it helps others who subsequently find the thread by searching.
Brian Bushnell is offline   Reply With Quote
Old 06-05-2015, 10:46 AM   #5
fwuffy
Junior Member
 
Location: Michigan, USA

Join Date: Jun 2015
Posts: 3
Default

As an update, this is only happening on chr X. Other chromosomes have proper QUAL values.
Edit- actually that was confusion on my part. It's happening on all chromosomes. Still no clue as to why, but bcftools call finds quality scores from mpileup output (presumably from data in the INFO fields) and populates that QUAL column.

Last edited by fwuffy; 06-05-2015 at 11:48 AM.
fwuffy is offline   Reply With Quote
Reply

Tags
bwa, illumina, samtools mpileup

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 02:46 PM.


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