IQTREE3-OMP(1)
iqtree3-omp(1) β User Commands. Multiprocessor version, bestie! π
π NAME π
iqtree3-omp β efficient phylogenetic software by maximum likelihood (multiprocessor version). No cap, she's iconic. π
π¦ SYNOPSIS π¦
iqtree [-s ALIGNMENT] [-p PARTITION] [-m MODEL] [-t TREE] ...
πΈ DESCRIPTION πΈ
IQ-TREE version 3.1.1 for Linux x86 64-bit built Jul 28 2026. 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. β¨
β¨ 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 independent runs (default: 1) π-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)
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
π 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 the Gamma model, with n rate categories bestie π-m ...+R[n]β FreeRate model, n categories (and yes,n=4if you don't say) β¨-m ...*R[n]β FreeRate but the model parameters are unlinked, they're doing their own thing fr-m ...+I+R[n]β invariable sites plus FreeRate, n categories, the whole combo meal π-m ...+Hnβ heterotachy model with n classes-m ...*Hnβ heterotachy with n classes and unlinked parameters--alpha-min NUMβ the minimum Gamma shape parameter for site rates (default:0.02)--gamma-medianβ use the median approximation for+Gsite rates (the default is mean, btw)--rateβ write the empirical Bayesian site rates out to a.ratefile π--mlrateβ write maximum likelihood site rates out to a.mlratefile
𧬠POLYMORPHISM AWARE MODELS (PoMo) π§¬
-s FILEβ the input counts file (peep the manual for this one)-m ...+Pβ a DNA substitution model (any of the ones above) used together with PoMo-m ...+N<POPSIZE>β set the 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 π
-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 to evaluate log-likelihoods on--test NUMβ number of replicates for the topology test--test-weightβ perform weighted KH and SH tests--test-auβ the approximately unbiased (AU) test (Shimodaira 2002)--au-epsilon NUMβ epsilon for the AU test: if|deltaL| < NUM, the tree gets 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 for an ancestral state (default: equil freq)
π TEST OF SYMMETRY π
--symtestβ perform three tests of symmetry--symtest-onlyβ run--symtestand then just exit--symtest-remove-badβ run--symtestand remove the bad partitions--symtest-remove-goodβ run--symtestand remove the good partitions--symtest-type MAR|INTβ use the MARginal or INTernal test when removing partitions--symtest-pval NUMERβ the p-value cutoff (default:0.05)--symtest-keep-zeroβ keep the NAs in the tests π
π¦ CONCORDANCE FACTOR ANALYSIS π¦
-t FILEβ reference tree used to assign concordance factor--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β same as--scfbut it uses likelihood (this one's recommended, no cap)-s FILEβ sequence alignment for--scf-p FILE|DIRβ 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 into.cf.quartet
π ALISIM: ALIGNMENT SIMULATOR π
Usage:
iqtree --alisim <OUTPUT_PREFIX> [-m MODEL] [-t TREE] ...
--alisim OUTPUT_ALIGNMENTβ switch AliSim ON and set 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 how many output datasets you want--seqtype STRINGβBIN,DNA,AA,CODON,MORPH{NUM_STATES}(default: auto-detect). For morphological data,0<NUM_STATES<=32--m MODEL_STRINGβ pick the evolutionary model. The Manual has more detail on this one--mdef FILEβ name of a NEXUS model file that defines new models (see the 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β turn on the feature that simulates a substitution-level mixture model--no-unalignedβ stop outputting a file of unaligned sequences when you're using indel models--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>β specify the option (MEAN(default), orSAMPLING, orMODEL) to mimic the site-frequencies for mixture models from the input alignment (see the Manual)--site-rate <OPTION>β specify the option (MEAN(default), orSAMPLING, orMODEL) to mimic the discrete rate heterogeneity from the input alignment (see the Manual)-t RANDOM{MODEL,NUM_TAXA}β specify the model and the number of taxa to generate a random tree-rlen MIN MEAN MAXβ give 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β hand over a definition file of distributions, which can be used to generate random model parameters--branch-distribution DISβ specify a distribution that the branch lengths of the input trees get randomly generated from and overridden with--branch-scale SCALEβ specify a value to scale all branch lengths--single-outputβ output all alignments into one single file--write-allβ turn on outputting internal sequences--seed NUMβ random seed number (default: CPU clock). Be careful β if you want your AliSim runs to be reproducible, you gotta specify the seed number, bestie-gzβ enable output compression, but it takes longer to run-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 π
--gentrius FILEβ the file has to 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 this stopping rule off. Default: 1MLN trees.-g_stop_i NUMβ stop after NUM intermediate trees were visited, or use0to turn this stopping rule off. Default: 10MLN trees.-g_stop_h NUMβ 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 π³-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β a limit on how many species-trees get written-g_print_inducedβ write the induced partition subtrees-g_print_mβ write the presence-absence matrix-g_rm_leaves NUMβ invoke reverse analysis for the complex datasets
π¦ GENOMIC EPIDEMIOLOGICAL ANALYSIS π¦
--pathogenβ apply the CMAPLE tree search algorithm if the sequence divergence is low, otherwise apply the IQ-TREE algorithm--pathogen-forceβ apply the CMAPLE tree search algorithm regardless of 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 the alternative SPRs and their SPRTA supports-T <num_thread>β specify the number of threads used to compute branch supports (SH-aLRT or SPRTA)
β TIME TREE RECONSTRUCTION β
--date FILEβ file containing dates of tips or ancestral nodes--date TAXNAMEβ pull the dates out of the 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 for computing the confidence interval--clock-sd NUMβ standard deviation for the lognormal relaxed clock (default:0.2)--date-no-outgroupβ exclude the outgroup from the time tree--date-outlier NUMβ Z-score cutoff used to remove outlier tips/nodes (e.g.3)--date-options ".."β extra options that get passed straight through to LSD2--dating STRINGβ dating method:LSDfor least square dating (this is the default)
and that's a wrap on the iqtree3-omp page, bestie π this one ships as iqtree3-omp 3.1.1+dfsg, dated July 2026, and it lives at IQTREE3-OMP(1) β go read the manual and slay your phylogenies πβ¨
π¦ 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