SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem calling consensus sequence with Samtools mpileup? Lspoor Bioinformatics 3 10-24-2011 07:43 AM
Compute the consensus quality, SNP quality in SAMTools lyz1030 Bioinformatics 0 04-13-2011 05:09 PM
Max Mapping / SNP / Consensus Quality in pileup files agc Bioinformatics 1 11-22-2010 06:42 AM
What cut off does GATK snp calling use for mapping quality score? foxyg Bioinformatics 0 09-25-2010 11:19 AM
SNP calling - is there an accepted Phred quality threshold? Francesco Lescai Bioinformatics 3 04-13-2010 11:51 AM

Reply
 
Thread Tools
Old 05-31-2012, 10:36 AM   #1
lg36
Member
 
Location: london

Join Date: Mar 2012
Posts: 12
Default Please Help - SNP Calling in relation to Consensus Sequence Quality

Dear All,

I've posted on this before and attracted no interest at all. So I'd be very grateful to get your opinions. I'm relatively new to NGS sequence analysis so apologies for any naive assumptions.

As I understand (perhaps misunderstand) if a SNP is deemed to be present (after all the applied SNP filters, scores etc) then (most of the time) it will appear in the consensus sequence???

Is it therefore right to think if the consensus sequence is of a 'higher quality', the SNP's too will be of a 'higher quality'??

I'm defining quality of the consensus by the ability to reflect the presence or absence of certain oligonucleotides that we know to be there from PCR studies in the lab.

What I'm doing at the moment is trying to set the parameters correctly (paired end width, etc) in order to reach a consensus sequence that is correct according to what we know in the lab. Then I was planning on calling the SNPs based on the mapped reads we used for the best consensus sequence. Is there any validity at all in this approach, or are the two things SNP calling and consensus sequence production too 'unlinked' to make it invalid????

Many thanks in advance, apologies for any misconceptions

lg36
lg36 is offline   Reply With Quote
Old 05-31-2012, 11:14 AM   #2
pag
Member
 
Location: CA, USA

Join Date: May 2012
Posts: 72
Default

Most base callers either call "n" at SNP location or pick one or the other (combination of sharpest peak and highest peak as scaled to average intensity for that base). Very few use IUPAC ambiguity codes to convey the potential difference. Quality for those locations, assuming an unbiased amplification level should be approximately 1/2 of the surrounding bases (I believe it's linearly scaled at least). So, if scripts don't already exist, you should be able to compare quality from one base to the neighboring two bases and see if the qual is significantly lower. Then you look at the trace itself (if available) to see if this is a potential SNP.

It's more difficult if it's a "chimera" or extended polymorphism or if there's a "frameshift."

I'd be interested to know if there are any chimera/SNP deconvolver programs or scripts out there for use that I can supply the reference sequence and it can let me know the alternate read. I was staring at at two overlapping reads that were significantly out-of-phase with each other yesterday and made about 50 base calls manually over an hour or two. Hardly the model of efficiency
pag is offline   Reply With Quote
Old 05-31-2012, 11:24 AM   #3
lg36
Member
 
Location: london

Join Date: Mar 2012
Posts: 12
Default

Thanks pag for your help. So, is it true to say if I can replicate what we know to be the reality of the consensus sequence (from our lab findings), then the SNPs that are output from the same mapped reads will be more reliable?
lg36 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 12:07 AM.


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