Describe the bug
Macsyfinder crashes when writing results for a model where all mandatory genes have multi_system= "true". It does it with any hmm profile and crashes even if there is only one gene in the model. If one or more of the mandatory genes don't have the multi_system option, then Macsyfinder runs correctly.
I've attached a profile definition I tested, to adapt with any hmm profile and the log file of one of the runs.
To Reproduce
Steps to reproduce the behavior:
-
The exact command lines you use 'macsyfinder --sequence-db complete_genomes.faa -m AGM aglQRS --db-type gembase --coverage-profile 0.5 --no-cut-ga -o test -w 0'
-
See error
Expected behavior
A clear and concise description of what you expected to happen.
Screenshots
exemple of a model that crashes:
error message:
Please complete the following information):
OS:
MacSyFinder Version:
MacSyFinder 2.1.6
using:
- Python 3.14.2 | packaged by conda-forge | (main, Dec 6 2025, 11:34:40) [Clang 20.1.8 ]
- MacSyLib 1.0.4
- NetworkX 3.6.1
- Pandas 3.0.0
Additional context
files attached:
macsyfinder.log
aglQRS.xml
Describe the bug
Macsyfinder crashes when writing results for a model where all mandatory genes have multi_system= "true". It does it with any hmm profile and crashes even if there is only one gene in the model. If one or more of the mandatory genes don't have the multi_system option, then Macsyfinder runs correctly.
I've attached a profile definition I tested, to adapt with any hmm profile and the log file of one of the runs.
To Reproduce
Steps to reproduce the behavior:
The exact command lines you use 'macsyfinder --sequence-db complete_genomes.faa -m AGM aglQRS --db-type gembase --coverage-profile 0.5 --no-cut-ga -o test -w 0'
See error
Expected behavior
A clear and concise description of what you expected to happen.
Screenshots
exemple of a model that crashes:
error message:
Please complete the following information):
OS:
MacSyFinder Version:
MacSyFinder 2.1.6
using:
Additional context
files attached:
macsyfinder.log
aglQRS.xml