bcftools mpileup before bcftools call
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”
• 32 views
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.
Traffic: 2564 users visited in the last hour
Read more here: Source link