Entering edit mode
4.4 years ago
zhufc
•
0
hello, everyone! I tried to predict horizontal genes transer in a bacterial genome using HGTector 2.0b1, but failed.
Using the default parameters, it seems to unable to dertermine the taxonomy. CMD: hgtector analyze -i test.tsv -o . -t /database/hgtdb/taxdump
Traceback (most recent call last):
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/util.py", line 407, in _get_taxon
return taxdump[tid]
KeyError: nan
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/analyze.py", line 332, in assign_taxonomy
self.data[sid], self.taxdump, self.input_cov)
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/analyze.py", line 387, in infer_genome_tax
for tid in get_lineage(bestid, taxdump):
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/util.py", line 639, in get_lineage
taxon = _get_taxon(cid, taxdump)
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/util.py", line 410, in _get_taxon
'TaxID {} is not found in taxonomy database.'.format(tid))
ValueError: TaxID nan is not found in taxonomy database.
enter code here
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/data/.conda/envs/hgtector/bin/hgtector", line 100, in <module>
main()
File "/data/.conda/envs/hgtector/bin/hgtector", line 39, in main
module(args)
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/analyze.py", line 134, in __call__
self.assign_taxonomy()
File "/data/zhufangchao/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/analyze.py", line 336, in assign_taxonomy
'Please specify manually.'.format(sid))
ValueError: Cannot auto-infer taxonomy for test. Please specify manually.
Then, the taxonomy was specified, but still returned a error.
Traceback (most recent call last):
File "/data/.conda/envs/hgtector/bin/hgtector", line 100, in <module>
main()
File "/data/.conda/envs/hgtector/bin/hgtector", line 39, in main
module(args)
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/analyze.py", line 134, in __call__
self.assign_taxonomy()
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/analyze.py", line 312, in assign_taxonomy
self.input_tax = dict_from_param(self.input_tax)
File "/data/.conda/envs/hgtector/lib/python3.7/site-packages/hgtector/util.py", line 197, in dict_from_param
return dict(x.split(':') for x in param.split(','))
ValueError: dictionary update sequence element #0 has length 1; 2 is required
Can anyone help Thank you.