SEQanswers (
-   Bioinformatics (
-   -   Samtools and intrasample diversity (

mstagliamonte 07-23-2014 08:19 AM

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.


swbarnes2 07-23-2014 12:16 PM

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 07-23-2014 02:34 PM

Thanks, I'll try :)


All times are GMT -8. The time now is 11:56 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.