IQTREE2-MPI(1)
π IQTREE2-MPI β the manual page, but make it make sense π
okay bestie, this is iqtree2-mpi(1), aka the multiprocessor version of IQ-TREE2 β efficient phylogenetic software by maximum likelihood π§¬β¨. It lives in the User Commands section of the manual, and it's lowkey a whole queen. Sit down, hydrate, and let's read her together π«Ά
πΈ NAME πΈ
iqtree2-mpi β efficient phylogenetic software by maximum likelihood, but make it multiprocessor π
(that's the MPI version, iykyk).
β¨ SYNOPSIS β¨
iqtree [-s ALIGNMENT] [-p PARTITION] [-m MODEL] [-t TREE] ...
yes bestie, that's the whole vibe. Options go in, trees come out π³
π DESCRIPTION π
This is IQ-TREE MPI multicore version 2.4.0 for Linux x86 64-bit, built Jul 28 2026 ποΈπΈ
Developed by Bui Quang Minh, Nguyen Lam Tung, Olga Chernomor, Heiko Schmidt, Dominik Schrempf, Michael Woodhams, Ly Trong Nhan, Thomas Wong β absolute icons, no cap π
π GENERAL OPTIONS: π
-h,--helpβ prints (even more) help usages, like she's not done explaining herself-s FILE[,...,FILE]β yourPHYLIP/FASTA/NEXUS/CLUSTAL/MSFalignment file(s) π-s DIRβ a whole directory of alignment files--seqtype STRINGβ pick your type:BIN,DNA,AA,NT2AA,CODON,MORPH(default: auto-detect, she'll figure it out herself)-t FILE|PARS|RANDβ starting tree (default:99parsimony andBIONJ)-o TAX[,...,TAX]β outgroup taxon (or a whole list of them) for writing.treefile--prefix STRINGβ prefix for all output files (default:aln/partition)--seed NUMβ random seed number, normally used for debugging purposes π--safeβ safe likelihood kernel, so you don't get numerical underflow (she's being careful π₯Ί)--mem NUM[G|M|%]β maximal RAM usage in GB | MB | % π§--runs NUMβ number of independent runs (default:1)-v,--verboseβ verbose mode, prints more messages to screen (she's a talker π¬)-V,--versionβ display version number--quietβ quiet mode, suppresses printing to screen (stdout) π€«-fconst f1,...,fNβ add constant patterns into alignment (N= no. states)--epsilon NUMβ likelihood epsilon for parameter estimate (default0.01)-T NUM|AUTOβ no. cores/threads, orAUTO-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 the tree search--undoβ revoke finished run, used when changing some options--cptime NUMβ minimum checkpoint interval (default60sec and adapt)
π¦ PARTITION MODEL: π¦
-p FILE|DIRβNEXUS/RAxMLpartition file, or a 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βNEXUSfile 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) β she's thorough π-g FILEβ (multifurcating) topological constraint tree file--fastβ fast search to resembleFastTree--polytomyβ collapse near-zero branches into polytomy--tree-fixβ fix the-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: π
-m TESTONLYβ standard model selection (likejModelTest,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)--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 definitionNEXUSfile (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 only1orallmodels for merging (default:1)--merge-model STR,...β comma-separated model list for merging--merge-rate 1|allβ use only1orallrate heterogeneity (default:1)--merge-rate STR,...β comma-separated rate list for merging--rcluster NUMβ percentage of partition pairs forrclusteralgorithm--rclusterf NUMβ percentage of partition pairs forrclusterfalgorithm--rcluster-max NUMβ max number of partition pairs (default:10*partitions)
π SUBSTITUTION MODEL π
π the -m flag herself
-m STRINGβ model name string (e.g.GTR+F+I+G)
and now, bestie, the model menu. it's giving Cheesecake Factory menu. every single one is here π
πΈ DNA:
HKY (default), JC, F81, K2P, K3P, K81uf, TN/TrN, TNef, TIM, TIMef, TVM, TVMef, SYM, GTR, or a 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_YY where XX is empirical and YY is 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 by RY, WS or MK)
π 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
π 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.+F1221means f_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 Gamma model with n categories-m ...+R[n]β FreeRate model with n categories (default n=4)-m ...*R[n]β FreeRate model with unlinked model parameters-m ...+I+R[n]β invariable sites plus FreeRate model with n categories-m ...+Hnβ Heterotachy model with n classes-m ...*Hnβ the heterotachy model, bestie, with n classes and parameters that are unlinked. Like, each class is living its own independent life π¦--alpha-minNUM β smallest Gamma shape parameter allowed for site rates (default:0.02)--gamma-medianβ use the median approximation for+Gsite rates (default is the mean, jsyk)--rateβ write the empirical Bayesian site rates out to a.ratefile π--mlrateβ write the maximum likelihood site rates out to a.mlratefile β¨
𧬠POLYMORPHISM AWARE MODELS (PoMo) π
-sFILE β the input counts file (see the manual, iykyk)-m ...+Pβ a DNA substitution model (the ones from up 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 (defaultn=4)
π COMPLEX MODELS (she's complicated, and that's iconic) π
-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 is detect)-m ...+ASCβ ascertainment bias correction (no cap, this one is required a lot of the time)--tree-freqFILE β an input tree, used to infer the site frequency model--site-freqFILE β an input site frequency model file--freq-maxβ use the posterior maximum instead of the mean approximation
π₯ TREE TOPOLOGY TEST π₯
--treesFILE β the set of trees whose log-likelihoods you wanna evaluate--testNUM β how many replicates for the topology test--test-weightβ do the weighted KH and SH tests--test-auβ the approximately unbiased (AU) test (Shimodaira 2002) π«--sitelhβ write the site log-likelihoods out to a.sitelhfile
π―οΈ ANCESTRAL STATE RECONSTRUCTION π―οΈ
--ancestralβ ancestral state reconstruction via empirical Bayes--asr-minNUM β minimum probability of an ancestral state (default:equil freq)
βοΈ TEST OF SYMMETRY βοΈ
--symtestβ run three tests of symmetry--symtest-onlyβ do--symtestthen exist (their spelling, not mine π she's doing her best)--symtest-remove-badβ do--symtestand yeet the bad partitions ποΈ--symtest-remove-goodβ do--symtestand yeet the good partitions (yes, good, that's what it says π)--symtest-type MAR|INTβ use the MARginal or INTernal test when removing partitions--symtest-pvalNUMER β the p-value cutoff (default:0.05)--symtest-keep-zeroβ keep theNAsin the tests π
πΈ CONCORDANCE FACTOR ANALYSIS πΈ
-tFILE β the reference tree used to assign the concordance factor--gcfFILE β the set of source trees for the gene concordance factor (gCF)--df-treeβ write the discordant trees associated withgDF1--scfNUM β the number of quartets for the site concordance factor (sCF)--scflNUM β like--scfbut using likelihood (recommended, bestie, this is the one π₯)-sFILE β the sequence alignment for--scf-pFILE|DIR β a partition file or directory for--scf--cf-verboseβ write CF per tree/locus tocf.stat_tree/_loci--cf-quartetβ write sCF for all the resampled quartets to.cf.quartet
π§ ALISIM: ALIGNMENT SIMULATOR π§
Usage, verbatim, mwah:
iqtree --alisim <OUTPUT_PREFIX> [-m MODEL] [-t TREE] ...
--alisimOUTPUT_ALIGNMENT β activate AliSim and name the output alignment filename β¨-tTREE_FILE β set the input tree file name--lengthLENGTH β set the length of the root sequence--num-alignmentsNUMBER β set how many output datasets you want--seqtypeSTRING βBIN,DNA,AA,CODON,MORPH{NUM_STATES}(default:auto-detect). For morphological data,0<NUM_STATES<=32--mMODEL_STRING β specify the evolutionary model. See the Manual for more detail--mdefFILE β the name of a NEXUS model file that defines new models (see Manual)--fundiTAXA_LIST,RHO β a list of taxa plus 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 simulating a substitution-level mixture model--no-unalignedβ disable outputting a file of unaligned sequences when using indel models--root-seqFILE,SEQ_NAME β specify the root sequence from an alignment-sFILE β specify the input sequence alignment--no-copy-gapsβ disable copying gaps from the input alignment (default:false)--site-freq<OPTION>β specify the option (MEAN(default), orSAMPLING, orMODEL) to mimic the site-frequencies for mixture models from the input alignment (see Manual)--site-rate<OPTION>β specify the option (MEAN(default), orSAMPLING, orMODEL) to mimic the discrete rate heterogeneity from the input alignment (see Manual)-tRANDOM{MODEL,NUM_TAXA} β specify the model and number of taxa to generate a random tree π³-rlenMIN MEAN MAX β three numbers pls: minimum, mean and maximum branch lengths when generating a random tree-pFILE β NEXUS/RAxML partition file, edge-linked proportional partition model-qFILE β like-pbut edge-linked equal partition model-QFILE β like-pbut edge-unlinked partition model--distributionFILE β supply a definition file of distributions, which could be used to generate random model parameters--branch-distributionDIS β specify a distribution from which the branch lengths of the input trees are randomly generated and overridden--branch-scaleSCALE β specify a value to scale all branch lengths--single-outputβ output all alignments into a single file--write-allβ enable outputting internal sequences π--seedNUM β 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)
π The User Manual is available at http://www.iqtree.org/doc/alisim
π ANALYSIS WITH GENTRIUS ALGORITHM π
--gentriusFILE β the file must contain either a single species-tree or a set of subtrees-pr_ab_matrixFILE β presence-absence matrix of loci coverage-sFILE βPHYLIP/FASTA/NEXUS/CLUSTAL/MSFalignment file(s)-pFILE β NEXUS/RAxML partition file-g_stop_tNUM β stop after NUM species-trees were generated, or use0to turn this stopping rule off. Default:1MLNtrees.-g_stop_iNUM β stop after NUM intermediate trees were visited, or use0to turn this stopping rule off. Default:10MLNtrees.-g_stop_hNUM β stop after NUM hours (CPU time), or use0to turn this stopping rule off. Default:7 days.-g_non_stopβ turn off all the stopping rules (chaos mode, go off queen π)-g_queryFILE β 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_limNUM β 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_leavesNUM β 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 regardless of sequence divergence, no matter what πͺ
β³ TIME TREE RECONSTRUCTION β³
--dateFILE β file containing dates of tips or ancestral nodes--dateTAXNAME β extract dates from taxon names after the last'|'--date-tipSTRING β tip dates as a real number orYYYY-MM-DD--date-rootSTRING β root date as a real number orYYYY-MM-DD--date-ciNUM β number of replicates to compute the confidence interval--clock-sdNUM β std-dev for the lognormal relaxed clock (default:0.2)--date-no-outgroupβ exclude the outgroup from the time tree--date-outlierNUM β Z-score cutoff to remove outlier tips/nodes (e.g.3)--date-options".."β extra options passed directly toLSD2--datingSTRING β dating method:LSDfor least square dating (default)
and that's the whole page, bestie, you made it to the end π go run some phylogenies, you absolute genius π«Ά
this was iqtree2-mpi version 2.4.0+dfsg, dated July 2026, and yes the man page itself is called IQTREE2-MPI(1) β love that for us ππΈβ¨
π¦ where this page comes from
- ubuntu stonking: iqtree 2.4.0+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