This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What Is The Difference Between Samtools Mpileup And Pileup

As of samtools 0.1.8 'mpileup' appears to have superseded 'pileup' for variant detection. The samtools website now has pileup as deprecated.

So i am wondering whether i should repeat my recent SNP/INDEL calling with mpileup? But does anyone know what the improvements are using mpileup?

For those who don't already know 'pileup' and 'mpileup' can be used to create a consensus sequence using NGS mapped reads relative to a reference genome, but can also highlight variants relative to the reference.

samtools pileup mpileup snp indel
• 20,190 views • link
updated by lh3 lh3 • written by Ian Ian

2 answers

For single-sample SNP calling, they differ little. Pileup uses BAQ by default, too. The major difference comes from indel calling. Mpileup implements a more advanced method. For indel calling, do not use pileup any more.

• 55 views • link
written by lh3 lh3

It seems pileup is going away in the next release. The pileup raw output was still very useful (read coverage among others). If that feature goes away, what tool can be used as a replacement? An alternative would be to write your version. I wanted to check before moving in that direction.

• 0 views • link
written by Drio Drio

Mpileup also outputs pileup.

• 0 views • link
written by lh3 lh3

Indeed, mpileup is a pileup that supports multiple samples and the old pileup is deprecated.

Mpileup by default uses a new algorithm (BAQ). The Snp calling is now done by BcfTools, which (I think) also made some improvements.

Since generating mpileup/bcf/vcf is not such a long process, I'd recommend re-calculating the calls.

• 0 views • link
written by Pablo Pablo

Log in to answer this question.

More posts like this