SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
TopHat samtools snp calling JasonR Bioinformatics 5 05-13-2013 09:43 AM
Samtools snp calling -- old version seems to be better for my dataset? wwmm933 Bioinformatics 9 12-09-2011 06:21 AM
Samtools SNP calling vidhya Bioinformatics 3 04-07-2011 07:17 AM
Samtools: SNP calling from only coding regions newbietonextgen Bioinformatics 5 12-10-2010 02:42 PM
SAMTOOLS SNP calling question harrb Bioinformatics 2 12-10-2010 07:37 AM

Reply
 
Thread Tools
Old 02-04-2010, 09:39 AM   #1
Jan
Junior Member
 
Location: Phoenix, AZ

Join Date: May 2009
Posts: 6
Default SAMtools and SNP calling

I'm a biologist by training and not a bioinformatics person so I'm hoping someone can help me understand the SNP calling parameters with SAMtools. We are aligning with BFAST and using SAMtools for our SNP caller. What I'm trying to understand is what the default settings for the SNP caller are (e.g. the minor allele must present X% for the call to be counted as a heterozygote) and how I can modify how the SNPs are called (I've studied the MAQ manual, but it hasn't helped me). We ran a SNP chip for validation - the majority of the discordant calls we observed between the SNP chip and our sequencing were hets not called by the SNP caller and when I go into IGV or our in-house visualization program I see the minor alleles present in a high enough percentages that I would have called them as heterozygotes - agreeing with the SNP chip. I'm working in tumor samples which can be "contaminated" with normal cells so I expect to see less than "textbook" heterozygotes presenting themselves and want to make sure I capture them, but not overly elevate my false positive rate either.
Jan is offline   Reply With Quote
Old 02-05-2010, 07:46 PM   #2
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

Samtools/maq are not designed for very deep coverage. In that case, please take the pileup output and design your own algorithm. For the SNP calling algorithm, it is not as simple as merely setting a threshold. You should read the maq/soapsnp paper for more details.
lh3 is offline   Reply With Quote
Old 09-16-2010, 02:01 PM   #3
plnagy
Junior Member
 
Location: New York

Join Date: Sep 2010
Posts: 2
Default

Even in samples where the variations are expected to be present in every cell , we experience that the program does not make the appropriate calls. We would like to verify what its settings are for calling heterozygousity and adjust them to have greater sensitivity while maintaining specificity. Is there a way to do that?
plnagy is offline   Reply With Quote
Reply

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 02:06 AM.


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