Question: hisat2 mm10 low transcriptome mapping rate and low uniquely mapping rate
0
gravatar for qianxu0517
6 months ago by
qianxu05170
qianxu05170 wrote:

I have rRNA-depleted RNA seq data from mouse. I wanted to map the reads to mm10 transcriptome using hisat2. My scripts are as follows: hisat2 -x /data/Hisat2Index/In \ --no-spliced-alignment \ --maxins 600 \ --mp 1,0 \ -1 UM_1_val_1.fq -2 /UM_2_val_2.fq \ -S UM_mp_1_0_hisat2.sam \ --summary-file UM_hisat2_mp_1_0_summary.txt

My mapping output is as follows: 19754290 reads; of these: 19754290 (100.00%) were paired; of these: 13774555 (69.73%) aligned concordantly 0 times 2766396 (14.00%) aligned concordantly exactly 1 time 3213339 (16.27%) aligned concordantly >1 times ---- 13774555 pairs aligned concordantly 0 times; of these: 570850 (4.14%) aligned discordantly 1 time ---- 13203705 pairs aligned 0 times concordantly or discordantly; of these: 26407410 mates make up the pairs; of these: 22991167 (87.06%) aligned 0 times 1623547 (6.15%) aligned exactly 1 time 1792696 (6.79%) aligned >1 times 41.81% overall alignment rate

I checked that all the adapters are removed. fastqc result is good Blast results shows there is no contamination. I also tried two different reference transcriptome: mm10 from genecode and esembl, the mapping rate is similar.

Can anyone gives any suggestion? Thanks very much!

snp rna-seq alignment • 320 views
ADD COMMENTlink modified 6 months ago • written 6 months ago by qianxu05170
0
gravatar for Arindam Ghosh
6 months ago by
Arindam Ghosh340
Finland
Arindam Ghosh340 wrote:

Have you checked the --rna-strandness option?

ADD COMMENTlink modified 6 months ago • written 6 months ago by Arindam Ghosh340

I have not try that, but I think I figured out the reason. By looking at the fastqc report I found that there are a severe GC bias which indicates a lot of duplicates in my fastq files. That is probably the reason for a low mapping rate.

ADD REPLYlink written 6 months ago by qianxu05170

Can you put up the FastQC GC content plot?

ADD REPLYlink written 6 months ago by Arindam Ghosh340

https://ibb.co/D5d2VHQ

ADD REPLYlink modified 6 months ago • written 6 months ago by qianxu05170

please go to this link: https://ibb.co/D5d2VHQ

ADD REPLYlink written 6 months ago by qianxu05170
0
gravatar for qianxu0517
6 months ago by
qianxu05170
qianxu05170 wrote:

Here is the figure for GC contents

ADD COMMENTlink written 6 months ago by qianxu05170
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1753 users visited in the last hour
_