SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools Segmentation error while using samtools index AnushaC Bioinformatics 3 06-18-2014 08:01 AM
How to calculate pi (nucleartide diversity) within subpopulation? AllenLee_china Bioinformatics 1 10-07-2012 12:32 AM
Library diversity estimate using qPCR Zapp Sample Prep / Library Generation 0 09-12-2012 02:32 AM
Low Diversity library ( 14 Ts) on HiSeq2000 mmpillai Illumina/Solexa 7 06-13-2012 06:03 AM

Reply
 
Thread Tools
Old 07-23-2014, 08:19 AM   #1
mstagliamonte
Member
 
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.

Thanks,
Max
mstagliamonte is offline   Reply With Quote
Old 07-23-2014, 12:16 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

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.
swbarnes2 is offline   Reply With Quote
Old 07-23-2014, 02:34 PM   #3
mstagliamonte
Member
 
Location: Florida

Join Date: Feb 2013
Posts: 33
Default

Thanks, I'll try

Regards,
Max
mstagliamonte is offline   Reply With Quote
Reply

Tags
heterozygous, snp analysis

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:55 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO