Question: Could HISAT2-build build Hierarchical Graph FM index (HGFM) for a reference genome plus transcripts for Rat (for example, Rn6)
0
gravatar for Peng Huang
2.5 years ago by
Peng Huang50
China/HangZhou/ZheJiang University
Peng Huang50 wrote:

Hello everyone! I am doing rat RNA-seq analysis, previously using tophat2, as HISAT2 is a successor for tophat2, I plan to switch to it. And no pre-built indexs exist for Rat, hence I am going to build manually but I am wondering is it necessary and available to build Hierarchical Graph FM index (HGFM) for a reference genome plus transcripts for rat? For example, with reference genome plus transcriptome sequence of rat as input for HISAT2-build ? 

transcriptome rna-seq hisat2 • 1.6k views
ADD COMMENTlink modified 2.5 years ago by Devon Ryan82k • written 2.5 years ago by Peng Huang50
1
gravatar for Devon Ryan
2.5 years ago by
Devon Ryan82k
Freiburg, Germany
Devon Ryan82k wrote:

You would only include the reference genome, not the transcriptome sequence. You can also provide the splice junctions and there's a python script that comes with hisat2 to generate that file for you.

ADD COMMENTlink written 2.5 years ago by Devon Ryan82k

Yes, i know what you mean, i just saw the description for the option "----known-splicesite-infile" before I wrote this post, because it said: 

"Note that it is better to use indexes built using annotated transcripts (such as genome_tran or genome_snp_tran), which works better than using this option"

so, is it still invalid to build such indexes like that?

or should i just use "-ss"、"--exon" options to provide a list of exons and splice sites exacted from GTF file to build indexes?

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by Peng Huang50

That's referring to building the index with the --ss and such options. That's better than building the index without them and then giving hisat2 the same options.

ADD REPLYlink written 2.5 years ago by Devon Ryan82k
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: 1599 users visited in the last hour