SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
loss of library DNA in the last purification, need help Anjani Sample Prep / Library Generation 3 06-27-2014 12:52 AM
mpileup -- bases missing Robby Bioinformatics 2 03-30-2012 02:48 AM
DNA loss on Bioanalyzer moedog45 Sample Prep / Library Generation 3 02-29-2012 04:03 AM
Loss of tiles with v3 chemistry hoisinjl Illumina/Solexa 8 10-05-2011 06:26 AM
Can I always get the net gain/loss of nts from an indel? sulicon Bioinformatics 2 08-02-2011 04:14 PM

Reply
 
Thread Tools
Old 06-30-2012, 12:08 PM   #1
raytobler
Junior Member
 
Location: Vienna

Join Date: Jun 2012
Posts: 1
Default loss of bases using mpileup

I have noticed a problem with mpileup (samtools) that results in a large loss of reads for certain base positions, with some bases being entirely absent. Briefly, I created a pileup file from a merged bam file which comprised 13 separate bams (using the merge option in samtools) and then created an mpileup from the same 13 bams and compared the results. The pileup contained 133M unique bases, while the mpileup contained only 111M unique bases, a reduction of 16%. Looking in more detail at the output, I noticed that not only were many bases completely absent, several others had massively reduced coverage. Moreover, in many cases the only remaining reads were those containing the alternate (non-reference) allele, all other possibilities were absent.

I was using samtools 0.1.9 in each case. All parameters were set to default for the pileup, and only -B -Q 0 were used for the mpileup. Prior to the creation of the different pileups the bam files had previously been filtered for base and mapping qualities below 20, and all ambiguously mapped reads were removed.

Any ideas on what might be causing this discrepency?
raytobler is offline   Reply With Quote
Old 07-02-2012, 12:40 AM   #2
Jane M
Senior Member
 
Location: Paris

Join Date: Aug 2011
Posts: 239
Default

You can try to add -A option I think (http://seqanswers.com/forums/showthread.php?t=21233)

Last edited by Jane M; 07-02-2012 at 12:42 AM.
Jane M is offline   Reply With Quote
Reply

Tags
merged bam, pileup, samtools

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 05:07 AM.


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