Incorrect gbk files when genes identifiers are numbers

[cfrioux]

• eggnog2gbk version: 0.1.0
• Python version: 3.7.7
• Operating System: MacOS 10.15.7

Description

Running emapper2gbk in genes mode with gene identifiers consisting of numbers does not create all the GBK features (translation etc.). There is no crash, a gbk is created but it lacks some important information.

What I Did

emapper2gbk genes -fn bin.fna -fp bin.faa -o bin.gbk -n "Prevotella" -a bin.tsv

LOCUS       _10007119               3225 bp    DNA              BCT 08-MAR-2022
DEFINITION  Prevotella genome.
ACCESSION   10007119
VERSION     10007119
KEYWORDS    Prevotella.
SOURCE      .
  ORGANISM  Prevotella
            Bacteria; Bacteroidetes; Bacteroidia; Bacteroidales; Prevotellaceae.
FEATURES             Location/Qualifiers
     source          1..3225
                     /scaffold="10007119"
                     /db_xref="taxon:838"
     gene            2..3225
                     /locus_tag="gene_10007119"
     CDS             2..3225
                     /locus_tag="gene_10007119"
ORIGIN
        1 atgaaagatc aaaatattaa gaaggtgttg ctcctcggct ccggtgcgtt gaagatcggt
       61 gaggccggcg agttcgacta ttccggttca caggcactca aggcgctgcg tgaggaaggc
      121 gtctacacgg tgctcatcaa tcctaatatc gccaccgtgc agacctccga gggcgtggcc
     [...]
//

When adding a prefix to all identifiers, a correct gbk is created:

LOCUS       g10007119               3225 bp    DNA              BCT 08-MAR-2022
DEFINITION  Prevotella genome.
ACCESSION   g10007119
VERSION     g10007119
KEYWORDS    Prevotella.
SOURCE      .
  ORGANISM  Prevotella
            Bacteria; Bacteroidetes; Bacteroidia; Bacteroidales; Prevotellaceae.
FEATURES             Location/Qualifiers
     source          1..3225
                     /scaffold="g10007119"
                     /db_xref="taxon:838"
     gene            2..3225
                     /locus_tag="g10007119"
     CDS             2..3225
                     /locus_tag="g10007119"
                     /gene="carB"
                     /EC_number="6.3.5.5"
                     /dbxref="KEGG:R00256"
                     /dbxref="KEGG:R00575"
                     /dbxref="KEGG:R01395"
                     /dbxref="KEGG:R10948"
                     /dbxref="KEGG:R10949"
                     /translation="MKDQNIKKVLLLGSGALKIGEAGEFDYSGSQALKALREEGVYTVL
                     INPNIATVQTSEGVADQIYFLP[...]"
ORIGIN
        1 atgaaagatc aaaatattaa gaaggtgttg ctcctcggct ccggtgcgtt gaagatcggt
       61 gaggccggcg agttcgacta ttccggttca caggcactca aggcgctgcg tgaggaaggc
      121 gtctacacgg tgctcatcaa tcctaatatc gccaccgtgc agacctccga gggcgtggcc
      [...]
\\

[cfrioux]

This case is partly accounted for already in
https://github.com/AuReMe/emapper2gbk/blob/master/emapper2gbk/genes_to_gbk.py#L131

[ArnaudBelcour]

This should be fixed in commit 9025e08 for emapper2gbk genes.

But there is still work to fix it for emapper2gbk genomes.

[ArnaudBelcour]

A first fix for genomes and genes have been made in 0.2.0.