This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Is it correct cammand for hisat indexing

/home/yog/software/hisat2-2.0.0-beta/hisat2-build -c Rs_1.0.chromosomes.fix.fasta Radish_index -p 8 --ss radish_splice --exon radish_exon

Where fasta is genome file

Radish_index is index base name

rna-seq
• 4,262 views • link
updated by Satyajeet Khare Satyajeet Khare • written by Bioinfonext Bioinfonext

Offtopic comment: what you mean is "command" instead of "cammand". I've tried to hint that in a few previous threads, but you probably missed it.

Ontopic comment: how to create a hisat index is clearly explained in the manual.

Slightlty offtpic comment: you are using an outdated version of hisat2

• 0 views • link
written by WouterDeCoster WouterDeCoster

Thanks for right suggest. trying to switch to latest version.

• 0 views • link
written by Bioinfonext Bioinfonext

1 answer

Small correction. You can remove -c

/home/yog/software/hisat2-2.0.0-beta/hisat2-build Rs_1.0.chromosomes.fix.fasta Radish_index -p 8 --ss radish_splice --exon radish_exon

Also, you may want to try without --ss and --exon options as these options will need large RAM.

• 0 views • link
written by Satyajeet Khare Satyajeet Khare

Error while running Hisat2:

[root@psgl hisat2]# /home/yog/software/hisat2-2.0.0-beta/hisat2 -p 8 -x Radish_index -1 /data/SNU_work/Analysis/unmapped/unmapped/unmapped_R1.fq -2 /data/SNU_work/Analysis/unmapped/unmapped/unmapped_R2.fq --max-intronlen 200079

(ERR): hisat2-align died with signal 11 (SEGV) (core dumped)

• 0 views • link
written by Bioinfonext Bioinfonext

Can you try without --max-intronlen parameter?

How was your index? Did hisat2-build run without error? If the issue persists, can you build index again and try alignment?

Also, though not at the time of indexing, you can provide --ss at the time of alignment such as...

hisat2 -p 10 --ss splicesites.txt -x genome_index -1 R1.fastq -2 R2.fastq 

hisat2 -p 10 --known-splicesite-infile splicesites.txt -x genome_index -1 R1.fastq -2 R2.fastq
• 0 views • link
written by Satyajeet Khare Satyajeet Khare

Thanks ...I installed the latest version and now it running well.

• 0 views • link
written by Bioinfonext Bioinfonext

Log in to answer this question.

More posts like this
  • Python scripts to map multiple pair end files to the genome
    written by Bioinfonext Bioinfonext
    I read python basic codes and I am fimiliar with for loops, But not able to understand how to make a python programme to run …
  • Error building database in snpEff
    written by Bioinfonext Bioinfonext
    I am getting some error while building database in SnpEff: java -jar /home/yog/software/snpEff/snpEff.jar build -gtf22 -v genome 00:00:00 SnpEff version SnpEff 4.3k (build 2017-03-29 17:16), …
  • How to get GFF file in genome guided assembly using StringTie
    written by Bioinfonext Bioinfonext
    For genome guide assembly, I used the below command: /home/yog/software/stringtie-1.3.2.Linux_x86_64/stringtie -p 8 -G gene.gtf -o genome.gtf -m 500 -v mapped.sorted.bam genome gtf file generated, please …
  • How to understand Hisat2 aligment result
    written by Bioinfonext Bioinfonext
    I map pair end read using Hisat2 to genome: [root@psgl hisat2]# /home/yog/software/hisat2-2.0.4/hisat2 -p 8 -x Radish_index -1 /data/SNU_work/Analysis/unmapped/unmapped/hisat2/unmap_Read1.fastq -2 /data/SNU_work/Analysis/unmapped/unmapped/hisat2/unmap_Read2.fastq -S unmap.sam Result: 101038635 reads; …
  • Genome guided assembly using cufflinks
    written by Bioinfonext Bioinfonext
    I am using this command to do assembly of mapped read on reference genome using the hisat2: /home/yog/software/cufflinks-2.2.1.Linux_x86_64/cufflinks -g gene.gtf -p 8 -o /data/SNU_work/Analysis/unmapped/unmapped/hisat2/assembly/ mapped.bam …
  • Assembly using cufflinks
    written by Bioinfonext Bioinfonext
    I am using this command to do assembly of mapped read on refernece using the hisat2: /home/yog/software/cufflinks-2.2.1.Linux_x86_64/cufflinks out_mapped.bam -G gene.gtf -p 8 -o /data/SNU_work/Analysis/unmapped/unmapped/hisat2/assembly/ But …
  • Error while doing indexing of fasta file using SAMTOOL faidx
    written by Bioinfonext Bioinfonext
    I am trying to do indexing of a CDS fasta file, I am geeting error: [root@psgl BWA]# /home/yog/software/samtools-1.3.1/samtools faidx Radish_cds.fasta [fai_build_core] different line length in …
  • To extract unmapped pair end reads
    written by Bioinfonext Bioinfonext
    Actually I map raw reads to CDS of Draft genome where all genes CDS is not present. I used below cammand lines for pair end …
  • Error while running BWA
    written by Bioinfonext Bioinfonext
    I build Indexing using this cammand: /home/yog/software/bwa-0.7.15/bwa index Rs_1.0.cds.fasta After that I am running BWA cammand for mapping raw reads to indexed transcripts using this …
  • GMAP: a genomic mapping and alignment program for mRNA and EST sequences
    written by Bioinfonext Bioinfonext
    Deal all, I installed GMAP, and build indexing of draft genome using gmap_build by using this cammand. /home/yog/software/gmap-2016-11-07/util/gmap_build -D /home/yog/database/ -d custard_genome Custard_1.0.chromosomes.fix.fasta Now want …