Provided by: macsyfinder_1.0.5-3build1_all 
NAME
macsyfinder - detection of macromolecular systems in protein datasets
SYNOPSIS
macsyfinder [-h] [--sequence-db SEQUENCE_DB] [--db-type
{unordered_replicon,ordered_replicon,gembase,unordered}] [--replicon-topology
{linear,circular}] [--topology-file TOPOLOGY_FILE] [--idx] [--inter-gene-max-space
INTER_GENE_MAX_SPACE INTER_GENE_MAX_SPACE] [--min-mandatory-genes-required
MIN_MANDATORY_GENES_REQUIRED MIN_MANDATORY_GENES_REQUIRED] [--min-genes-required
MIN_GENES_REQUIRED MIN_GENES_REQUIRED] [--max-nb-genes MAX_NB_GENES MAX_NB_GENES]
[--multi-loci MULTI_LOCI] [--hmmer HMMER_EXE] [--index-db INDEX_DB_EXE] [--e-value-search
E_VALUE_RES] [--i-evalue-select I_EVALUE_SEL] [--coverage-profile COVERAGE_PROFILE] [-d
DEF_DIR] [-o OUT_DIR] [-r RES_SEARCH_DIR] [--res-search-suffix RES_SEARCH_SUFFIX]
[--res-extract-suffix RES_EXTRACT_SUFFIX] [-p PROFILE_DIR] [--profile-suffix
PROFILE_SUFFIX] [-w WORKER_NB] [-v] [--log LOG_FILE] [--config CFG_FILE] [--previous-run
PREVIOUS_RUN] systems [systems ...]
DESCRIPTION
MacSyFinder is a program to model and detect macromolecular systems, genetic pathways...
in protein datasets. In prokaryotes, these systems have often evolutionarily conserved
properties: they are made of conserved components, and are encoded in compact loci
(conserved genetic architecture). The user models these systems with MacSyFinder to
reflect these conserved features, and to allow their efficient detection
OPTIONS
positional arguments:
systems
The systems to detect. This is an obligatory option with no keyword associated to
it. To detect all the protein secretion systems and related appendages: set to
"all" (case insensitive). Otherwise, a single or multiple systems can be specified.
For example: "T2SS T4P".
optional arguments:
-h, --help
show this help message and exit
Input dataset options:
--sequence-db SEQUENCE_DB
Path to the sequence dataset in fasta format.
--db-type {unordered_replicon,ordered_replicon,gembase,unordered}
The type of dataset to deal with. "unordered_replicon" corresponds to a
non-assembled genome, "unordered" to a metagenomic dataset, "ordered_replicon" to
an assembled genome, and "gembase" to a set of replicons where sequence identifiers
follow this convention: ">RepliconName SequenceID".
--replicon-topology {linear,circular}
The topology of the replicons (this option is meaningful only if the db_type is
'ordered_replicon' or 'gembase'.
--topology-file TOPOLOGY_FILE
Topology file path. The topology file allows one to specify a topology (linear or
circular) for each replicon (this option is meaningful only if the db_type is
'ordered_replicon' or 'gembase'. A topology file is a tabular file with two
columns: the 1st is the replicon name, and the 2nd the corresponding topology:
"RepliconA linear"
--idx Forces to build the indexes for the sequence dataset even if they were presviously
computed and present at the dataset location (default = False)
Systems detection options:
--inter-gene-max-space INTER_GENE_MAX_SPACE INTER_GENE_MAX_SPACE
Co-localization criterion: maximum number of components non-matched by a profile
allowed between two matched components for them to be considered contiguous. Option
only meaningful for 'ordered' datasets. The first value must match to a system, the
second to a number of components. This option can be repeated several times:
"--inter-gene-max-space T2SS 12 --inter-gene-max-space Flagellum 20"
--min-mandatory-genes-required MIN_MANDATORY_GENES_REQUIRED MIN_MANDATORY_GENES_REQUIRED
The minimal number of mandatory genes required for system assessment. The first
value must correspond to a system name, the second value to an integer. This option
can be repeated several times: "--minmandatory-genes-required T2SS 15
--min-mandatorygenes-required Flagellum 10"
--min-genes-required MIN_GENES_REQUIRED MIN_GENES_REQUIRED
The minimal number of genes required for system assessment (includes both
'mandatory' and 'accessory' components). The first value must correspond to a
system name, the second value to an integer. This option can be repeated several
times: "--min-genesrequired T2SS 15 --min-genes-required Flagellum 10"
--max-nb-genes MAX_NB_GENES MAX_NB_GENES
The maximal number of genes required for system assessment. The first value must
correspond to a system name, the second value to an integer. This option can be
repeated several times: "--max-nb-genes T2SS 5 --max-nb-genes Flagellum 10
--multi-loci MULTI_LOCI
Allow the storage of multi-loci systems for the specified systems. The systems are
specified as a comma separated list (--multi-loci sys1,sys2) default is False
Options for Hmmer execution and hits filtering:
--hmmer HMMER_EXE
Path to the Hmmer program.
--index-db INDEX_DB_EXE
The indexer to be used for Hmmer. The value can be either 'makeblastdb' or
'formatdb' or the path to one of these binary (default = makeblastb)
--e-value-search E_VALUE_RES
Maximal e-value for hits to be reported during Hmmer search. (default = 1)
--i-evalue-select I_EVALUE_SEL
Maximal independent e-value for Hmmer hits to be selected for system detection.
(default = 0.001)
--coverage-profile COVERAGE_PROFILE
Minimal profile coverage required in the hit alignment to allow the hit selection
for system detection. (default = 0.5)
Path options:
-d DEF_DIR, --def DEF_DIR
Path to the systems definition files.
-o OUT_DIR, --out-dir OUT_DIR
Path to the directory where to store results. if outdir is specified res-search-dir
will be ignored.
-r RES_SEARCH_DIR, --res-search-dir RES_SEARCH_DIR
Path to the directory where to store MacSyFinder search results directories
(default current working directory).
--res-search-suffix RES_SEARCH_SUFFIX
The suffix to give to Hmmer raw output files.
--res-extract-suffix RES_EXTRACT_SUFFIX
The suffix to give to filtered hits output files.
-p PROFILE_DIR, --profile-dir PROFILE_DIR
Path to the profiles directory.
--profile-suffix PROFILE_SUFFIX
The suffix of profile files. For each 'Gene' element, the corresponding profile is
searched in the 'profile_dir', in a file which name is based on the Gene name + the
profile suffix. For instance, if the Gene is named 'gspG' and the suffix is
'.hmm3', then the profile should be placed at the specified location and be named
'gspG.hmm3'
General options:
-w WORKER_NB, --worker WORKER_NB
Number of workers to be used by MacSyFinder. In the case the user wants to run
MacSyFinder in a multithread mode. (0 mean all cores will be used, default 1)
-v, --verbosity
Increases the verbosity level. There are 4 levels: Error messages (default),
Warning (-v), Info (-vv) and Debug.(-vvv)
--log LOG_FILE
Path to the directory where to store the 'macsyfinder.log' log file.
--config CFG_FILE
Path to a putative MacSyFinder configuration file to be used.
--previous-run PREVIOUS_RUN
Path to a previous MacSyFinder run directory. It allows one to skip the Hmmer
search step on same dataset, as it uses previous run results and thus parameters
regarding Hmmer detection. The configuration file from this previous run will be
used. (conflict with options --config, --sequence-db, --profile-suffix,
--resextract-suffix, --e-value-res, --db-type, --hmmer)
For more details, visit the MacSyFinder website and see the MacSyFinder documentation.