help for HGTector use
0
0
Entering edit mode
5.8 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.

software error genome • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 4124 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6