Question: RNA-Seq reads maped against DNA or cDNA??
0
gravatar for sangram_keshari
2.2 years ago by
sangram_keshari230 wrote:

I am trying to map some scRNA-Seq data from Arabidopsis Thaliana using HISAT2. The aim is to DEGs in further.

As I read from documentation index need to be built at the start to proceed with mapping.

So my doubt is should I build the index from the ensemble DNA file or cDNA file?

rna-seq • 1.1k views
ADD COMMENTlink modified 2.2 years ago by swbarnes27.5k • written 2.2 years ago by sangram_keshari230
2
gravatar for WouterDeCoster
2.2 years ago by
Belgium
WouterDeCoster43k wrote:

Since HISAT2 performs spliced alignment you typically use the genome DNA file for building your index.

ADD COMMENTlink written 2.2 years ago by WouterDeCoster43k
2
gravatar for swbarnes2
2.2 years ago by
swbarnes27.5k
United States
swbarnes27.5k wrote:

It will be faster to align to a list of all the putative cDNAs, but today's mappers are capable of aligning to a genome, and looking for splicing, especially if it's possible to pass a gtf file which contains coordinates of known exons. But you might need a different piece of software to take your coordinates of mapped reads and convert that to gene counts. STAR aligner will do this for you, or you can use something like htseq-count.

ADD COMMENTlink written 2.2 years ago by swbarnes27.5k
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: 1770 users visited in the last hour