Question: How to compare the alignment to genome and transcriptome
gravatar for Xiaokang ZH
11 months ago by
Xiaokang ZH40
Xiaokang ZH40 wrote:

I'm analysing RNA-Seq data of cod (Gadus morhua) which is a non-model species, meaning that the quality of its genome and annotation is not as good as model species such as human and mouse. The downstream analysis of the RNA-Seq alignment is to do differential expression analysis and pathway enrichment analysis. I see that for model species, the mainstream is to align to genome (a little bit confused here, since an answer under this post said that the mainstream of model genome is to align to transcriptome, but I found the opposite based on some papers and published RNA-Seq workflows). But how about the non-model species?

When I aligned the RAN-Seq reads using HISAT2, I tried aligning to both genome (download DNA fasta file and annotation gtf file from ENSEMBL) and transcriptome (download cDNA fasta file from ENSEMBL). I found that the alignment rate to genome was over 90% but the one to transcriptome was only around 40%. Does that mean that it's better to use genome instead of transcriptome? I also tried using Qualimap to analyse the BAM files generated from two alignments, but I didn't see much difference.

ADD COMMENTlink modified 10 months ago by ATpoint36k • written 11 months ago by Xiaokang ZH40
Please log in to add an answer.


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