Entering edit mode
5.7 years ago
enxxx23
▴
280
What is the fastest (non-pseudo) aligner for RNA-seq illumina seq data today in year 2019?
What is the fastest (non-pseudo) aligner for RNA-seq illumina seq data today in year 2019?
The better choices are:
Both are very fast and highly accurate. Unless you have large differences in GC content there is no better tool - then they perform very similar. If you have a high GC-content Salmon is probably the better option due to its gcBias algorithm.
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
fastest
is a relative metric unless you do an apples-to-apples comparison on hardware you will eventually use. Choice of alignment options can have a significant impact on speed. It would be difficult to get those to align well among aligners.bbmap.sh
from BBMap suite will stand with the best of them on any given day.bwa mem
may be the smallest memory footprint aligner (~6-7G for human genome).BWA MEM is not a RNA-seq aligner by design.
Sure. While the statement is true if you are looking for a splice-aware aligner it is not applicable in this case. Though you could use it if you were dealing with bacterial RNAseq data.
True, I should have mentioned the target organisms, which in this case are eukaryotes. Retroviruses do not have RNA so doing RNA-seq on retroviruses is the only choice! ;-)
what do you mean by
I don't believe that speed only is a valuable concept. I could write an EXTREMELY fast aligner (with terrible accuracy).
Challenge accepted !
ok great. Now submit to Nature Methods
They'd accept it, for sure.
At last try to :
if read in reference: print('Aligned!')
STAR is for me, the best compromise between running time and efficiency
https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5792058/
I think STAR and HiSAT are probably similar. STAR requires less tuning, but HiSAT can be tuned to give very similar accuracy performance. The benefit of HiSAT is that it uses much less memory.