This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How are the variants called by samtools?

Sorry if this question sounds ridiculous but I am very much puzzled as how variants are called by samtools? I mean what is the correct path?

samtools mpileup created the mpileup file or vcf file which have the variants? What does BAQ does; I suppose it is bundled with samtools, is that correct?

Also, what does bcf tools call function does. I am totally confused at this point; please help

samtools baq
• 2,550 views • link
updated by Ram Ram • written by lakhujanivijay lakhujanivijay

Did you mean this?

samtools mpileup -d $coverage_max -L $coverage_max -E -uf reference.fasta input.bam -P ILLUMINA -Q 30 | bcftools view -bvgc - > calls.bcf
bcftools view calls.bcf | vcfutils varFilter -d 30 -a 10 -1 0.001 -2 0.001 -3 0 -4 0.0001 - > calls.vcf
• 0 views • link
updated by Ram Ram • written by Max Ivon Max Ivon

Yes, can you explain it a little bit like what parameter and what tool is doing what?

• 0 views • link
written by lakhujanivijay lakhujanivijay

0 answers

No answers yet.

Log in to answer this question.

More posts like this