Hi,
I aligned five 5x mate-pair libraries with different insert sizes to a 300M genome and started Pindel without the dispersed duplicates option on all BAM files. It's running for 10 days now on 30 cores and processeced around half of the genome.
Is there a way to speed up Pindel? E.g. de-duplication? Is that the expected runtime?
I read about processing the chromosomes individually. How do you deal with interchromosomal duplication in that case?
chris
2 answers
Newest Pindel code? You could split by smaller segments (say 10MB) per job with 4 cores and this will not affect interchr prediction.
If on exome data, I sometimes saw coverage spikes (several kx) in several narrow regions, this causes Pindel to slow down and have memory issue. You might use -J to exclude those regions. I have tried to solve this but have not found a smart way to handle huge coverage variation in the data.
kai
If you are interested in indels, I suggest trying BBMap. It allows alignment across long indels so you can detect them in a simple mapping+pileup rather than using an expensive realignment. It also allows an arbitrarily large number of indels or substitutions per reads, so has very high sensitivity to multiple events.
Note - while it can detect deletions of arbitrary length (depending on the maxindel flag), it cannot detect insertions more than ~50% of read length, since it only uses information from individual reads. So it's not a complete replacement for Pindel.
Log in to answer this question.
Same here.