SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > SOLiD



Similar Threads
Thread Thread Starter Forum Replies Last Post
maq map problem prm36 Bioinformatics 0 03-12-2010 09:12 AM
MAQ - colorspace alignment troubles Jonathan Bioinformatics 1 01-20-2010 10:16 PM
maq map command... hcpark79 Bioinformatics 2 12-05-2009 05:29 AM
maq map -H file baohua100 Bioinformatics 1 08-07-2009 08:06 AM
maq map/match settings in colorspace der_eiskern Bioinformatics 0 07-28-2009 01:10 PM

Reply
 
Thread Tools
Old 11-29-2009, 02:31 AM   #1
Mr Mutundes
Member
 
Location: Sydney, Australia

Join Date: Jan 2009
Posts: 17
Default maq map -n in colorspace

Hi,

I have been investigating the use of MAQ 0.7.1 for aligning SOLiD reads to a reference and calling SNPs. I am wondering about the behaviour of the -n option which sets maximum mismatches allowed when using maq's map command. MAQ allows only values of 1,2 or 3 for this option. Since in colorspace a SNP is comprised of two adjacent color changes does this mean that reads with more than one SNP cannot be aligned with maq map?

Thanks for any help.

Mr M
Mr Mutundes is offline   Reply With Quote
Old 11-29-2009, 05:36 AM   #2
Chipper
Senior Member
 
Location: Sweden

Join Date: Mar 2008
Posts: 324
Default

Not sure about MAQ, but bwa works well with CS. You can set it to allow 2 mismatches in the first (good part) of the read and several in the last part and still get good speed (and more unique alignments than mapreads). But my impression after some extensive testing is that BFAST is the best option since it is able to place more low-quality reads uniquely as well as using the first color. Have not yet tested if the low-quality alignments are correct though...
Chipper is offline   Reply With Quote
Old 12-01-2009, 02:22 PM   #3
Mr Mutundes
Member
 
Location: Sydney, Australia

Join Date: Jan 2009
Posts: 17
Default

Thanks Chipper
Mr Mutundes 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 08:13 PM.


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