# STAR-long parameters for aligning RNA ONT reads to genome

Are there any suggested parameters to align ONT reads to the reference genome using STAR-long? For now, I used the parameters suggested here, but I noticed a weird behaviour.

I have RNA reads (D. melanogaster) from R7 and R9 flowcells, separately. I only selected to analyze 2D reads in pass category.

I have, respectively, 113249 reads for R7 and 40318 reads for R9. I aligned those reads and get (only!) 150 uniquely mapped reads for R7 data and 8017 uniquely mapped reads for R9 data. I tried to run again the same command on a different server with fresh compilation, but the output file is consistent with these 150 reads.

However, if I align the same with GMAP, I get 78016 uniquely mapped reads for R7 and 33523 uniquely mapped reads for R9, so I suspect that something went wrong in the alignment run.

I am aware that the two mappers behave very differently, STAR-long being more precise and preferring to report mappings of fewer reads but at better loci, and GMAP being overall less precise, trying to map the most of the reads but at not-so-good loci.

I was wondering if some of you had experience with this and could suggest me the best parameters for RNA reads from ONT?

• I guess my reply would be more appropriate as a comment - but I can't. I'm not sure if that setting of SE makes sense, but whatever. This is not necessarily a solution for the behaviour you see, but I would suggest trying github.com/lh3/minimap2 for cDNA-seq alignment. Sep 18 '17 at 8:51
• At my hand, STAR-long doesn't work well with noisy reads. It worked for iso-seq mostly because the error rate of iso-seq reads is much lower. Jan 17 '18 at 15:48

I've had great results using minimap2, particularly when combined with a pre-treatment of Canu for error correction (using minimap2 for the read-to-read mapping):

# correct reads
~/install/canu/canu-1.6/Linux-amd64/bin/canu overlapper=minimap \
-p 4T1_BC06 -d 4T1_BC06 \
-nanopore-raw workspace/pass/barcode06/fastq_runid_*.fastq