SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bowtie overhang to ref sequence johanna.s Bioinformatics 0 01-17-2012 10:50 AM
how does bowtie handle n's in reads? ref sequence? nashp Bioinformatics 0 01-10-2011 04:09 PM
1000 genome data and other human ref sequence differ johnadam33 Bioinformatics 4 01-05-2011 02:16 AM
So many software for alignment!!! found Bioinformatics 5 03-03-2009 06:35 AM
SOCS: Efficient mapping of Applied Biosystems SOLiD sequence data to a ref genome... ECO Literature Watch 0 10-20-2008 07:53 PM

Reply
 
Thread Tools
Old 01-15-2010, 10:25 AM   #1
mlee
Junior Member
 
Location: sf

Join Date: Jan 2010
Posts: 2
Question 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!
mlee is offline   Reply With Quote
Old 01-16-2010, 03:34 PM   #2
spenthil
Member
 
Location: San Francisco

Join Date: Sep 2009
Posts: 44
Default

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
__________________
--
Senthil Palanisami
spenthil is offline   Reply With Quote
Old 01-18-2010, 02:20 AM   #3
dnusol
Senior Member
 
Location: Spain

Join Date: Jul 2009
Posts: 133
Default

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

Last edited by dnusol; 01-18-2010 at 02:23 AM.
dnusol is offline   Reply With Quote
Old 01-18-2010, 04:25 AM   #4
BaCh
Member
 
Location: Germany

Join Date: May 2008
Posts: 79
Default

Quote:
Originally Posted by mlee View Post
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.
BaCh 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 07:46 AM.


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