CONSTAX Options¶
To visualize CONSTAX options:
gian@gian-Z390-GY:~/tutorial$ constax --help
This is what CONSTAX will display on the terminal
# constax --help
usage: constax [-h] [-c CONF] [-n NUM_THREADS] [-m MHITS]
[-e EVALUE] [-p P_IDEN] [-d DB] [-f TRAINFILE]
[-i INPUT] [-o OUTPUT] [-x TAX] [-t] [-b]
[--select_by_keyword SELECT_BY_KEYWORD] [--msu_hpcc]
[-s] [--consistent] [--make_plot] [--check]
[--mem MEM] [--sintax_path SINTAX_PATH]
[--utax_path UTAX_PATH] [--rdp_path RDP_PATH]
[--constax_path CONSTAX_PATH] [--pathfile PATHFILE]
[--isolates ISOLATES]
[--isolates_query_coverage ISOLATES_QUERY_COVERAGE]
[--isolates_percent_identity ISOLATES_PERCENT_IDENTITY]
[--high_level_db HIGH_LEVEL_DB]
[--high_level_query_coverage HIGH_LEVEL_QUERY_COVERAGE]
[--high_level_percent_identity HIGH_LEVEL_PERCENT_IDENTITY]
[--combine_only] [-v]
optional arguments:
-h, --help show this help message and exit
-c CONF, --conf CONF Classification confidence threshold (default: 0.8)
-n NUM_THREADS, --num_threads NUM_THREADS
directory to for output files (default: 1)
-m MHITS, --mhits MHITS
Maximum number of BLAST hits to use, for use with -b
option (default: 10)
-e EVALUE, --evalue EVALUE
Maximum expect value of BLAST hits to use, for use
with -b option (default: 1.0)
-p P_IDEN, --p_iden P_IDEN
Minimum proportion identity of BLAST hits to use, for
use with -b option (default: 0.0)
-d DB, --db DB Database to train classifiers, in FASTA format
(default: )
-f TRAINFILE, --trainfile TRAINFILE
Path to which training files will be written (default:
./training_files)
-i INPUT, --input INPUT
Input file in FASTA format containing sequence records
to classify (default: otus.fasta)
-o OUTPUT, --output OUTPUT
Output directory for classifications (default:
./outputs)
-x TAX, --tax TAX Directory for taxonomy assignments (default:
./taxonomy_assignments)
-t, --train Complete training if specified (default: False)
-b, --blast Use BLAST instead of UTAX if specified (default:
False)
--select_by_keyword SELECT_BY_KEYWORD
Takes a keyword argument and --input FASTA file to
produce a filtered database with headers containing
the keyword with name --output (default: False)
--msu_hpcc If specified, use executable paths on Michigan State
University HPCC. Overrides other path arguments
(default: False)
-s, --conservative If specified, use conservative consensus rule (2 False
= False winner) (default: False)
--consistent If specified, show if the consensus taxonomy is
consistent with the real hierarchical taxonomy
(default: False)
--make_plot If specified, run R script to make plot of classified
taxa (default: False)
--check If specified, runs checks but stops before training or
classifying (default: False)
--mem MEM Memory available to use for RDP, in MB. 32000MB
recommended for UNITE, 128000MB for SILVA (default:
32000)
--sintax_path SINTAX_PATH
Path to USEARCH/VSEARCH executable for SINTAX
classification (default: False)
--utax_path UTAX_PATH
Path to USEARCH executable for UTAX classification
(default: False)
--rdp_path RDP_PATH Path to RDP classifier.jar file (default: False)
--constax_path CONSTAX_PATH
Path to CONSTAX scripts (default: False)
--pathfile PATHFILE File with paths to SINTAX, UTAX, RDP, and CONSTAX
executables (default: pathfile.txt)
--isolates ISOLATES FASTA formatted file of isolates to use BLAST against
(default: False)
--isolates_query_coverage ISOLATES_QUERY_COVERAGE
Threshold of sequence query coverage to report isolate
matches (default: 75)
--isolates_percent_identity ISOLATES_PERCENT_IDENTITY
Threshold of aligned sequence percent identity to
report isolate matches (default: 1)
--high_level_db HIGH_LEVEL_DB
FASTA database file of representative sequences for
assignment of high level taxonomy (default: False)
--high_level_query_coverage HIGH_LEVEL_QUERY_COVERAGE
Threshold of sequence query coverage to report high-
level taxonomy matches (default: 75)
--high_level_percent_identity HIGH_LEVEL_PERCENT_IDENTITY
Threshold of aligned sequence percent identity to
report high-level taxonomy matches (default: 1)
--combine_only Only combine taxonomy without rerunning classifiers
(default: False)
-v, --version Display version and exit (default: False)
Options details¶
-c, --conf=0.8
Classification confidence threshold, used by each classifier (0,1]. Increase for improved specificity, reduced sensitivity.
-n, --num_threads=1
Number of threads to use for parallelization. Maximum classification speed at about 32 threads. Training only uses 1 thread.
-m, --max_hits=10
Maximum number of BLAST hits to use, for use with -b option. When classifying with BLAST, this many hits are kept. Confidence for a given taxa is based on the proportion of these hits agree with that taxa. 5 works well for UNITE, 20 with SILVA (standard, not NR).
-e, --evalue=1
Maximum expect value of BLAST hits to use, for use with -b option. When classifying with BLAST, only hits under this expect value threshold are used. Decreasing will increase specificity, but decrease sensitivity at high taxonomic ranks.
-p, --p_iden=0.8
Minimum proportion identity of BLAST hits to use, for use with -b option. Minimum proportion of conserve bases to keep hit.
-d, --db
Database to train classifiers. UNITE and SILVA formats are supported. See Datasets.
-f, --trainfile=./training_files
Path to which training files will be written.
-i, --input=otus.fasta
Input file in FASTA format containing sequence records to classify.
-o, --output=./outputs
Output directory for classifications.
-x, --tax=./taxonomy_assignments
Directory for taxonomy assignments.
-t, --train
Complete training if specified. Cannot run classification without training files present, so this option is necessary at least at the first time you run CONSTAX or you changed the taxonomic referenced sequence database.
-b, --blast
Use BLAST instead of UTAX if specified. If installed with conda, this in the option that will work by default. UTAX is available from USEARCH. BLAST classification generally performs better with faster training, similar classification speed, and greater accuracy.
--msu_hpcc
If specified, use executable paths on Michigan State University HPCC. Overrides other path arguments.
--s, conservative
If specified, use conservative consensus rule (2 null = null winner. For example, if BLAST is the only algorithm that classifies OTU_135 to Family Strophariaceae while SINTAX and RDP give no classification, then no classification is reported at the rank of Family for OTU_135 in the CONSTAX taxonomy). According to our tests, works better for SILVA database to use this option.
--consistent
If specified, show if the consensus taxonomy is consistent with the real hierarchical taxonomy. In this case, a 1 indicates that all subtaxa are contained within each parent taxa. For example, the genus assigned is within the family assigned.
--make_plot
If specified, run R script to make plot of classified taxa. The plot compares how many OTUs were classifies at each rank for RDP, SINTAX, BLAST, and CONSTAX.
--check
If specified, runs checks but stops before training or classifying.
--mem
Memory available to use for RDP, in MB. 32000MB recommended for UNITE, 128000MB for SILVA. This is necessary for training the referenced databases.
--sintax_path
Path to USEARCH/VSEARCH executable for SINTAX classification. Can also be vsearch
if already on path.
--utax_path
Path to USEARCH executable for UTAX classification.
--rdp_path
Path to RDP classifier.jar
file, or classifier
if on path from RDPTools conda install.
--constax_path
Path to CONSTAX scripts.
--pathfile
File with paths to SINTAX, UTAX, RDP, and CONSTAX executables. This useful in your local CONSTAX installation, please the tutorial for how to set a pathifile up in your system.
--isolates
FASTA formatted file of isolates to use BLAST against.
--isolates_query_coverage
Threshold of sequence query coverage to report isolate matches, in percent.
--isolates_percent_identity
Threshold of aligned sequence percent identity to report isolate matches.
--high_level_db
FASTA database file of representative sequences for assignment of high level taxonomy. For this option you can use the SILVA NR99 database for SSU/16S/18S sequences or the the UNITE database for Eukaryotic ITS/28S sequences. This option is useful to match your OTUs representative sequences to a reference using a lower cutoff so you can identify for example, which sequences are Fungi and which ones are not.
--high_level_query_coverage
Threshold of sequence query coverage to report high-level taxonomy matches, in percent.
--high_level_percent_identity
Threshold of aligned sequence percent identity to report high-level taxonomy matches.
--combine_only
If specified, only reruns combine taxonomy without rerunning classifiers. Allows for changing parameters
including: -c, --conf
, -e, --evalue
, -p, --p_iden
, -s, --conservative
, --isolates_query_coverage
, --isolates_percent_identity
,
--high_level_query_coverage
, and high_level_percent_identity
.