Error when running Ecotag
I am getting an error when I try to run ecotag, the error goes like:
MacBook: taniavc$ ecotag -t ~/data/TAXA/nodes -R ~/data/DABA/Teleostei_mito.fasta S1_seeds_nonsingleton.fasta > S1.ecotag.fasta
Reading taxonomy dump file...
List all taxonomy rank...
Indexing taxonomy...
Indexing parent and rank...
Adding scientific name...
Adding taxid alias...
Adding deleted taxid...
Reading reference DB ... : 304600
Traceback (most recent call last):
File "/Users/taniavc/miniconda3/bin/ecotag", line 346, in <module>
taxonlink[seqid]=int(seq['taxid'])
File "src/obitools/_obitools.pyx", line 263, in obitools._obitools.BioSequence.__getitem__
File "src/obitools/_obitools.pyx", line 217, in obitools._obitools.BioSequence.getKey
KeyError: 'taxid'Info: Inside ~/data/TAXA/nodes I have: citations.dmp, division.dmp, merged.dmp, nodes.dmp, delnodes.dmp, gencode.dmp, names.dmp. I constructed my DB from Genbank, and its a Fasta file. My S1_seeds_nonsingleton.fasta is the result of the previous Obitools steps.
Any ideas for what I am doing wrong? Thanks!