Istvan Albert I tried to use bio gff to convert a genebank file into a gff, it produced this error: TypeError: '>' not supported between instances of 'NoneType' and 'int'
The gene bank accession number is EF489041.1 Thanks!
How do I create a gff(3) file from a genebank file?
Hello
Background I'm trying to analyse a RNA-seq experiment of Bacillus subtilis PY79. As part of that, I need to create an ensemble database using ensembldb (https://bioconductor.org/packages/release/bioc/html/ensembldb.html). For this I need a gff file of the genome. I tried to download it from NCBI, however, I get an error because that gff file lacks "gene_id". Since I cannot find any other gff file of that subspecies, I am now trying to generate it from the gb file.
The Problem I have a genebank (gb) file which I have downloaded from NCBI (https://www.ncbi.nlm.nih.gov/nuccore/NC_022898.1?report=genbank then send to -> file -> GeneBank (full)). I wish to convert it to a gff3 file. I have attempted several things, but no succeeded.
What I've Tried
- I tried the julia package BioJulia/GenomicAnnotations.jl, but run into an error (https://github.com/BioJulia/GenomicAnnotations.jl/issues/7). I have searched earlier cases of the same query (https://www.biostars.org/p/134589/).
- I have managed to download bp_genbank2gff3.pl, but found no instructions that helped me to run it properly.
- I have tried to find gff tools. I go to https://bmi.inf.ethz.ch/supplements/gff-tools, but only find a link to a place where I can run it online, but I'm not actually able to add any input to the online interface.
• 4,453 views
•
link
updated
by
Istvan Albert
•
written
by
torkel.loman •
•
written
by
torkel.loman •
1 answer
This usecase is one of the many for which I wrote the bio package. Get the file:
bio fetch NC_022898
convert the Genbank file to GFF like so:
bio convert NC_022898 --gff > annotations.gff
See more here:
https://www.bioinfo.help/bio-gff.html
note how much nicer and prettier the gene models made with bio are (plus are fully compatible with featureCounts):
Disclaimer: the package is still under heavy development and has not been sufficiently tested
• 0 views
•
link
written
by
Istvan Albert
open an issue or a discussion on the source code website so that this can be fixed ther
• 0 views
•
link
written
by
Istvan Albert
Log in to answer this question.
More posts like this
-
Convert whole genome mRNA file from NCBI to gff file.
written by nish_laa •
-
Converting genbank to GTF?
written by Vanish007 •
-
How to get enome feature annotation through NCBI api ?
written by demolidd77 •
-
Howcomplete STAR annotation with NCBI Genebank
written by emilie.derisoud •
-
Converting gff hints to genbank format for augustus training
written by praasu •
-
Where to download a complete Homo sapiens reference genome in gene bank (.gbk) and its.gff annotati…
written by cambridge_hockey •
-
Reading Genbank files
written by torkel.loman •
-
Sample .sam file
written by torkel.loman •
-
Possible Gene Unflattening Error In Bp_Genbank2Gff3.Pl Bioperl Script
written by Panos
-
Is This A Valid Gff3 File For Htseq?
written by biotech
The GFF file for this strain does have
geneidentifier. You should be able to use that for your counting usingfeatureCounts. This is bacterial RNAseq so things are simpler. Align with aligner of your choice and then usefeatureCountswith-g geneoption. If you choose this file then be sure to get the corresponding genome fasta file to create your indexes. That way all identifiers will match.