RNA-Seq read alignment genomic and transcriptomic
Hi,
I am trying to align my RNA-seq datasets to the human reference genome and transcriptome. However, I am not sure if it is needed and if yes, what is the purpose of --exon parameter in HISAT2 command for index building.
Also, how can I build index for human transcriptome (considering cDNA for index build) and align my reads to it, using HISAT2?
Any help would be appreciated. Thanks.
• 1,619 views
•
link
written
by
imaparna27 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
cfRNA-seq data
written by frarodmar17 •
-
Bulk RNAseq Salmon index building which transcriptome to use
written by Orange •
-
Do indexes made with cDNA fasta file contain only mRNAs?
written by bioinfo •
-
Strandedness of RNA-seq results
written by jamesymtang •
-
Building a new index for HISAT2 for RNA-Seq data analysis
written by Archit •
-
RNA read alignment using transcriptome
written by evelyn
-
What is the differences in index building by hisat with or without exon, splicing variant files?
written by mathavanbioinfo
-
Novel lncRNA identification and quantification RNA seq?
written by mathavanbioinfo
-
HISAT2 Indexing using annotation for Rattus_norvegicus
written by neranjan •
-
hisat2-build index of Arabidopsis_thaliana, met the problem "Profiling timer expired"
written by rain.sea.he •
Hi, i have not used HISAT2, but i aligned my data using STAR. it's RNASeq aligner but it's required 16GB of RAM or above. first step is generating Genome Index "assuming it's Human data download Reference data with it's annotation and i highly recommend you to use GeneCod database 'https://www.gencodegenes.org"
*genomeDir= is where to placed indexed genome *sjdbOverhang= specified the length of the genomics sequenced around annotated junctions for construction the splicing junction database it's correspond to the read length of you RNASeq data (value of 100 will usually 'Not allows' work with most cases) * this process will take around 30 Min. "for me" and it's dependent on machine specifications
second step Mapping: *For PE:
*For advance option see STAR manual 'https://physiology.med.cornell.edu/faculty/skrabanek/lab/angsd/lecture_notes/STARmanual.pdf'
Now you are done with STAR and you can proceed with GATK for Alignment post processing and variant calling, i highly recommend you to follow GATK best practices 'https://gatk.broadinstitute.org/hc/en-us/articles/360035531192-RNAseq-short-variant-discovery-SNPs-Indels-'
i hope it's helpful!
@a.alnawfal.1992 Thanks, you suggestion has been so much helpful. Now, I am doing my hands-on on STAR.
However, I was able solve it through Hisat2 as well. I just downloaded the cDNA fasta files from Ensembl and built the index using
hisat2-build(same as genome index building). Transcriptome assembly before quantification is an important step before quantification, when differential expression is the goal.