View Single Post
Old 07-23-2014, 08:19 AM   #1
Location: Florida

Join Date: Feb 2013
Posts: 33
Default Samtools and intrasample diversity

Hi, Seqanswers community,

I need some help with my variant calling.

I am dealing with a monocellular organisms, and I am expecting intrasample diversity due to multiclonality. This could mean that the real allele frequency is not necessarily 50% each (in case of biallelic snp)
After calling variants with samtools/bcftools, I would like to determine intrasample allele frequencies for each individual sample. Is there a way I can quantify the support for each allele with samtools (probably using the sam or bam files) or any similar program?

Hope what I've written makes sense, please let me know if you need any more information.

mstagliamonte is offline   Reply With Quote