This is a test version of Biostars. For the public version, visit https://www.biostars.org.
direct RNA long-reads alignment against reference genome

I am analysing Nanopore direct RNA-seq data and I am a bit confused. I have to align the RNA reads against a reference genome with minimap2 but direct RNA fastq files contain RNA reads and not cDNA. Can I align these reads against a reference genome directly? I guess it depends on library type, if reads are forward-stranded or reverse-stranded.

nanopore sequencing
• 1,036 views • link
written by frarodmar17 • frarodmar17

1 answer

Can I align these reads against a reference genome directly? I

Yes you can with the following pre-set as noted in minimap2 quick guide.

./minimap2 -ax splice -uf -k14 ref.fa reads.fa > aln.sam  # noisy Nanopore direct RNA-seq

More details here:https://github.com/lh3/minimap2?tab=readme-ov-file#map-long-splice

• 0 views • link
written by GenoMax GenoMax

Okay, thanks! I removed reads.fa and I used the index created with minimap2 in .mmi format. Is this correct?

• 0 views • link
written by frarodmar17 • frarodmar17

Log in to answer this question.

More posts like this