IQTREE2-OMP(1)
πΈ IQTREE2-OMP(1) πΈ
okay so bestie, this is the manual page from the IQTREE2-OMP section of the manual, and it's all about a little program called iqtree2-omp. it's long, so we're doing it in parts β this is part 1, don't panic π
β¨ NAME β¨
iqtree2-omp β efficient phylogenetic software by maximum likelihood (multiprocessor version) π¦
β¨ SYNOPSIS β¨
iqtree [-s ALIGNMENT] [-p PARTITION] [-m MODEL] [-t TREE] ...
(yes the binary in the synopsis is literally iqtree, that's just how it is, iykyk π
)
β¨ DESCRIPTION β¨
IQ-TREE 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 legends, no cap.
π GENERAL OPTIONS:
-h,--helpβ prints (more) help usages. your friendly lil "what do i even do" button π«Ά-s FILE[,...,FILE]β PHYLIP/FASTA/NEXUS/CLUSTAL/MSF alignment file(s). feed it your data bestie πΈ-s DIRβ a whole directory of alignment files π--seqtype STRINGβ the sequence type:BIN,DNA,AA,NT2AA,CODON,MORPH(default: auto-detect π)-t FILE|PARS|RANDβ the starting tree (default:99parsimony and BIONJ π±)-o TAX[,...,TAX]β outgroup taxon (or a whole list of them) for writing.treefile--prefix STRINGβ prefix for all your output files (default:aln/partition) π--seed NUMβ random seed number, normally used for debugging purpose π²--safeβ safe likelihood kernel to avoid numerical underflow (safety first diva π )--mem NUM[G|M|%]β maximal RAM usage in GB | MB | %--runs NUMβ number of indepedent runs (default:1)-v,--verboseβ verbose mode, prints more messages to screen π’-V,--versionβ display version number--quietβ quiet mode, suppresses printing to screen (stdout) π€«-fconst f1,...,fNβ add constant patterns into the 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 (default:60sec and adapt)
π PARTITION MODEL:
-p FILE|DIRβ NEXUS/RAxML partition 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 number = the 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 the.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 the-ttree (no tree search performed) π--treelsβ write locally optimal trees into the.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 (must be>=1000) π-J,--ufjack NUMβ replicates for ultrafast jackknife (must be>=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 the.ufbootfile (default: none)--wbtlβ like--boot-treesbut also writing branch lengths--nmax NUMβ maximum number of iterations (default:1000)--nstep NUMβ iterations for the 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 (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 the 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 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 the rcluster algorithm--rclusterf NUMβ percentage of partition pairs for the 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)
here's the big menu of model names, bestie, take your time π
- 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, - Protein:
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 a file containing user-model parameters π
π STATE FREQUENCY:
-m ...+Fβ empirically counted frequencies from the 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 withncategories (defaultn=4)-m ...*G[n]β discrete Gamma model with unlinked model parameters-m ...+I+G[n]β invariable sites plus Gamma model withncategories-m ...+R[n]β FreeRate model withncategories (defaultn=4)-m ...*R[n]β FreeRate model with unlinked model parameters-m ...+I+R[n]β invariable sites plus FreeRate model withncategories-m ...+Hnβ heterotachy model withnclasses-m ...*Hnβ the heterotachy model, bestie, with n classes and unlinked parameters β every class gets to keep its own settings, no sharing, no borrowing, main character energy only π--alpha-min NUMβ the minimum Gamma shape parameter for site rates (default:0.02)--gamma-medianβ use the median approximation for+Gsite rates (default is mean, jsyk π)--rateβ write out empirical Bayesian site rates into a.ratefile β¨--mlrateβ write out maximum likelihood site rates into a.mlratefile π₯
𧬠POLYMORPHISM AWARE MODELS (PoMo) π§¬
-s FILEβ the input counts file (go read the manual for this one, it's giving special)-m ...+Pβ the DNA substitution model (the one from up above) that you use alongside 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β use the posterior maximum instead of the mean approximation
π TREE TOPOLOGY TEST π
--trees FILEβ a whole set of trees whose log-likelihoods get evaluated--test NUMβ how many replicates for the topology test--test-weightβ perform weighted KH and SH tests--test-auβ the approximately unbiased (AU) test (Shimodaira 2002)--sitelhβ write site log-likelihoods into a.sitelhfile
π ANCESTRAL STATE RECONSTRUCTION π
--ancestralβ ancestral state reconstruction, empirical Bayes style--asr-min NUMβ minimum probability of an ancestral state (default: equil freq)
π¦ TEST OF SYMMETRY π¦
--symtestβ run three tests of symmetry--symtest-onlyβ do--symtestthen exist (yes, "exist", that's literally what it says and honestly, iconic)--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 you're 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 factors--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 π)-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 resampled quartets to.cf.quartet
π§ ALISIM: ALIGNMENT SIMULATOR π§
Usage, periodt:
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 how long the root sequence is--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β specify the evolutionary model. See the Manual for more detail--mdef FILEβ name of a NEXUS model file to define new models (see Manual)--fundi TAXA_LIST,RHOβ specify a list of taxa, and Rho (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 to simulate substitution-level mixture model--no-unalignedβ disable 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β 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)-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 the AliSim reproducible, users should specify the seed number-gzβ enable output compression but taking 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 π
--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 afterNUMspecies-trees were generated, or use0to turn off this stopping rule. Default: 1MLN trees-g_stop_i NUMβ stop afterNUMintermediate trees were visited, or use0to turn off this stopping rule. Default: 10MLN trees-g_stop_h NUMβ stop afterNUMhours (CPU time), or use0to turn off this stopping rule. Default: 7 days-g_non_stopβ turn off all 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β limit on the number of species-trees to be written-g_print_inducedβ write induced partition subtrees-g_print_mβ write 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 regardless of sequence divergence
π 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 confidence interval--clock-sd NUMβ std-dev for lognormal relaxed clock (default:0.2)--date-no-outgroupβ exclude outgroup from 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)
And that's the whole thing, bestie π«Ά that was iqtree2-omp version 2.4.0+dfsg, dated July 2026, straight from the IQTREE2-OMP(1) page. Go forth and phylogeny β¨πΈ
π¦ 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