SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
annotation file for maq snp calling johnsequence Bioinformatics 0 10-26-2011 11:39 AM
Any software for snp calling other than MAQ? ericfit Bioinformatics 1 01-21-2011 10:40 PM
MAQ vs BWA/SAMTools - differences in SNP calling algorithm? mard Bioinformatics 2 10-20-2010 03:17 PM
percentage of SNP calling in MAQ johnsequence Bioinformatics 1 01-08-2010 10:54 AM
MAQ and SNP calling michy Bioinformatics 3 12-11-2009 09:23 AM

Reply
 
Thread Tools
Old 09-09-2008, 10:50 AM   #1
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default MAQ SNP calling

> Hi,
>
> From my understanding MAQ is tested for 20-50x coverage regions for SNP
> calling. However, we routinely have projects where depth of coverage
> from Solexa reads is ~5000x.
>
> I see quite a few 'SNPs' that are not called by MAQ.
> These are the duplicate samples, where MAQ does NOT call a SNP
>
> Pos Depth %a %c %g %t
> 64 6632 0.11 92.01 0.69 7.19
> 64 6687 0.07 90.88 0.40 8.64
>
> While for these duplicates it makes a SNP call!
> 64 8594 0.05 91.91 0.54 7.49
> 64 7889 0.13 91.15 0.49 8.23
>
> Any suggestions on what threshold/filter to use to make consistent
> calling with high depth of coverage regions?
bioinfosm is offline   Reply With Quote
Old 09-10-2008, 11:09 AM   #2
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default SNP calling in high depth/coverage solexa data

Any other options to use apart from MAQ in calling SNPs on such data?
bioinfosm is offline   Reply With Quote
Old 09-10-2008, 10:03 PM   #3
apfejes
Senior Member
 
Location: Vancouver, Canada

Join Date: Feb 2008
Posts: 236
Default

Sort of on-topic....

I'm not sure what the answer is to your question about using MAQ, but I have had the pleasure of doing SNP calls on a MAQ .map file, for which I wrote my own SNP caller. Doing your own SNP calls is actually very easy, but leaves you with the same questions about filters (though, no consistency issues, that I'm aware of.).
__________________
The more you know, the more you know you don't know. —Aristotle
apfejes is offline   Reply With Quote
Old 09-11-2008, 08:13 AM   #4
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default

Yes, MAQ SNP calling is amazingly good! I just wish it scaled to the high-depth data I have. The other alternative which I will try today is take a random subset of my data and artificially get manageable depth
bioinfosm is offline   Reply With Quote
Old 10-31-2008, 03:15 AM   #5
Rao
Member
 
Location: India

Join Date: Oct 2008
Posts: 36
Default

How Maq works for small indels?
Rao is offline   Reply With Quote
Old 10-31-2008, 02:08 PM   #6
dvh
Member
 
Location: london, uk

Join Date: Jul 2008
Posts: 35
Default

Quote:
Originally Posted by Rao View Post
How Maq works for small indels?
from the PDF manual/published paper: maq for indels needs paired end reads. maps read 1, then does gapped alignment of read 2. maq indelpe will list the indels called from the read2 gapped alignment.

try novocraft software if you want more indels from single end or paired end reads.
dvh is offline   Reply With Quote
Old 10-31-2008, 09:49 PM   #7
Rao
Member
 
Location: India

Join Date: Oct 2008
Posts: 36
Default

Quote:
Originally Posted by dvh View Post
from the PDF manual/published paper: maq for indels needs paired end reads. maps read 1, then does gapped alignment of read 2. maq indelpe will list the indels called from the read2 gapped alignment.

try novocraft software if you want more indels from single end or paired end reads.
Thanks dvh,
How about SOAP
Rao is offline   Reply With Quote
Old 11-03-2008, 08:04 AM   #8
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default

I used SOAP to detect upto 4bp indels in the solexa data. I usually use it on the reads that MAQ is not able to align to reference.
bioinfosm is offline   Reply With Quote
Old 11-05-2008, 08:43 AM   #9
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default

Another case of MAQ SNP calling that I would like to discuss here:

cns.final.snp
PKD1 37187 T C 255 255 1.12 63 62
PKD1 37189 T C 255 255 1.56 63 62

These are so close to each other, and look very likely to be False Positives. I was of the idea that MAQ throws such away, but not so..
bioinfosm is offline   Reply With Quote
Old 05-11-2009, 11:37 AM   #10
bioinfosm
Senior Member
 
Location: USA

Join Date: Jan 2008
Posts: 482
Default Missed SNP call; obvious from pileup output

I see this obvious SNP from the pileup file with 382 reads mapping the position of which 99.74% say its a C, as opposed to the reference T.

Still, the position is not reported in cns.snp file! As I understand, the maq cns2snp command on consensus.cns does not do any filtering and reports all SNPs. This ought to have been reported..

Anything I am missing here?
bioinfosm is offline   Reply With Quote
Old 01-03-2012, 01:03 AM   #11
phatjoe
Member
 
Location: Asia

Join Date: Aug 2011
Posts: 13
Talking MAQ SNP calling doubts

Quote:
Originally Posted by bioinfosm View Post
Another case of MAQ SNP calling that I would like to discuss here:

cns.final.snp
PKD1 37187 T C 255 255 1.12 63 62
PKD1 37189 T C 255 255 1.56 63 62

These are so close to each other, and look very likely to be False Positives. I was of the idea that MAQ throws such away, but not so..
Hi Bioinfoism,

This is probably 4 years late, but have you resolve the SNP calling issue with MAQ? I am actually facing the same scenario here where neighbouring SNPs were called. Any ideas on what is causing this? Also, since neighbouring SNPs were called, can the alignment be trusted?

Rookie question. But thanks in advance
phatjoe is offline   Reply With Quote
Reply

Tags
depth, maq, snp, solexa

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:12 PM.


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