Mapping reads

From QuigleyWiki
Revision as of 05:26, 4 July 2015 by IanQuigley (Talk | contribs)

Jump to: navigation, search

Once you have some sequence, you'll want to map it to some reference. 100M reads that are 50 bases long each really aren't gonna help anybody until you put them in context. There are lots of mappers out there. Some are easier to use than others, and for the most part, any one of them will get the job done.


Taejoon Kwon did a nice test of how the different mappers performed on X. laevis, which is important because the pseudogenomes look very similar and that can confuse mappers sometimes. He found that RNA-STAR is probably the best, but the commands to run it are very clunky (see below).


To keep it simple, I am going to recommend first-timers go with bowtie2. It's a fine mapper, even if it's not as good as RNA-STAR. I'll also give RNA-STAR commands for you to try if you're ambitious.