Samtools and intrasample diversity

[mstagliamonte]

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

[swbarnes2]

I'd use samtools to make a pileup file, and then parse the line of the pileup for how many ref and alt letters you see.

[mstagliamonte]

Thanks, I'll try

Regards,
Max