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.
Thanks,
Max
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.
Thanks,
Max
Comment