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.7 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.7 years ago by Devon Ryan85k • written 2.7 years ago by Peng Huang50
1
gravatar for Devon Ryan
2.7 years ago by
Devon Ryan85k
Freiburg, Germany
Devon Ryan85k 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.7 years ago by Devon Ryan85k

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 28 days ago by RamRS18k • written 2.7 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.7 years ago by Devon Ryan85k
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: 1718 users visited in the last hour