This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Low Reads Aligned to hg38 snp tran index of hisat2

I have a single-ended RNA-seq dataset. I have aligned the dataset to hisat2 snp and tran index for gh38 using the following code:

hisat2 --novel-splicesite-outfile SRR11594433_splicesite.txt --novel-splicesite-infile SRR11594433_splicesite.txt --time --un SRR11594433_unpair.fastq --summary-file SRR11594433_summary.txt --met-file SRR11594433_met.txt --no-unal -x <path to hisat2 transcriptomic index> -U SRR11594433.fastq -S SRR11594433.sam

The output shows that only 15% of the reads are aligned to the transcriptome by hisat2 about which I'm very concerned.

Please let me know what can be the reasons for such low portion of aligned reads and how I can solve this problem.

As indicated in the code, the dataset that I'm analyzing is SRR11594433.fastq. I hope that this can be helpful.

hisat2 rna-seq transcriptome alignment rsem

0 answers

No answers yet.

Log in to answer this question.