This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Hisat2 index creating

I tried to run 

hisat2_extract_splice_sites.py genes.gtf > genes.ss
hisat2_extract_exons.py genes.gtf > genes.exon

Unfortunately, both files are empty.

What did I miss?

alignment rna-seq
• 1,150 views • link
written by m.t.lorenc • m.t.lorenc

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • masurca invalid forward file for PE library
    written by m.t.lorenc • m.t.lorenc
    Hi, I ran into the following problem: ``` > singularity run --bind /scratch/masurca/:/data /pkg/containers/masurca_4.1.0.sif masurca /data/NRotund.cfg Error line 17 of configuration file '/data/NRotund.cfg': invalid forward …
  • jellyfish histo is empty
    written by m.t.lorenc • m.t.lorenc
    I ran Jellyfish ``` > zcat V*_R1.fastp.fq.gz | jellyfish count -t 2 -C -m 21 -s 5G --quality-start=33 -o reads.jf > cat reads.jf 000001535{"alignment":8,"canonical":true,"cmdline":["count","-t","2","-C","-m","21","-s","5G","--quality-start=33","-o","reads.jf"],"counter_len":4,"exe_path":"/mnt/hpccs01/work/x/miniconda2/envs/jellyfish/bin/jellyfish","format":"binary/sorted","hostname":"cl5n014","key_len":42,"matrix1":{"c":42,"columns":[7114389000,2192333556,7375600293,3957177567,1753639420,5791419737,7660672264,3955008356,2531272710,4865888693,6095688361,6974692332,2128614574,2258809230,885814170,1703365744,5945426963,375951490,7217126865,6145592432,7418978358,2071372677,3810558711,7150470371,5456096994,5475478,7949231332,5100019741,6318414384,7141154869,3543742516,6024410566,5605745954,575704984,8124532068,1053004640,2675922405,1674046033,5535382637,1853889752,1425936082,3907867046],"identity":false,"r":33},"max_reprobe":126,"pwd":"/mnt/hpccs01/scratch/x/y/fastpDNA","reprobes":[1,1,3,6,10,15,21,28,36,45,55,66,78,91,105,120,136,153,171,190,210,231,253,276,300,325,351,378,406,435,465,496,528,561,595,630,666,703,741,780,820,861,903,946,990,1035,1081,1128,1176,1225,1275,1326,1378,1431,1485,1540,1596,1653,1711,1770,1830,1891,1953,2016,2080,2145,2211,2278,2346,2415,2485,2556,2628,2701,2775,2850,2926,3003,3081,3160,3240,3321,3403,3486,3570,3655,3741,3828,3916,4005,4095,4186,4278,4371,4465,4560,4656,4753,4851,4950,5050,5151,5253,5356,5460,5565,5671,5778,5886,5995,6105,6216,6328,6441,6555,6670,6786,6903,7021,7140,7260,7381,7503,7626,7750,7875,8001],"size":8589934592,"time":"Wed Nov …
  • Jellyfish problem with Failed to open input file 'reads.jf'
    written by m.t.lorenc • m.t.lorenc
    Hi, I can't get jellyfish running with two different ways: Method 1 jellyfish count -t 10 -C -m 21 -s 10G --quality-start=33 -o reads.jf <(zcat …
  • gtf file without exons
    written by poecile.pal • poecile.pal
    Good morning, I have `gtf` file with genome annotation (source - AUGUSTUS). There are only genes, transcripts, stop_codons, CDSs and start_codons as features. But I …
  • intersecting two GFF3 missing data from the second file
    written by Ric Ric
    Hi, I have both StringTie (*stringtie_merged.gtf.fasta.transdecoder.genome.gff*3) and BRAKER2 (*augustus.hints_utr.gff3*) GFF3 files which I would like to intersect with bedtools. Whenever StringTie covers the genome region …
  • samtools sort: truncated file. Aborting
    written by m.t.lorenc • m.t.lorenc
    Hi, I got an error by using this the below command: samtools sort -O SAM -n -@ 2 -o OZBenth2_.fastp.fq.gz.name-sorted.sam OZBenth2_.fastp.fq.gz.sam [W::sam_read1] Parse error at …
  • hisat2-align exited with value 1
    written by m.t.lorenc • m.t.lorenc
    Hi, I have got the following error samblaster: Version 0.1.24 samblaster: Inputting from stdin samblaster: Outputting to stdout Error reading _rstarts[] array: 16256, 20484 Error: …
  • vcf to venn diagram
    written by Ric Ric
    Hi, I tried to reproduce this [article](https://justinbagley.rbind.io/2018/07/25/how-to-find-common-variants-in-multiple-vcf-files/) but this command failed: $ bcftools isec -n +2 calls-16c-gt-q20.vcf calls-SA-gt-q20.vcf > isec_file1-v-2_out.vcf [main] Unrecognized command. I am …
  • eliminate empty reads in R1 and R2 fastq file
    written by cabraham03 • cabraham03
    Hi, I have a genome in two fastq files, an I have tried to determinate the deep coverage using, bwa, samtools, and BAMstats, the problem …
  • abyss-pe 2.0.0 unitig.fa, contig.fa and scaffolds.fa files empty
    written by pepijn.kooij • pepijn.kooij
    I just deployed abyss 2.0.0 and tried to run the abyss-pe pipeline, using both my own data as well as the provided test data, but …