SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
samtools depth vs mpileup output question Kennels Bioinformatics 3 02-08-2012 05:48 PM
Maximum read Depth in Samtools Anjali Bioinformatics 2 01-16-2012 03:15 AM
How to generate allele depth at a specific position meher Bioinformatics 11 12-15-2011 01:48 PM
samtools SNP Quality Scores in Depth alejo69 Bioinformatics 1 09-08-2011 09:34 AM
What is read depth -D100 means in samtools? ketan_bnf Bioinformatics 1 07-28-2011 07:54 PM

Reply
 
Thread Tools
Old 05-06-2011, 12:15 PM   #1
rdeborja
Member
 
Location: Toronto

Join Date: Aug 2008
Posts: 42
Default SAMTools showing incorrect depth for a given position

I have a bam file and output from both samtools 0.1.12a and bedtools 2.11.2. The depth provided by samtools pileup is ~4000, samtools mpileup is 592 and bedtools coveragebed 12,995. The bedtools value has been confirmed when performing a samtools view for the given chromosome and position of interest.

Has anyone else encountered a descrepancy between these values as described above. Someone has mentioned something similar on the Samtools mailing list but no responses have been provided.

Cheers,
Richard
rdeborja is offline   Reply With Quote
Old 05-07-2011, 04:25 PM   #2
rdeborja
Member
 
Location: Toronto

Join Date: Aug 2008
Posts: 42
Default

Looks like it was a simple solution. Needed to upgrade to the latest version of Samtools (we were using 0.1.12) . Seems like that version of Samtools will use only properly paired reads in the BAM file. If you you use 0.1.13 or later, Samtools mpileup has a -A flag which will use anomalous read pairs in SNP/Indel calling.
rdeborja 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 01:52 AM.


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