-
Notifications
You must be signed in to change notification settings - Fork 253
Description
Expected Behavior
Expected the following two commands to result in the same database:
mmseqs filtertaxseqdb swissProtSomeClasses test1 --taxon-list '4891||1075807||147549'
and
mmseqs filtertaxseqdb swissProtSomeClasses test2 --taxon-list '4891,1075807,147549'
Also expected the following command to fail/throw warning:
mmseqs filtertaxseqdb swissProtSomeClasses test3 --taxon-list '489!@!@1075807****147549'
Current Behavior
The first two commands result in two different databases:
wc -l test1
15447 test1
wc -l test2
32 test2
The third command runs without issuing any warning (it effectively does nothing to the database).
Steps to Reproduce (for bugs)
Please make sure to execute the reproduction steps with newly recreated and empty tmp folders.
- Download a small NCBI-like taxonomy database swissProtSomeClasses from here
- Run the commands above
MMseqs Output (for bugs)
Please make sure to also post the complete output of MMseqs. You can use gist.github.com for large output.
filtertaxseqdb metaeuk-regression-master/sacc_tax/swissProtSomeClasses test1 --taxon-list 4891||1075807||147549
MMseqs Version: e2510e8f6911e4340c62989aa9ba1b9c58e18d60
Compressed 0
Selected taxa 4891||1075807||147549
Subdb mode 0
Threads 8
Verbosity 3
Loading NCBI taxonomy
Loading nodes file ... Done, got 13938 nodes
Loading merged file ... Done, added 0 merged nodes.
Loading names file ... Done
Making matrix ... Done
Init RMQ ...Done
and
filtertaxseqdb metaeuk-regression-master/sacc_tax/swissProtSomeClasses test2 --taxon-list 4891,1075807,147549
MMseqs Version: e2510e8f6911e4340c62989aa9ba1b9c58e18d60
Compressed 0
Selected taxa 4891,1075807,147549
Subdb mode 0
Threads 8
Verbosity 3
Loading NCBI taxonomy
Loading nodes file ... Done, got 13938 nodes
Loading merged file ... Done, added 0 merged nodes.
Loading names file ... Done
Making matrix ... Done
Init RMQ ...Done
Context
Providing context helps us come up with a solution and improve our documentation for the future.
The help for modules that use --taxon-list allows for comma separated values:
--taxon-list STR Taxonomy ID, possibly multiple values separated by ',' []
Your Environment
Include as many relevant details about the environment you experienced the bug in.
- Git commit used (The string after "MMseqs Version:" when you execute MMseqs without any parameters): e2510e8
- Which MMseqs version was used (Statically-compiled, self-compiled, Homebrew, etc.): self-compiled
- For self-compiled and Homebrew: Compiler and Cmake versions used and their invocation: cmake version 3.5.1, c++ (Ubuntu 5.4.0-6ubuntu1~16.04.12) 5.4.0 20160609
- Server specifications (especially CPU support for AVX2/SSE and amount of system memory):
- Operating system and version: ubuntu1~16.04.12