Splice aware or not for RNA sequence data
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?
1 answer
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.
Log in to answer this question.
More posts like this
-
hisat2 error
written by aj123
-
Hisat2
written by Pegasus •
-
Problem in MarkDuplicates of picard
written by Nemo •
-
Not producing sam file in hisat2
written by Nemo •
-
Adding @RG header line to bam and sam files in hisat2
written by Nemo •
-
(ERR): hisat2-align died with signal 6 (ABRT) (core dumped)
written by zhichaoli •
-
How to align multiple fastq files with Hisat2 ?
written by sarai •
-
Adding tag XS to HISAT2 sam files perior to StringTie
written by Farbod
-
hisat2 --sra-acc with paired reads producing single read output
written by avp25 •
-
Alignment stats summary by STAR aligner, please help!
written by Jackie
