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"
•
written
by
Egelbets •
1 answer
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.
Log in to answer this question.
-
Variant calling using samtools
written by Chris •
-
Consensus sequence calling with normalisation of indels
written by michael.flower.14 •
-
Does bcftools call ignore duplicate markings?
written by Egelbets •
-
bcftools command not found
written by michael.flower.14 •
-
bcftools filter
written by Gregor Rot
-
How to force bcftools to call all variants
written by Hans
-
Bcftools Call method not detecting single nucleotide deletions?
written by kennyworkman •
-
bcftools mpileup no output (empty calls.bcf file)
written by RobertUt •
-
priors for calculting genotype likelihood using RNAseq data
written by prasundutta87
-
samtools mpileup VCF output
written by igor
