bcftools mpileup before bcftools call
1
1
Entering edit mode
3.1 years ago
Egelbets ▴ 30

I want to variant call using bcftools call. However, when researching, I see a lot of people running bcftools mpileup before doing the actual variant calling with call. For example (from here):

bcftools mpileup -f reference.fa alignments.bam | bcftools call -mv -Ob -o calls.bcf

Could someone explain to me why this is done? I also don't understand what is meant with "mpileup generates genotype likelihoods at each genomic position with coverage"

bcftools mpileup call variant-calling • 1.2k views
ADD COMMENT
2
Entering edit mode
3.1 years ago

The answer is exactly what you quoted, though the instructions could be more explicit.

A variant is simply a base with multiple, reasonably likely alternatives.

To find what varies we need to look at "every" base. Thus the 'pileup' will generate the likelihood at every base out of which the "call" decides which base passes the various thresholds and quality measures.

ADD COMMENT

Login before adding your answer.

Traffic: 1046 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6