This is a test version of Biostars. For the public version, visit https://www.biostars.org.
help for HGTector use

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.

software error genome

0 answers

No answers yet.

Log in to answer this question.