question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

GFF3 with CSI index shows "Error: Not a TBI file" in linear genome view

See original GitHub issue

Describe the bug

I’ve tried multiple GFF3 files which require CSI index and all of them show the same behaviour. Specifically, the track appears to be added to the configuration correctly as shown below but when I turn it on through the Linear Genome view track selector I see the following error message:

Screen Shot 2020-07-20 at 3 37 52 PM

Here is the config file:

{
  "assemblies": [
    {
      "name": "Pisum sativum Cameor v1",
      "sequence": {
        "type": "ReferenceSequenceTrack",
        "trackId": "Pisum sativum Cameor v1-ReferenceSequenceTrack",
        "adapter": {
          "type": "BgzipFastaAdapter",
          "fastaLocation": {
            "uri": "https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.P4FG/pissa.Cameor.gnm1.P4FG.genome_main.fna.gz"
          },
          "faiLocation": {
            "uri": "https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.P4FG/pissa.Cameor.gnm1.P4FG.genome_main.fna.gz.fai"
          },
          "gziLocation": {
            "uri": "https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.P4FG/pissa.Cameor.gnm1.P4FG.genome_main.fna.gz.gzi"
          }
        }
      },
      "aliases": [
        "pissa.Cameor"
      ]
    }
  ],
  "configuration": {},
  "connections": [],
  "defaultSession": {
    "name": "New Session"
  },
  "tracks": [
    {
      "type": "BasicTrack",
      "trackId": "genes",
      "name": "Genes",
      "category": [
        "Genes"
      ],
      "assemblyNames": [
        "Pisum sativum Cameor v1"
      ],
      "adapter": {
        "type": "Gff3TabixAdapter",
        "gffGzLocation": {
          "uri": "https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.ann1.7SZR/pissa.Cameor.gnm1.ann1.7SZR.gene_models_main.gff3.gz"
        },
        "index": {
          "location": {
            "uri": "https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.ann1.7SZR/pissa.Cameor.gnm1.ann1.7SZR.gene_models_main.gff3.gz.csi"
          },
          "indexType": "CSI"
        }
      }
    }
  ]
}

To Reproduce

  1. Execute the following commands in your web root:
jbrowse create jbrowse2
cd jbrowse2
jbrowse add-assembly https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.P4FG/pissa.Cameor.gnm1.P4FG.genome_main.fna.gz --name "Pisum sativum Cameor v1" --alias pissa.Cameor --type bgzipFasta
jbrowse add-track https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.ann1.7SZR/pissa.Cameor.gnm1.ann1.7SZR.gene_models_main.gff3.gz.csi --name="Genes" --trackId=genes --category="Genes"
  1. Navigate to http://localhost/jbrowse2.
  2. Add a linear Genome View (File > Add > Linear Genome View)
  3. Click the “Select Tracks” button and enable the “Genes” track. It will add the track but the error shown in the above screenshot should be visible instead of the gene glyphs.

NOTE: I’m using a docker image for easy reproducibility. You can access it here including instructions on it’s usage. If using the docker, execute the following commands after starting your container rather then those listed in step 1 above.

docker exec useJBrowse jbrowse add-assembly https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.P4FG/pissa.Cameor.gnm1.P4FG.genome_main.fna.gz --name "Pisum sativum Cameor v1" --alias pissa.Cameor --type bgzipFasta
docker exec useJBrowse jbrowse add-track https://v1.legumefederation.org/data/public/Pisum_sativum/Cameor.gnm1.ann1.7SZR/pissa.Cameor.gnm1.ann1.7SZR.gene_models_main.gff3.gz.csi --name="Genes" --trackId=genes --category="Genes"

I’ve also tested this on my debian web server with the same results.

Expected behavior

I expect to see the gene glyphs since #797 implies CSI is supported for GFF3 files (confirmed by config being correctly added).

Version:

Browser: Firefox 77.0.1 and Chrome 83.0.4103.116

$ jbrowse --version
@gmod/jbrowse-cli/0.0.1-beta.16 linux-x64 node-v10.21.0
$ node --version
v10.21.0
$ npm --version
6.14.4
$ cat /etc/os-release
NAME="CentOS Linux"
VERSION="7 (Core)"
ID="centos"
ID_LIKE="rhel fedora"
VERSION_ID="7"
PRETTY_NAME="CentOS Linux 7 (Core)"
ANSI_COLOR="0;31"
CPE_NAME="cpe:/o:centos:centos:7"
HOME_URL="https://www.centos.org/"
BUG_REPORT_URL="https://bugs.centos.org/"

CENTOS_MANTISBT_PROJECT="CentOS-7"
CENTOS_MANTISBT_PROJECT_VERSION="7"
REDHAT_SUPPORT_PRODUCT="centos"
REDHAT_SUPPORT_PRODUCT_VERSION="7"

Additional context

The public example files used are from Pisum sativum (field pea) Cameor and available for usage here: https://v1.legumefederation.org/data/public/Pisum_sativum/. I’ve also tested this with Lentil but those files are not yet public.

Issue Analytics

  • State:closed
  • Created 3 years ago
  • Comments:6

github_iconTop GitHub Comments

1reaction
cmdcolincommented, Jul 21, 2020

Super, glad that works

I might also recommend using the SvgFeatureRenderer. I used this command to load a GFF track before

jbrowse add-track https://s3.amazonaws.com/jbrowse.org/genomes/hg19/GRCh37_latest_genomic.sort.gff.gz --config '{"renderer": {"type": "SvgFeatureRenderer"}}'

With regards to that error running jbrowse upgrade with this, we will try and fix that too. I would love to use that technique in future issues if people are trying to debug their latest instance, hopefully we can get that upgrade path smoothed out. We saw that EISDIR error in an older version of our CLI but if it is still happening we will check into it

0reactions
cmdcolincommented, Jul 21, 2020

I can confirm I have locally beta 0.0.16 and you also have 0.0.16 of the cli tools so it would appear it still happens under some circumstances

Read more comments on GitHub >

github_iconTop Results From Across the Web

JBrowse 2 guide
First, open a Linear Genome View using the navigation bar (File > Add > Linear Genome View), and click on the “Select Tracks”...
Read more >
Annotating Genomes with GFF3 or GTF files - NCBI - NIH
This page describes how to create an annoated genome submission from GFF3 or GTF files, using the beta version of our process.
Read more >
samtools-announce Mailing List for SAM tools - SourceForge
Former error message would claim that there was a "truncated file or corrupt BAM index file" with no real justification. Also reset errno...
Read more >
Tabix With Gff3 - Biostar
So, I tried to sort and index it with tabix, but its not working! The sort with the line bellow seems to go...
Read more >
DRAGEN Bio-IT Platform 3.7 User Guide (document ...
The DRAGEN DNA Pipeline accelerates the secondary analysis of NGS data. For example, the time taken to process an entire human genome at...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found