Fastest (non-pseudo) aligner for RNA-seq illumina seq data (year 2019)
1
0
Entering edit mode
5.1 years ago
enxxx23 ▴ 280

What is the fastest (non-pseudo) aligner for RNA-seq illumina seq data today in year 2019?

RNA-Seq • 2.7k views
ADD COMMENT
1
Entering edit mode

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).

ADD REPLY
0
Entering edit mode

BWA MEM is not a RNA-seq aligner by design.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

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! ;-)

ADD REPLY
0
Entering edit mode

what do you mean by

Retroviruses do not have RNA

ADD REPLY
0
Entering edit mode

I don't believe that speed only is a valuable concept. I could write an EXTREMELY fast aligner (with terrible accuracy).

ADD REPLY
1
Entering edit mode

Challenge accepted !

ADD REPLY
3
Entering edit mode
for read in SeqIO.parse("reads.fastq.gz"):
    pass
ADD REPLY
0
Entering edit mode

ok great. Now submit to Nature Methods

ADD REPLY
0
Entering edit mode

They'd accept it, for sure.

ADD REPLY
0
Entering edit mode

At last try to :

read.seq == reference
ADD REPLY
0
Entering edit mode

if read in reference: print('Aligned!')

ADD REPLY
0
Entering edit mode

STAR is for me, the best compromise between running time and efficiency

https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5792058/

ADD REPLY
1
Entering edit mode

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.

ADD REPLY
1
Entering edit mode
5.1 years ago

The better choices are:

  1. Salmon
  2. Kallisto

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.

ADD COMMENT
0
Entering edit mode

OP explicitly mentioned being looking for aligners, not pseudo-aligners...

ADD REPLY
0
Entering edit mode

Woops - Missed the -non :D

ADD REPLY

Login before adding your answer.

Traffic: 2436 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6