This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Analize RNA-seq data using Hisat2

Hi all, I'm trying to analyze RNA-seq data. I performed paired-end sequencing and to prepare the libraries I have used Illumina Stranded mRNA Prep, so my reads are stranded. After quality check I used Hisat2 to perform the alignment and I selected the option unstranded. I also tried to perform the aligment using the strand information (forward/reverse) but nothing change. So I was wondering am I missing something? Did I perform the aligment correctly?

Thank you for everyone that will help me

hisat2 rna-seq alignment

I think providing the actual command you used could help to assess if this was done correctly :)

I performed it on the Galaxy platform. Selected parameters are those in the attached picture Outline of the commands on Galaxy

Thank you for your help

It would significantly affect your output if your reads were "reverse-stranded" and you chose "forward-stranded" for alignment instead. In this scenario, it would be much more difference in the output because your reads were reverse and you chose forward stranded as alignment so, it couldn't align the right reads.

But now I think it is expected that it was not a lot of change, because it aligned the reads regardless of the orientation. but of course, for downstream analysis, it is better to choose the correct parameters as you know.

ps: if your library is reverse stranded choose RF; and if it is forward-stranded choose FR in Galaxy.

0 answers

No answers yet.

Log in to answer this question.