IQTREE3-MPI(1)
π¦ iqtree3-mpi(1) β User Commands π¦
okay bestie, so this is the manual page for iqtree3-mpi, section 1 of the User Commands. it's giving efficient phylogenetic software by maximum likelihood, but it's the multiprocessor version. slay. π
πΈ NAME πΈ
iqtree3-mpi β efficient phylogenetic software by maximum likelihood (multiprocessor version). periodt. π
β¨ SYNOPSIS β¨
iqtree [-s ALIGNMENT] [-p PARTITION] [-m MODEL] [-t TREE] ...
π« DESCRIPTION π«
IQ-TREE MPI version 3.1.1 for Linux x86 64-bit built Jul 28 2026. no cap.
Developed by Bui Quang Minh, Thomas Wong, Nhan Ly-Trong, Huaiyan Ren.
Contributed by Lam-Tung Nguyen, Dominik Schrempf, Chris Bielow, Olga Chernomor, Michael Woodhams, Diep Thi Hoang, Heiko Schmidt. iconic. π
β¨ GENERAL OPTIONS β¨
-h,--helpβ print (more) help usages. π-s FILE[,...,FILE]β PHYLIP/FASTA/NEXUS/CLUSTAL/MSF alignment file(s).-s DIRβ directory of alignment files.--seqtype STRINGβBIN,DNA,AA,NT2AA,CODON,MORPH(default: auto-detect).-t FILE|PARS|RANDβ starting tree (default: 99 parsimony and BIONJ).-o TAX[,...,TAX]β outgroup taxon (list) for writing.treefile.--prefix STRINGβ prefix for all output files (default: aln/partition).--seed NUMβ random seed number, normally used for debugging purpose.--safeβ safe likelihood kernel to avoid numerical underflow.--mem NUM[G|M|%]β maximal RAM usage in GB | MB | %.--runs NUMβ number of indepedent runs (default: 1). (yes, "indepedent" is in the original, keep it)-v,--verboseβ verbose mode, printing more messages to screen.-V,--versionβ display version number.--quietβ quiet mode, suppress printing to screen (stdout).-fconst f1,...,fNβ add constant patterns into alignment (N=no. states).--epsilon NUMβ likelihood epsilon for parameter estimate (default 0.01).-T NUM|AUTOβ no. cores/threads or AUTO-detect (default: 1).--threads-max NUMβ max number of threads for-T AUTO(default: all cores).
π CHECKPOINT π
--redoβ redo both ModelFinder and tree search.--redo-treeβ restore ModelFinder and only redo tree search.--undoβ revoke finished run, used when changing some options.--cptime NUMβ minimum checkpoint interval (default: 60 sec and adapt).
π§© PARTITION MODEL π§©
-p FILE|DIRβ NEXUS/RAxML partition file or directory with alignments. Edge-linked proportional partition model.-q FILE|DIRβ like-pbut edge-linked equal partition model.-Q FILE|DIRβ like-pbut edge-unlinked partition model.-S FILE|DIRβ like-pbut separate tree inference.--subsample NUMβ randomly sub-sample partitions (negative for complement).--subsample-seed NUMβ random number seed for--subsample.
πΊοΈ LIKELIHOOD/QUARTET MAPPING πΊοΈ
--lmap NUMβ number of quartets for likelihood mapping analysis.--lmclust FILEβ NEXUS file containing clusters for likelihood mapping.--quartetlhβ print quartet log-likelihoods to.quartetlhfile.
π³ TREE SEARCH ALGORITHM π³
--ninit NUMβ number of initial parsimony trees (default: 100).--ntop NUMβ number of top initial trees (default: 20).--nbest NUMβ number of best trees retained during search (default: 5).-n NUMβ fix number of iterations to stop (default: OFF).--nstop NUMβ number of unsuccessful iterations to stop (default: 100).--perturb NUMβ perturbation strength for randomized NNI (default: 0.5).--radius NUMβ radius for parsimony SPR search (default: 6).--allnniβ perform more thorough NNI search (default: OFF).-g FILEβ (multifurcating) topological constraint tree file.--fastβ fast search to resemble FastTree.--polytomyβ collapse near-zero branches into polytomy.--tree-fixβ fix-ttree (no tree search performed).--treelsβ write locally optimal trees into.treelsfile.--show-lhβ compute tree likelihood without optimisation.--terraceβ check if the tree lies on a phylogenetic terrace.
β‘ ULTRAFAST BOOTSTRAP/JACKKNIFE β‘
-B,--ufboot NUMβ replicates for ultrafast bootstrap (>=1000).-J,--ufjack NUMβ replicates for ultrafast jackknife (>=1000).--jack-prop NUMβ subsampling proportion for jackknife (default: 0.5).--sampling STRINGβGENE|GENESITEresampling for partitions (default:SITE).--boot-treesβ write bootstrap trees to.ufbootfile (default: none).--wbtlβ like--boot-treesbut also writing branch lengths.--nmax NUMβ maximum number of iterations (default: 1000).--nstep NUMβ iterations for UFBoot stopping rule (default: 100).--bcor NUMβ minimum correlation coefficient (default: 0.99).--beps NUMβ RELL epsilon to break tie (default: 0.5).--bnniβ optimize UFBoot trees by NNI on bootstrap alignment.
π NON-PARAMETRIC BOOTSTRAP/JACKKNIFE π
-b,--boot NUMβ replicates for bootstrap + ML tree + consensus tree.-j,--jack NUMβ replicates for jackknife + ML tree + consensus tree.--jack-prop NUMβ subsampling proportion for jackknife (default: 0.5).--bcon NUMβ replicates for bootstrap + consensus tree.--bonly NUMβ replicates for bootstrap only.--tbeβ transfer bootstrap expectation.
πΏ SINGLE BRANCH TEST πΏ
--alrt NUMβ replicates for SH approximate likelihood ratio test.--alrt 0β parametric aLRT test (Anisimova and Gascuel 2006).--abayesβ approximate Bayes test (Anisimova et al. 2011).--lbp NUMβ replicates for fast local bootstrap probabilities.
π MODEL-FINDER π
--use-nn-modelβ use neural network for tree inference.--nn-path-modelβ neural network file for substitution model (onnx format).--nn-path-ratesβ neural network file for alpha value (onnx format).-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 ...+LMβ additionally test Lie Markov models.-m ...+LMRYβ additionally test Lie Markov models with RY symmetry.-m ...+LMWSβ additionally test Lie Markov models with WS symmetry.-m ...+LMMKβ additionally test Lie Markov models with MK symmetry.-m ...+LMSSβ additionally test strand-symmetric models.--mset STRINGβ restrict search to models supported by other programs (raxml,phyml,mrbayes,beast1orbeast2). Ifmrbayesis selected, will output a MrBayes Block File if Data Type is supported.--mset STR,...β comma-separated model list (e.g.-mset WAG,LG,JTT).--msub STRINGβ amino-acid model source (nuclear,mitochondrial,chloroplastorviral).--mfreq STR,...β list of state frequencies.--mrate STR,...β list of rate heterogeneity among sites (e.g.-mrate E,I,G,I+G,Ris used for-m MF).--cmin NUMβ min categories for FreeRate model[+R](default: 2).--cmax NUMβ max categories for FreeRate model[+R](default: 10).--merit AIC|AICc|BICβ Akaike|Bayesian information criterion (default:BIC).--mtreeβ perform full tree search for every model.--madd STR,...β list of mixture models to consider.--mdef FILEβ model definition NEXUS file (see Manual).--modelomaticβ find best codon/protein/DNA models (Whelan et al. 2015).
π§© PARTITION-FINDER π§©
--mergeβ merge partitions to increase model fit.--merge greedy|rcluster|rclusterfβ set merging algorithm (default:rclusterf).--merge-model 1|allβ use only 1 or all models for merging (default: 1).--merge-model STR,...β comma-separated model list for merging.--merge-rate 1|allβ use only 1 or all rate heterogeneity (default: 1).--merge-rate STR,...β comma-separated rate list for merging.--rcluster NUMβ percentage of partition pairs for rcluster algorithm.--rclusterf NUMβ percentage of partition pairs for rclusterf algorithm.--rcluster-max NUMβ max number of partition pairs (default: 10*partitions).
𧬠SUBSTITUTION MODEL π§¬
-m STRINGβ model name string (e.g.GTR+F+I+G). π now here's the tea on whatSTRINGcan be:- 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,FLAVI,Q.LG,Q.pfam,Q.pfam_gb,Q.bird,Q.mammal,Q.insect,Q.plant,Q.yeast - 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_YYwhereXXis empirical andYYis mechanistic model - Morphology/SNP:
MK(default),ORDERED,GTR - Lie Markov DNA:
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(optionally prefixed byRY,WSorMK) - Non-reversible:
STRSYM(strand symmetric model, equiv.WS6.6),NONREV,UNREST(unrestricted model, equiv.12.12)NQ.pfam,NQ.bird,NQ.mammal,NQ.insect,NQ.plant,NQ.yeast - Otherwise: Name of file containing user-model parameters
- DNA:
ποΈ STATE FREQUENCY ποΈ
-m ...+Fβ empirically counted frequencies from alignment.-m ...+FOβ optimized frequencies by maximum-likelihood.-m ...+FQβ equal frequencies.-m ...+FRYβ for DNA,freq(A+G)=1/2=freq(C+T).-m ...+FWSβ for DNA,freq(A+T)=1/2=freq(C+G).-m ...+FMKβ for DNA,freq(A+C)=1/2=freq(G+T).-m ...+Fabcdβ 4-digit constraint on ACGT frequency (e.g.+F1221meansf_A=f_T,f_C=f_G).-m ...+FUβ amino-acid frequencies given protein matrix.-m ...+F1x4β equal NT frequencies over three codon positions.-m ...+F3x4β unequal NT frequencies over three codon positions.
π₯ RATE HETEROGENEITY AMONG SITES π₯
-m ...+Iβ a proportion of invariable sites.-m ...+G[n]β discrete Gamma model with n categories (default n=4).-m ...*G[n]β discrete Gamma model with unlinked model parameters.-m ...+I+G[n]β invariable sites PLUS the Gamma model, with n rate categories, bestie π-m ...+R[n]β FreeRate model with n categories (default n=4) π-m ...*R[n]β FreeRate model, but the model parameters are unlinked π-m ...+I+R[n]β invariable sites + FreeRate model with n categories, the whole combo meal β¨-m ...+Hnβ Heterotachy model with n classes π¦-m ...*Hnβ Heterotachy model with n classes AND unlinked parameters π--alpha-min NUMβ minimum Gamma shape parameter for site rates (default: 0.02)--gamma-medianβ median approximation for+Gsite rates (default is mean, fyi)--rateβ writes the empirical Bayesian site rates into a.ratefile π--mlrateβ writes maximum likelihood site rates into a.mlratefile π
𧬠POLYMORPHISM AWARE MODELS (PoMo) π§¬
-s FILEβ input counts file (go read the manual for this one, bestie)-m ...+Pβ a DNA substitution model (see above β¬οΈ) used together with PoMo-m ...+N<POPSIZE>β virtual population size (default: 9)-m ...+WB|WH|S]β weighted binomial sampling-m ...+WHβ weighted hypergeometric sampling-m ...+Sβ sampled sampling-m ...+G[n]β discrete Gamma rate with n categories (default n=4) π
π COMPLEX MODELS π
-m "MIX{m1,...,mK}"β mixture model with K components-m "FMIX{f1,...fK}"β frequency mixture model with K components--mix-optβ optimize the mixture weights (default: detect) βοΈ-m ...+ASCβ ascertainment bias correction--tree-freq FILEβ input tree used to infer the site frequency model π³--site-freq FILEβ input site frequency model file--freq-maxβ posterior maximum instead of the mean approximation
π₯ TREE TOPOLOGY TEST π₯
--trees FILEβ a set of trees whose log-likelihoods get evaluated--test NUMβ number of replicates for the topology test--test-weightβ perform weighted KH and SH tests--test-auβ approximately unbiased (AU) test (Shimodaira 2002)--au-epsilon NUMβ epsilon for the AU test: if|deltaL| < NUM, the tree is kept no matter what the p-value says (default: 0.001)--sitelhβ write site log-likelihoods into a.sitelhfile
πΈ ANCESTRAL STATE RECONSTRUCTION πΈ
--ancestralβ ancestral state reconstruction via empirical Bayes π°οΈ--asr-min NUMβ minimum probability of an ancestral state (default: equil freq)
π TEST OF SYMMETRY π
--symtestβ perform three tests of symmetry--symtest-onlyβ do--symtestthen exist--symtest-remove-badβ do--symtestand remove the bad partitions ποΈ--symtest-remove-goodβ do--symtestand remove the good partitions ποΈ--symtest-type MAR|INTβ use the MARginal/INTernal test when removing partitions--symtest-pval NUMERβ p-value cutoff (default: 0.05)--symtest-keep-zeroβ keep the NAs in the tests π«Ά
π¦ CONCORDANCE FACTOR ANALYSIS π¦
-t FILEβ reference tree to assign the concordance factor to--gcf FILEβ set of source trees for the gene concordance factor (gCF)--df-treeβ write the discordant trees associated with gDF1--scf NUMβ number of quartets for the site concordance factor (sCF)--scfl NUMβ like--scfbut using likelihood (recommended, fr fr β)-s FILEβ sequence alignment for--scf-p FILE|DIRβ partition file or directory for--scf--cf-verboseβ write CF per tree/locus intocf.stat_tree/_loci--cf-quartetβ write the sCF for all resampled quartets into.cf.quartet
π§ ALISIM: ALIGNMENT SIMULATOR π§
Usage: iqtree --alisim <OUTPUT_PREFIX> [-m MODEL] [-t TREE] ...
--alisim OUTPUT_ALIGNMENTβ activate AliSim and pick the output alignment filename β¨-t TREE_FILEβ set the input tree file name--length LENGTHβ set the length of the root sequence--num-alignments NUMBERβ set the number of output datasets--seqtype STRINGβBIN,DNA,AA,CODON,MORPH{NUM_STATES}(default: auto-detect). For morphological data,0<NUM_STATES<=32π--m MODEL_STRINGβ specify the evolutionary model. See the Manual for way more detail--mdef FILEβ name of a NEXUS model file used to define new models (see Manual)--fundi TAXA_LIST,RHOβ specify a list of taxa, and Rho (the Fundi weight) for the FunDi model--indel <INS>,<DEL>β set the insertion and deletion rate of the indel model, relative to the substitution rate--indel-size <INS_DIS>,<DEL_DIS>β set the insertion and deletion size distributions--sub-level-mixtureβ enable the feature that simulates a substitution-level mixture model--no-unalignedβ stop outputting a file of unaligned sequences when indel models are in play--root-seq FILE,SEQ_NAMEβ specify the root sequence from an alignment-s FILEβ specify the input sequence alignment--no-copy-gapsβ stop copying gaps from the input alignment (default: false)--site-freq <OPTION>β pick the option (MEAN(default), orSAMPLING, orMODEL) to mimic the site-frequencies for mixture models from the input alignment (see Manual)--site-rate <OPTION>β pick the option (MEAN(default), orSAMPLING, orMODEL) to mimic the discrete rate heterogeneity from the input alignment (see Manual)-t RANDOM{MODEL,NUM_TAXA}β specify the model and the number of taxa to generate a random tree π³-rlen MIN MEAN MAXβ specify three numbers: minimum, mean and maximum branch lengths when generating a random tree-p FILEβ NEXUS/RAxML partition file, edge-linked proportional partition model-q FILEβ like-pbut edge-linked equal partition model-Q FILEβ like-pbut edge-unlinked partition model--distribution FILEβ supply a definition file of distributions, which could be used to generate random model parameters--branch-distribution DISβ specify a distribution, from which branch lengths of the input trees are randomly generated and overridden--branch-scale SCALEβ specify a value to scale all branch lengths π--single-outputβ output all alignments into a single file--write-allβ enable outputting internal sequences--seed NUMβ random seed number (default: CPU clock). Be careful: to make AliSim reproducible, users should specify the seed number π±-gzβ enable output compression, but it takes longer running time π-af phy|fastaβ set the output format (default: phylip)
User Manual is available at http://www.iqtree.org/doc/AliSim π
π ANALYSIS WITH GENTRIUS ALGORITHM π
--gentrius FILEβ the file must contain either a single species-tree or a set of subtrees-pr_ab_matrix FILEβ presence-absence matrix of loci coverage-s FILEβ PHYLIP/FASTA/NEXUS/CLUSTAL/MSF alignment file(s)-p FILEβ NEXUS/RAxML partition file-g_stop_t NUMβ stop after NUM species-trees were generated, or use0to turn off this stopping rule. Default: 1MLN trees.-g_stop_i NUMβ stop after NUM intermediate trees were visited, or use0to turn off this stopping rule. Default: 10MLN trees.-g_stop_h NUMβ stop after NUM hours (CPU time), or use0to turn off this stopping rule. Default: 7 days.-g_non_stopβ turn off all stopping rules. Chaos! Go nuts πββοΈ-g_query FILEβ species-trees to test for an identical set of subtrees-g_printβ write all generated species-trees. WARNING: there might be millions of trees! π-g_print_lim NUMβ limit on the number of species-trees to be written-g_print_inducedβ write induced partition subtrees-g_print_mβ write the presence-absence matrix-g_rm_leaves NUMβ invoke reverse analysis for complex datasets
π¦ GENOMIC EPIDEMIOLOGICAL ANALYSIS π¦
--pathogenβ apply the CMAPLE tree search algorithm if sequence divergence is low, otherwise apply the IQ-TREE algorithm--pathogen-forceβ apply the CMAPLE tree search algorithm no matter how divergent the sequences are--alrt <num_rep>β specify the number of replicates to compute SH-aLRT--sprtaβ compute SPRTA (DeMaio et al., 2024) branch supports--sprta-zero-branchβ compute SPRTA supports for zero-length branches--sprta-other-placesβ output alternative SPRs and their SPRTA supports-T <num_thread>β specify the number of threads used for computing branch supports (SH-aLRT or SPRTA) π§΅
β³ TIME TREE RECONSTRUCTION β³
--date FILEβ file containing dates of tips or ancestral nodes--date TAXNAMEβ extract dates from taxon names after the last'|'--date-tip STRINGβ tip dates as a real number orYYYY-MM-DD--date-root STRINGβ root date as a real number orYYYY-MM-DD--date-ci NUMβ number of replicates to compute the confidence interval--clock-sd NUMβ std-dev for the lognormal relaxed clock (default: 0.2)--date-no-outgroupβ exclude the outgroup from the time tree--date-outlier NUMβ Z-score cutoff to remove outlier tips/nodes (e.g.3)--date-options ".."β extra options passing directly to LSD2--dating STRINGβ dating method:LSDfor least square dating (default)
okay bestie that's the whole thing, you're officially a phylogenetic icon now π this was iqtree3-mpi version 3.1.1+dfsg, dropped in July 2026, and yes it's serving IQTREE3-MPI(1) energy the entire time πβ¨ go build that tree fr fr π³π«Ά
π¦ where this page comes from
- ubuntu stonking: iqtree3 3.1.1+dfsg-2build1 Β· GPL-2.0-or-later
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash on : a language model's rewrite of the original, under its terms