SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   alignment software and ref sequence (http://seqanswers.com/forums/showthread.php?t=3657)

mlee 01-15-2010 11:25 AM

alignment software and ref sequence
 
We are currently using MAQ to align our Solexa reads to our reference sequence. Does anyone know if MAQ can also tell me which parts of my reference sequence have no reads aligned to them? If not, can anyone suggest a program that can tell me this info? Thanks!

spenthil 01-16-2010 04:34 PM

Pretty sure MAQ can't, and a quick scan of the man page didn't yield anything.

Have you tried one of the visualization tools? Only one I know off the top of my head is tablet, but unfortunately haven't had a chance to play with it myself:
http://bioinf.scri.ac.uk/tablet/index.shtml

dnusol 01-18-2010 03:20 AM

Actually, if I am not wrong, Maq`s pileup output has a column with the number of reads aligning at a given position, so you can get those with a 0 coverage out

HTH,

Dave

BaCh 01-18-2010 05:25 AM

Quote:

Originally Posted by mlee (Post 12470)
We are currently using MAQ to align our Solexa reads to our reference sequence. Does anyone know if MAQ can also tell me which parts of my reference sequence have no reads aligned to them? If not, can anyone suggest a program that can tell me this info? Thanks!

MIRA does. Can even give back a list of genes and/or intergenic regions which are hit by this (if references sequence has annotations like, e.g., a GenBank). If you're working with 20 to 30m reads, have a look at it.

Beware, being a de-novo assembler, it's a lot more memory intensive than other mapping programs when used for mapping only.

Regards,
B.

Disclaimer: I wrote MIRA. It's useful to me but it might not be for you. Feedback appreciated. ;)


All times are GMT -8. The time now is 06:49 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.