I'm calling the pileup like this:
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
Code:
samtools pileup -f ref.fa aln_reads.bam -c > out.pileup
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
Comment