Question: Interpreting Hisat2 alignment output
1
gravatar for caranlove
14 months ago by
caranlove10
caranlove10 wrote:

I am am aligning my RNAseq data to a reference genome using Hisat2 for the first time and I have what I am sure is a basic question. However I am still confused after reading a number of online resources.

Broadly, my pipeline goes from FastQC to rCorrector to Trimomatic to Hisat2, and I am not certain exactly how to interpret my results.

From aligning my trimmed reads I get an output something like:

23113803 reads; of these:
    23113803 (100.00%) were paired; of these:
    21488690 (92.97%) aligned concordantly 0 times
    753270 (3.26%) aligned concordantly exactly 1 time
    871843 (3.77%) aligned concordantly >1 times
    ----
    21488690 pairs aligned concordantly 0 times; of these:
      5618651 (26.15%) aligned discordantly 1 time
    ----
    15870039 pairs aligned 0 times concordantly or discordantly; of these:
      31740078 mates make up the pairs; of these:
        2583394 (8.14%) aligned 0 times
        14947960 (47.09%) aligned exactly 1 time
        14208724 (44.77%) aligned >1 times
94.41% overall alignment rate
  

I am a bit confused as to how to interpret these outputs and wonder if there is an 'ideal' percentage of reads that have been aligned 0, exactly 1, and > 1 time?

As well as how to interpret high overall alignment rates with high percentages of paired reads that aligned concordantly 0 times. Thank you in advance for any help!

rna-seq alignment hisat • 1.1k views
ADD COMMENTlink modified 14 months ago by lakhujanivijay5.3k • written 14 months ago by caranlove10

Model organism? What is the read length, is the reference genome of high quality and is this any kind of lowest input RNA-seq? Why did you use this corrector and trimmomatic? Standard RNA-seq typically does not require any pre-processing prior to alignment.

ADD REPLYlink modified 14 months ago • written 14 months ago by ATpoint40k

Thank you for the response! This is not a model organism but a closely related wild canid referenced to the CanFam3 genome. I used rcorrector and trimomatic to perform kmer filtering and remove a small amount of adapter contamination. The Illumina RNA-seq library prep was performed at a genomics core and the QC before and after sequencing did not suggest any sample or protocol issues. FastQC results suggest that read quality is good, though there is some evidence of repetitive reads. Libraries were run on a NextSeq 2x 75bp Mid Output Flow Cell.

I have also run Hisat alignment on the raw data before trimming and am pasting the alignment summary below:

26411123 reads; of these: 26411123 (100.00%) were paired; of these: 5170797 (19.58%) aligned concordantly 0 times 14261798 (54.00%) aligned concordantly exactly 1 time 6978528 (26.42%) aligned concordantly >1 times ---- 5170797 pairs aligned concordantly 0 times; of these: 51392 (0.99%) aligned discordantly 1 time ---- 5119405 pairs aligned 0 times concordantly or discordantly; of these: 10238810 mates make up the pairs; of these: 7175529 (70.08%) aligned 0 times 2030231 (19.83%) aligned exactly 1 time 1033050 (10.09%) aligned >1 times 86.42% overall alignment rate

ADD REPLYlink written 14 months ago by caranlove10
3
gravatar for lakhujanivijay
14 months ago by
lakhujanivijay5.3k
India/Ahmedabad
lakhujanivijay5.3k wrote:

ADD COMMENTlink modified 14 months ago • written 14 months ago by lakhujanivijay5.3k

Thank you, this is a helpful breakdown!

ADD REPLYlink written 14 months ago by caranlove10
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: 1014 users visited in the last hour