Splice aware or not for RNA sequence data
1
0
Entering edit mode
6 weeks ago
Nemo • 0

Hello,

I have paired fastq rna sequences. I am using hisat2 in order to align the reads to my reference based on the following command:

hisat2 --no-spliced-alignment --no-unal -p 2 -x ./reference.fasta -1 ./Data/D10_R1.fastq.gz -2 ./Data/D10_R2.fastq.gz -S D10.sam

Though I am not sure should I keep --no-spliced-alignment tag or not. For RNA data should I consider splices?

splice RNA hisat2 sequence alignment • 229 views
ADD COMMENT
1
Entering edit mode
6 weeks ago
Rob 5.5k

What is in reference.fasta? Is it the reference genome, or the collection of extracted and spliced transcripts. Using HISAT2, it is most common to align against the genome directly. In this case, you should absolutely remove --no-spliced-alignment, as you are aligning RNA-seq reads, and you expect any read crossing a junction to align in a spliced fashion to the underlying genome.

ADD COMMENT
0
Entering edit mode

Thanks Rob for your quick reply. The reference.fasta is Wuhan-Hu-1.fasta related to wuhan genome.

ADD REPLY

Login before adding your answer.

Traffic: 1195 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