GitLab

Hi,

When running obi build_ref_db --taxonomy Coll01_r140_stripped/taxonomy/my_tax Coll01_r140_stripped/Coll01_r140_refs_clean Coll01_r140_stripped/Coll01_r140_db_99 I get:

2023-08-21 10:57:00,647 [build_ref_db : INFO ]  obi build_ref_db
Aligning sequences...
Building kmer tables... : Done

Error opening the taxid column when building a reference database
Traceback (most recent call last):
  File "/kyukon/home/gent/462/vsc46214/obitools3/obi3-env/bin/obi", line 4, in <module>
    __import__('pkg_resources').run_script('OBITools3==3.0.1b24', 'obi')
  File "/kyukon/home/gent/462/vsc46214/obitools3/obi3-env/lib/python3.10/site-packages/pkg_resources/__init__.py", line 722, in run_script
    self.require(requires)[0].run_script(script_name, ns)
  File "/kyukon/home/gent/462/vsc46214/obitools3/obi3-env/lib/python3.10/site-packages/pkg_resources/__init__.py", line 1561, in run_script
    exec(code, namespace, namespace)
  File "/kyukon/home/gent/462/vsc46214/obitools3/obi3-env/lib/python3.10/site-packages/OBITools3-3.0.1b24-py3.10-linux-x86_64.egg/EGG-INFO/scripts/obi", line 62, in <module>
    config[root_config_name]['module'].run(config)
  File "python/obitools3/commands/build_ref_db.pyx", line 90, in obitools3.commands.build_ref_db.run
Exception: Error building a reference database

I used the standard ncbi taxdump, and the fasta I used as import for the reference db sequences looks like this:

>AM263556
taatttttattaattctttttaaggattatttacctaacggtagttattttatatggagt
tttagatgctccatggttgccccaacaaaatagtatctcaatctatatctacaataagat
tgagagtgctataaaac
>AM263557
taatttttattaattctttttaaggattatttacctaacggtagttattttatatggagt
tttagatgctccatggttgccccaacaaaatagtacctcaatctatatctacaataagat
tgagagtgctataaaac
>AM263558
taatttttattaattctttttaaggattatttacctaacggtagttattttatatggagt
tttagatgctccatggttgccccaacaaaatagtatctcaatctatatctacaataagat
tgagagtgctataaaac
>AM263559
taatttttattaattctttttaaggattatttacctaacggtagttattttatatggagt
tttagatgctccatggttgccccaacaaaatagtacctcaatctatatctacaataagat
tgagagtgctataaaac

Any idea what I am doing wrong?

Kind regards, Sam