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!
-
-
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:
--
Senthil Palanisami -
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,
DaveLast edited by dnusol; 01-18-2010, 03:23 AM.Comment
-
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.Originally posted by mlee View Post
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.
Comment
-
The immune system’s power comes from its genetic diversity, allowing myriad threats to be neutralized through first recognizing foreign antigens. That diversity is also what makes the immune system so difficult to study. Recent advances in sequencing technology and computational biology, however, are giving researchers new tools to understand immune responses and immune-related diseases in greater detail.
This convergence of genetics, immunology, and computation...09-01-2026, 05:41 AM
Comment