SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
PCR Jumping (Template Switching) in Amplicon sequencing eladSeq 454 Pyrosequencing 2 04-02-2012 07:10 AM
problem about stampy komais Bioinformatics 2 03-27-2012 06:30 AM
Stampy on 454 colindaven Bioinformatics 3 01-10-2012 05:22 AM
Stampy alignment bair Bioinformatics 1 12-21-2010 04:57 AM
Issues in switching from GAII PE to Mate Pair? plattsa General 2 07-14-2010 07:10 AM

Reply
 
Thread Tools
Old 04-10-2012, 11:24 AM   #1
groundcover
Junior Member
 
Location: East coast Ground

Join Date: Apr 2012
Posts: 3
Default Slow Pileup After Switching to Stampy

I'm calling the pileup like this:

Code:
samtools pileup -f ref.fa aln_reads.bam -c > out.pileup
I used to use BWA to generate the bam file and the pileup step above was quite fast (under 30 minutes).

Then I switched to using stampy to generate the bam file and this pileup step takes about 24 hours.

Do you guys have any insight into what causes pileup to run slower or faster? Is there anything I can do to my bam file to make it faster? Can I index anything? Would it help to upgrade samtools, I'm on 0.1.7 (r510).

Thanks,

Greg
groundcover is offline   Reply With Quote
Old 04-11-2012, 04:07 AM   #2
colindaven
Senior Member
 
Location: Germany

Join Date: Oct 2008
Posts: 415
Default

Definitely upgrade samtools, the latest is 0.1.18 I believe.

This is a bit strange, as I've used both aligners with Samtools extensively. Is one Stampy aligning many more reads (I doubt it ) ?
colindaven is offline   Reply With Quote
Old 04-11-2012, 04:25 AM   #3
groundcover
Junior Member
 
Location: East coast Ground

Join Date: Apr 2012
Posts: 3
Default

I think Stampy is mapping more reads. So that could be part of the issue. But it's certainly not mapping 100x more reads. Perhaps pileup isn't O(n) though?

I have to be careful upgrading samtools as they're dropped support for pileup.

But perhaps I can try upgrading to the latest version that still has pileup.
groundcover is offline   Reply With Quote
Old 04-11-2012, 07:11 AM   #4
Rocketknight
Member
 
Location: Ireland

Join Date: Sep 2011
Posts: 86
Default

The alignments Stampy gets that BWA misses often involve several insertions or deletions, sometimes consecutively. Perhaps alignments that aren't quite as 'clean' are slowing samtools down somehow?

Also, are you aware that pileup has been replaced by mpileup in newer versions? That functionality isn't gone entirely, it's just been renamed. Or do you need some function that's only available with the original pileup?
Rocketknight is offline   Reply With Quote
Old 04-11-2012, 07:28 AM   #5
groundcover
Junior Member
 
Location: East coast Ground

Join Date: Apr 2012
Posts: 3
Default

I need the 10 column -c (consensus) output of pileup. Sadly the new mpileup command doesn't support that.

I'm now wondering if there's a way to parelize the pileup run? Maybe split it up by chromosomes and then recombine the resulting pileups.
groundcover is offline   Reply With Quote
Reply

Tags
bam, mpileup, pileup, samtools, stampy

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 06:02 AM.


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