alignment software and ref sequence

mlee · 01-15-2010, 11:25 AM

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

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.

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

01-15-2010, 11:25 AM	#1
mlee Junior Member Location: sf Join Date: Jan 2010 Posts: 2	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!

01-16-2010, 04:34 PM	#2
spenthil Member Location: San Francisco Join Date: Sep 2009 Posts: 44	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

01-18-2010, 03:20 AM	#3
dnusol Senior Member Location: Spain Join Date: Jul 2009 Posts: 133	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 03:23 AM.