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?
•
written
by
sangram_keshari
2 answers
Since HISAT2 performs spliced alignment you typically use the genome DNA file for building your index.
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.
Log in to answer this question.
-
Feature count no Result
written by Muhammad •
-
Building a new index for HISAT2 for RNA-Seq data analysis
written by Archit •
-
RNA-Seq read alignment genomic and transcriptomic
written by imaparna27 •
-
Building Salmon Index
written by SeroroO •
-
Set the max coverage wanted with Hisat2 ?
written by Chvatil
-
Arabidopsis thaliana RNA-Seq analysis: Is 68% transcript annotation acceptable/expected with Ensemb…
written by arctic •
-
hisat2-build --snp build personal index but mapping result doesn't change
written by liying •
-
Low read mapping efficiency with HISAT2
written by givans •
-
Arabidopsis thaliana whole genome
written by sangram_keshari
-
Altanalyze for Arabidopsis thaliana
written by mail2steff
