Provided by: iqtree_1.6.12+dfsg-1build1_amd64
NAME
iqtree-omp - efficient phylogenetic software by maximum likelihood (multiprocessor version)
SYNOPSIS
iqtree-omp -s <alignment> [OPTIONS]
DESCRIPTION
IQ-TREE multicore version 1.6.12 for Linux 64-bit built Mar 23 2020 Developed by Bui Quang Minh, Nguyen Lam Tung, Olga Chernomor, Heiko Schmidt, Dominik Schrempf, Michael Woodhams. GENERAL OPTIONS: -? or -h Print this help dialog -version Display version number -s <alignment> Input alignment in PHYLIP/FASTA/NEXUS/CLUSTAL/MSF format -st <data_type> BIN, DNA, AA, NT2AA, CODON, MORPH (default: auto-detect) -q <partition_file> Edge-linked partition model (file in NEXUS/RAxML format) -spp <partition_file> Like -q option but allowing partition-specific rates -sp <partition_file> Edge-unlinked partition model (like -M option of RAxML) -t <start_tree_file> or -t BIONJ or -t RANDOM Starting tree (default: 99 parsimony tree and BIONJ) -te <user_tree_file> Like -t but fixing user tree (no tree search performed) -o <outgroup_taxon> Outgroup taxon name for writing .treefile -pre <PREFIX> Prefix for all output files (default: aln/partition) -nt <num_threads> Number of cores/threads or AUTO for automatic detection -ntmax <max_threads> Max number of threads by -nt AUTO (default: #CPU cores) -seed <number> Random seed number, normally used for debugging purpose -v, -vv, -vvv Verbose mode, printing more messages to screen -quiet Quiet mode, suppress printing to screen (stdout) -keep-ident Keep identical sequences (default: remove & finally add) -safe Safe likelihood kernel to avoid numerical underflow -mem RAM Maximal RAM usage for memory saving mode --runs NUMBER Number of indepedent runs (default: 1) CHECKPOINTING TO RESUME STOPPED RUN: -redo Redo analysis even for successful runs (default: resume) -cptime <seconds> Minimum checkpoint time interval (default: 60 sec) LIKELIHOOD MAPPING ANALYSIS: -lmap <#quartets> Number of quartets for likelihood mapping analysis -lmclust <clustfile> NEXUS file containing clusters for likelihood mapping -wql Print quartet log-likelihoods to .quartetlh file NEW STOCHASTIC TREE SEARCH ALGORITHM: -ninit <number> Number of initial parsimony trees (default: 100) -ntop <number> Number of top initial trees (default: 20) -nbest <number> Number of best trees retained during search (defaut: 5) -n <#iterations> Fix number of iterations to stop (default: auto) -nstop <number> Number of unsuccessful iterations to stop (default: 100) -pers <proportion> Perturbation strength for randomized NNI (default: 0.5) -sprrad <number> Radius for parsimony SPR search (default: 6) -allnni Perform more thorough NNI search (default: off) -g <constraint_tree> (Multifurcating) topological constraint tree file -fast Fast search to resemble FastTree ULTRAFAST BOOTSTRAP: -bb <#replicates> Ultrafast bootstrap (>=1000) -bsam GENE|GENESITE Resample GENE or GENE+SITE for partition (default: SITE) -wbt Write bootstrap trees to .ufboot file (default: none) -wbtl Like -wbt but also writing branch lengths -nm <#iterations> Maximum number of iterations (default: 1000) -nstep <#iterations> #Iterations for UFBoot stopping rule (default: 100) -bcor <min_corr> Minimum correlation coefficient (default: 0.99) -beps <epsilon> RELL epsilon to break tie (default: 0.5) -bnni Optimize UFBoot trees by NNI on bootstrap alignment -j <jackknife> Proportion of sites for jackknife (default: NONE) STANDARD NON-PARAMETRIC BOOTSTRAP: -b <#replicates> Bootstrap + ML tree + consensus tree (>=100) -bc <#replicates> Bootstrap + consensus tree -bo <#replicates> Bootstrap only SINGLE BRANCH TEST: -alrt <#replicates> SH-like approximate likelihood ratio test (SH-aLRT) -alrt 0 Parametric aLRT test (Anisimova and Gascuel 2006) -abayes approximate Bayes test (Anisimova et al. 2011) -lbp <#replicates> Fast local bootstrap probabilities MODEL-FINDER: -m TESTONLY Standard model selection (like jModelTest, ProtTest) -m TEST Standard model selection followed by tree inference -m MF Extended model selection with FreeRate heterogeneity -m MFP Extended model selection followed by tree inference -m TESTMERGEONLY Find best partition scheme (like PartitionFinder) -m TESTMERGE Find best partition scheme followed by tree inference -m MF+MERGE Find best partition scheme incl. FreeRate heterogeneity -m MFP+MERGE Like -m MF+MERGE followed by tree inference -rcluster <percent> Percentage of partition pairs (relaxed clustering alg.) -rclusterf <perc.> Percentage of partition pairs (fast relaxed clustering) -rcluster-max <num> Max number of partition pairs (default: 10*#partitions) -mset program Restrict search to models supported by other programs (raxml, phyml or mrbayes) -mset <lm-subset> Restrict search to a subset of the Lie-Markov models Options for lm-subset are: liemarkov, liemarkovry, liemarkovws, liemarkovmk, strandsymmetric -mset m1,...,mk Restrict search to models in a comma-separated list (e.g. -mset WAG,LG,JTT) -msub source Restrict search to AA models for specific sources (nuclear, mitochondrial, chloroplast or viral) -mfreq f1,...,fk Restrict search to using a list of state frequencies (default AA: -mfreq FU,F; codon: -mfreq ,F1x4,F3x4,F) -mrate r1,...,rk Restrict search to a list of rate-across-sites models (e.g. -mrate E,I,G,I+G,R is used for -m MF) -cmin <kmin> Min #categories for FreeRate model [+R] (default: 2) -cmax <kmax> Max #categories for FreeRate model [+R] (default: 10) -merit AIC|AICc|BIC Optimality criterion to use (default: all) -mtree Perform full tree search for each model considered -mredo Ignore model results computed earlier (default: reuse) -madd mx1,...,mxk List of mixture models to also consider -mdef <nexus_file> A model definition NEXUS file (see Manual) SUBSTITUTION MODEL: -m <model_name> DNA: HKY (default), JC, F81, K2P, K3P, K81uf, TN/TrN, TNef, TIM, TIMef, TVM, TVMef, SYM, GTR, or 6-digit model specification (e.g., 010010 = HKY) Protein: LG (default), Poisson, cpREV, mtREV, Dayhoff, mtMAM, JTT, WAG, mtART, mtZOA, VT, rtREV, DCMut, PMB, HIVb, HIVw, JTTDCMut, FLU, Blosum62, GTR20, mtMet, mtVer, mtInv Protein mixture: C10,...,C60, EX2, EX3, EHO, UL2, UL3, EX_EHO, LG4M, LG4X Binary: JC2 (default), GTR2 Empirical codon: KOSI07, SCHN05 Mechanistic codon: GY (default), MG, MGK, GY0K, GY1KTS, GY1KTV, GY2K, MG1KTS, MG1KTV, MG2K Semi-empirical codon: XX_YY where XX is empirical and YY is mechanistic model Morphology/SNP: MK (default), ORDERED, GTR Lie Markov DNA: One of the following, optionally prefixed by RY, WS or MK: 1.1, 2.2b, 3.3a, 3.3b, 3.3c, 3.4, 4.4a, 4.4b, 4.5a, 4.5b, 5.6a, 5.6b, 5.7a, 5.7b, 5.7c, 5.11a,5.11b,5.11c,5.16, 6.6, 6.7a, 6.7b, 6.8a, 6.8b, 6.17a, 6.17b,8.8, 8.10a,8.10b, 8.16, 8.17, 8.18, 9.20a,9.20b,10.12, 10.34,12.12 Non-reversible: STRSYM (strand symmetric model, synonymous with WS6.6) Non-reversible: UNREST (most general unrestricted model, functionally equivalent to 12.12) Models can have parameters appended in brackets. e.g. '-mRY3.4{0.2,-0.3}+I' specifies parameters for RY3.4 model but leaves proportion of invariant sites unspecified. '-mRY3.4{0.2,-0.3}+I{0.5} gives both. When this is done, the given parameters will be taken as fixed (default) or as start point for optimization (if -optfromgiven option supplied) Otherwise: Name of file containing user-model parameters (rate parameters and state frequencies) STATE FREQUENCY: Append one of the following +F... to -m <model_name> +F Empirically counted frequencies from alignment +FO (letter-O) Optimized frequencies by maximum-likelihood +FQ Equal frequencies +FRY, +FWS, +FMK For DNA models only, +FRY is freq(A+G)=1/2=freq(C+T), +FWS is freq(A+T)=1/2=freq(C+G), +FMK is freq(A+C)=1/2=freq(G+T). +F#### where # are digits - for DNA models only, for basis in ACGT order, digits indicate which frequencies are constrained to be the same. E.g. +F1221 means freq(A)=freq(T), freq(C)=freq(G). +FU Amino-acid frequencies by the given protein matrix +F1x4 (codon model) Equal NT frequencies over three codon positions +F3x4 (codon model) Unequal NT frequencies over three codon positions MIXTURE MODEL: -m "MIX{model1,...,modelK}" Mixture model with K components -m "FMIX{freq1,...freqK}" Frequency mixture model with K components -mwopt Turn on optimizing mixture weights (default: none) RATE HETEROGENEITY AMONG SITES: -m modelname+I A proportion of invariable sites -m modelname+G[n] Discrete Gamma model with n categories (default n=4) -m modelname*G[n] Discrete Gamma model with unlinked model parameters -m modelname+I+G[n] Invariable sites plus Gamma model with n categories -m modelname+R[n] FreeRate model with n categories (default n=4) -m modelname*R[n] FreeRate model with unlinked model parameters -m modelname+I+R[n] Invariable sites plus FreeRate model with n categories -m modelname+Hn Heterotachy model with n classes -m modelname*Hn Heterotachy model with n classes and unlinked parameters -a <Gamma_shape> Gamma shape parameter for site rates (default: estimate) -amin <min_shape> Min Gamma shape parameter for site rates (default: 0.02) -gmedian Median approximation for +G site rates (default: mean) --opt-gamma-inv More thorough estimation for +I+G model parameters -i <p_invar> Proportion of invariable sites (default: estimate) -wsr Write site rates to .rate file -mh Computing site-specific rates to .mhrate file using Meyer & von Haeseler (2003) method POLYMORPHISM AWARE MODELS (PoMo): -s <counts_file> Input counts file (see manual) -m <MODEL>+P DNA substitution model (see above) used with PoMo +N<POPSIZE> Virtual population size (default: 9) +[WB|WH|S] Sampling method (default: +WB), WB: Weighted binomial, WH: Weighted hypergeometric S: Sampled sampling +G[n] Discrete Gamma rate model with n categories (default n=4) ASCERTAINMENT BIAS CORRECTION: -m modelname+ASC Correction for absence of invariant sites in alignment SINGLE TOPOLOGY HETEROTACHY MODEL: -m <model_name>+H[k] Heterotachy model mixed branch lengths with k classes -m "MIX{m1,...mK}+H" -nni-eval <m> Loop m times for NNI evaluation (default m=1) SITE-SPECIFIC FREQUENCY MODEL: -ft <tree_file> Input tree to infer site frequency model -fs <in_freq_file> Input site frequency model file -fmax Posterior maximum instead of mean approximation CONSENSUS RECONSTRUCTION: -t <tree_file> Set of input trees for consensus reconstruction -minsup <threshold> Min split support in range [0,1]; 0.5 for majority-rule consensus (default: 0, i.e. extended consensus) -bi <burnin> Discarding <burnin> trees at beginning of <treefile> -con Computing consensus tree to .contree file -net Computing consensus network to .nex file -sup <target_tree> Assigning support values for <target_tree> to .suptree -suptag <name> Node name (or ALL) to assign tree IDs where node occurs ROBINSON-FOULDS DISTANCE: -rf_all Computing all-to-all RF distances of trees in <treefile> -rf <treefile2> Computing all RF distances between two sets of trees stored in <treefile> and <treefile2> -rf_adj Computing RF distances of adjacent trees in <treefile> TREE TOPOLOGY TEST: -z <trees_file> Evaluating a set of user trees -zb <#replicates> Performing BP,KH,SH,ELW tests for trees passed via -z -zw Also performing weighted-KH and weighted-SH tests -au Also performing approximately unbiased (AU) test ANCESTRAL STATE RECONSTRUCTION: -asr Ancestral state reconstruction by empirical Bayes -asr-min <prob> Min probability of ancestral state (default: equil freq) GENERATING RANDOM TREES: -r <num_taxa> Create a random tree under Yule-Harding model -ru <num_taxa> Create a random tree under Uniform model -rcat <num_taxa> Create a random caterpillar tree -rbal <num_taxa> Create a random balanced tree -rcsg <num_taxa> Create a random circular split network -rlen <min_len> <mean_len> <max_len> min, mean, and max branch lengths of random trees MISCELLANEOUS: -wt Write locally optimal trees into .treels file -blfix Fix branch lengths of user tree passed via -te -blscale Scale branch lengths of user tree passed via -t -blmin Min branch length for optimization (default 0.000001) -blmax Max branch length for optimization (default 100) -wsr Write site rates and categories to .rate file -wsl Write site log-likelihoods to .sitelh file -wslr Write site log-likelihoods per rate category -wslm Write site log-likelihoods per mixture class -wslmr Write site log-likelihoods per mixture+rate class -wspr Write site probabilities per rate category -wspm Write site probabilities per mixture class -wspmr Write site probabilities per mixture+rate class -wpl Write partition log-likelihoods to .partlh file -fconst f1,...,fN Add constant patterns into alignment (N=#nstates) -me <epsilon> LogL epsilon for parameter estimation (default 0.01) --no-outfiles Suppress printing output files --eigenlib Use Eigen3 library -alninfo Print alignment sites statistics to .alninfo -czb Collapse zero branches in final tree --show-lh Compute tree likelihood without optimisation