Question: Snp Detection In Viruses: Does A Haploid Caller Like Freebayes Need To Be Used Instead Of Something Like Samtools/Bcftools?
gravatar for JacobS
5.2 years ago by
Cleveland, Ohio
JacobS890 wrote:

I've done variant calling for many eukaryotic organisms, and always followed a simple, highly vetted pipeline of mapping my reads to a reference, generating pileup with samtools, and calling variants with bcftools. Now I'm interested in SNP calling for viral genomes, and although I know of several published papers that use this pipeline in virology, I heard that a haploid-specific algorithm, like that used in FreeBayes may be a more accurate analysis route.

Is this really the case? Does anyone know what negative effects will occur if using bcftools with viral sequence? Are there better options than FreeBayes?

Thanks for any advice.

variant-calling snp • 2.7k views
ADD COMMENTlink modified 5.2 years ago by Andreas2.4k • written 5.2 years ago by JacobS890
gravatar for Andreas
5.2 years ago by
Andreas2.4k wrote:

Well, basically samtools' (and therefore bcftools') genotyping function assumes a diploid genome. SNV are usually limited/clipped at frequencies around 20% for that same reason, I believe. For viral genomes it's better to use programs like LoFreq, Breseq and SNVer and the like. These programs are especially useful if you're interested in low frequency variants.


ADD COMMENTlink written 5.2 years ago by Andreas2.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2409 users visited in the last hour