Question: RNA-Seq reads maped against DNA or cDNA??
0
gravatar for sangram_keshari
16 months ago by
sangram_keshari130 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 • 769 views
ADD COMMENTlink modified 16 months ago by swbarnes25.3k • written 16 months ago by sangram_keshari130
2
gravatar for WouterDeCoster
16 months ago by
Belgium
WouterDeCoster38k wrote:

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

ADD COMMENTlink written 16 months ago by WouterDeCoster38k
2
gravatar for swbarnes2
16 months ago by
swbarnes25.3k
United States
swbarnes25.3k 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 16 months ago by swbarnes25.3k
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: 918 users visited in the last hour