Mercurial > repos > leomrtns > iqtree
annotate iqtree.xml @ 3:a396ccf6079c draft default tip
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
| author | leomrtns |
|---|---|
| date | Wed, 18 Dec 2019 13:46:22 +0000 |
| parents | 87daf702e477 |
| children |
| rev | line source |
|---|---|
| 0 | 1 <tool id="iqtree" name="IQ-TREE" version="@TOOL_VERSION@" > |
| 2 <description>Phylogenetic tree construction from multiple sequences and analyses via maximum likelihood</description> | |
| 3 <macros> | |
| 4 <import>iqtree_macros.xml</import> | |
| 5 </macros> | |
| 6 <expand macro="requirements" /> | |
| 7 <command detect_errors='exit_code' ><![CDATA[ | |
| 8 iqtree | |
| 9 -pre PREF | |
| 10 -nt \${GALAXY_SLOTS:-10} | |
| 11 -redo | |
| 12 | |
| 13 ## file | |
| 14 #if $general_options.s | |
| 15 -s '$general_options.s' | |
| 16 #end if | |
| 17 | |
| 18 ## file | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
19 #if str($initial_tree.mode.tree) == "file" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
20 -t '$initial_tree.mode.t' |
| 0 | 21 $tree_parameters.constructing_consensus.con |
| 22 $tree_parameters.constructing_consensus.net | |
| 23 #if str($tree_parameters.constructing_consensus.bi) != '' | |
| 24 -bi '$tree_parameters.constructing_consensus.bi' | |
| 25 #end if | |
| 26 | |
| 27 ## file | |
| 28 #if $tree_parameters.constructing_consensus.sup | |
| 29 -sup '$tree_parameters.constructing_consensus.sup' | |
| 30 #end if | |
| 31 | |
| 32 #if str($tree_parameters.constructing_consensus.suptag) != '' | |
| 33 -suptag '$tree_parameters.constructing_consensus.suptag' | |
| 34 #end if | |
| 35 | |
| 36 ## file | |
| 37 #if $tree_parameters.computing_robinson_foulds.rf | |
| 38 -rf '$tree_parameters.computing_robinson_foulds.rf' | |
| 39 $tree_parameters.computing_robinson_foulds.rf_all | |
| 40 $tree_parameters.computing_robinson_foulds.rf_adj | |
| 41 #end if | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
42 #else if str($initial_tree.mode.tree) == "fixed" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
43 -te '$initial_tree.mode.te' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
44 #else if str($initial_tree.mode.tree) == "default" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
45 #if str($initial_tree.mode.ninit) != '' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
46 -ninit '$initial_tree.mode.ninit' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
47 #end if |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
48 #if str($initial_tree.mode.sprrad) != '' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
49 -sprrad '$initial_tree.mode.sprrad' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
50 #end if |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
51 #else if str($initial_tree.mode.tree) == "fast" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
52 -fast -nbest 1 |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
53 #else if str($initial_tree.mode.tree) == "pars" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
54 -t PARS -ninit 2 -nbest 1 |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
55 #else |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
56 -t BIONJ -ninit 2 -nbest 1 |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
57 #end if |
| 0 | 58 |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
59 #if str($initial_tree.mode.tree) in ["fixed", "default", "file"]: |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
60 #if str($tree_parameters.tree_search.nbest) != '' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
61 -nbest '$tree_parameters.tree_search.nbest' |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
62 #end if |
| 0 | 63 #end if |
| 64 | |
| 65 #if str($general_options.st) != '' | |
| 66 -st '$general_options.st' | |
| 67 #end if | |
| 68 | |
| 69 #if str($general_options.seed) != '' | |
| 70 -seed '$general_options.seed' | |
| 71 #end if | |
| 72 | |
| 73 $general_options.keep_ident | |
| 74 $general_options.safe | |
| 75 | |
| 76 #if str($likelihood_mapping.lmap) != '' | |
| 77 -lmap '$likelihood_mapping.lmap' | |
| 78 #end if | |
| 79 | |
| 80 ## file | |
| 81 #if $likelihood_mapping.lmclust | |
| 82 -lmclust '$likelihood_mapping.lmclust' | |
| 83 #end if | |
| 84 | |
| 85 $likelihood_mapping.wql | |
| 86 | |
| 87 #if str($modelling_parameters.automatic_model.cond_model.m) != '' | |
| 88 -m '$modelling_parameters.automatic_model.cond_model.m' | |
| 89 #end if | |
| 90 | |
| 91 #if str($modelling_parameters.automatic_model.rcluster) != '' | |
| 92 -rcluster '$modelling_parameters.automatic_model.rcluster' | |
| 93 #end if | |
| 94 | |
| 95 #if str($modelling_parameters.automatic_model.mset) != '' | |
| 96 -mset '$modelling_parameters.automatic_model.mset' | |
| 97 #end if | |
| 98 | |
| 99 #if str($modelling_parameters.automatic_model.msub) != '' | |
| 100 -msub '$modelling_parameters.automatic_model.msub' | |
| 101 #end if | |
| 102 | |
| 103 #if str($modelling_parameters.automatic_model.mfreq) != '' | |
| 104 -mfreq '$modelling_parameters.automatic_model.mfreq' | |
| 105 #end if | |
| 106 | |
| 107 #if str($modelling_parameters.automatic_model.mrate) != '' | |
| 108 -mrate '$modelling_parameters.automatic_model.mrate' | |
| 109 #end if | |
| 110 | |
| 111 #if str($modelling_parameters.automatic_model.cmin) != '' | |
| 112 -cmin '$modelling_parameters.automatic_model.cmin' | |
| 113 #end if | |
| 114 | |
| 115 #if str($modelling_parameters.automatic_model.cmax) != '' | |
| 116 -cmax '$modelling_parameters.automatic_model.cmax' | |
| 117 #end if | |
| 118 | |
| 119 #if str($modelling_parameters.automatic_model.merit) != '' | |
| 120 -merit '$modelling_parameters.automatic_model.merit' | |
| 121 #end if | |
| 122 | |
| 123 $modelling_parameters.automatic_model.mtree | |
| 124 | |
| 125 #if str($modelling_parameters.automatic_model.madd) != '' | |
| 126 -madd '$modelling_parameters.automatic_model.madd' | |
| 127 #end if | |
| 128 | |
| 129 ## file | |
| 130 #if $modelling_parameters.automatic_model.mdef | |
| 131 -mdef '$modelling_parameters.automatic_model.mdef' | |
| 132 #end if | |
| 133 | |
| 134 $modelling_parameters.specifying_substitution.mwopt | |
| 135 | |
| 136 #if str($modelling_parameters.rate_heterogeneity.a) != '' | |
| 137 -a '$modelling_parameters.rate_heterogeneity.a' | |
| 138 #end if | |
| 139 | |
| 140 $modelling_parameters.rate_heterogeneity.gmedian | |
| 141 | |
| 142 #if str($modelling_parameters.rate_heterogeneity.i) != '' | |
| 143 -i '$modelling_parameters.rate_heterogeneity.i' | |
| 144 #end if | |
| 145 | |
| 146 $modelling_parameters.rate_heterogeneity.opt_gamma_inv | |
| 147 $modelling_parameters.rate_heterogeneity.wsr | |
| 148 | |
| 149 ## file | |
| 150 #if $modelling_parameters.partition_model.q | |
| 151 -q '$modelling_parameters.partition_model.q' | |
| 152 #end if | |
| 153 | |
| 154 $modelling_parameters.partition_model.spp | |
| 155 | |
| 156 ## file | |
| 157 #if $modelling_parameters.partition_model.sp | |
| 158 -sp '$modelling_parameters.partition_model.sp' | |
| 159 #end if | |
| 160 | |
| 161 ## file | |
| 162 #if $modelling_parameters.site_specific_frequency.ft | |
| 163 -ft '$modelling_parameters.site_specific_frequency.ft' | |
| 164 #end if | |
| 165 | |
| 166 #if str($modelling_parameters.site_specific_frequency.fs) != '' | |
| 167 -fs '$modelling_parameters.site_specific_frequency.fs' | |
| 168 #end if | |
| 169 | |
| 170 $modelling_parameters.site_specific_frequency.fmax | |
| 171 | |
| 2 | 172 |
| 0 | 173 #if str($tree_parameters.tree_search.ntop) != '' |
| 174 -ntop '$tree_parameters.tree_search.ntop' | |
| 175 #end if | |
| 176 | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
177 |
| 0 | 178 |
| 179 #if str($tree_parameters.tree_search.nstop) != '' | |
| 180 -nstop '$tree_parameters.tree_search.nstop' | |
| 181 #end if | |
| 182 | |
| 183 #if str($tree_parameters.tree_search.n) != '' | |
| 184 -n '$tree_parameters.tree_search.n' | |
| 185 #end if | |
| 186 | |
| 187 #if str($tree_parameters.tree_search.pers) != '' | |
| 188 -pers '$tree_parameters.tree_search.pers' | |
| 189 #end if | |
| 190 | |
| 191 $tree_parameters.tree_search.allnni | |
| 192 $tree_parameters.tree_search.djc | |
| 193 | |
| 194 ## file | |
| 195 #if $tree_parameters.tree_search.g | |
| 196 -g '$tree_parameters.tree_search.g' | |
| 197 #end if | |
| 198 | |
| 199 #if str($tree_parameters.single_branch.alrt) != '' | |
| 200 -alrt '$tree_parameters.single_branch.alrt' | |
| 201 #end if | |
| 202 | |
| 203 $tree_parameters.single_branch.abayes | |
| 204 | |
| 205 #if str($tree_parameters.single_branch.lbp) != '' | |
| 206 -lbp '$tree_parameters.single_branch.lbp' | |
| 207 #end if | |
| 208 | |
| 209 ## file | |
| 210 #if $tree_parameters.tree_topology.z | |
| 211 -z '$tree_parameters.tree_topology.z' | |
| 212 #end if | |
| 213 | |
| 214 #if str($tree_parameters.tree_topology.zb) != '' | |
| 215 -zb '$tree_parameters.tree_topology.zb' | |
| 216 #end if | |
| 217 | |
| 218 $tree_parameters.tree_topology.zw | |
| 219 $tree_parameters.tree_topology.au | |
| 220 | |
| 221 #if str($tree_parameters.constructing_consensus.minsup) != '' | |
| 222 -minsup '$tree_parameters.constructing_consensus.minsup' | |
| 223 #end if | |
| 224 | |
| 225 | |
| 226 #if str($tree_parameters.generating_random.r) != '' | |
| 227 -r '$tree_parameters.generating_random.r' | |
| 228 #end if | |
| 229 | |
| 230 $tree_parameters.generating_random.ru | |
| 231 $tree_parameters.generating_random.rcat | |
| 232 $tree_parameters.generating_random.rbal | |
| 233 $tree_parameters.generating_random.rcsg | |
| 234 | |
| 235 #if str($tree_parameters.generating_random.rlen) != '' | |
| 236 -rlen '$tree_parameters.generating_random.rlen' | |
| 237 #end if | |
| 238 | |
| 239 #if str($bootstrap_parameters.ultrafast_bootstrap.bb) != '' | |
| 240 -bb '$bootstrap_parameters.ultrafast_bootstrap.bb' | |
| 241 #end if | |
| 242 | |
| 243 $bootstrap_parameters.ultrafast_bootstrap.wbt | |
| 244 $bootstrap_parameters.ultrafast_bootstrap.wbtl | |
| 245 | |
| 246 #if str($bootstrap_parameters.ultrafast_bootstrap.nm) != '' | |
| 247 -nm '$bootstrap_parameters.ultrafast_bootstrap.nm' | |
| 248 #end if | |
| 249 | |
| 250 #if str($bootstrap_parameters.ultrafast_bootstrap.bcor) != '' | |
| 251 -bcor '$bootstrap_parameters.ultrafast_bootstrap.bcor' | |
| 252 #end if | |
| 253 | |
| 254 #if str($bootstrap_parameters.ultrafast_bootstrap.nstep) != '' | |
| 255 -nstep '$bootstrap_parameters.ultrafast_bootstrap.nstep' | |
| 256 #end if | |
| 257 | |
| 258 #if str($bootstrap_parameters.ultrafast_bootstrap.beps) != '' | |
| 259 -beps '$bootstrap_parameters.ultrafast_bootstrap.beps' | |
| 260 #end if | |
| 261 | |
| 262 #if str($bootstrap_parameters.ultrafast_bootstrap.bspec) != '' | |
| 263 -bspec '$bootstrap_parameters.ultrafast_bootstrap.bspec' | |
| 264 #end if | |
| 265 | |
| 266 $bootstrap_parameters.ultrafast_bootstrap.bnni | |
| 267 | |
| 268 #if str($bootstrap_parameters.nonparametric_bootstrap.b) != '' | |
| 269 -b '$bootstrap_parameters.nonparametric_bootstrap.b' | |
| 270 #end if | |
| 271 | |
| 272 $bootstrap_parameters.nonparametric_bootstrap.bc | |
| 273 $bootstrap_parameters.nonparametric_bootstrap.bo | |
| 274 | |
| 275 #if str($miscellaneous_options.fconst) != '' | |
| 276 -fconst '$miscellaneous_options.fconst' | |
| 277 #end if | |
| 278 | |
| 279 ]]> | |
| 280 </command> | |
| 281 <inputs> | |
| 282 <section name="general_options" expanded="True" title="General options"> | |
| 283 <param argument="-s" type="data" format="txt" optional="true" label="Specify input alignment file in PHYLIP, FASTA, NEXUS, CLUSTAL or MSF format."/> | |
| 284 <param argument="-st" type="select" label="Specify sequence type as either of DNA, AA, BIN, MORPH, CODON or NT2AA for DNA, amino-acid, binary, morphological, codon or DNA-to-AA-translated sequences"> | |
| 285 <option value="DNA">DNA</option> | |
| 286 <option value="AA">AA</option> | |
| 287 <option value="BIN">BIN</option> | |
| 288 <option value="MORPH">MORPH</option> | |
| 289 <option value="CODON">CODON</option> | |
| 290 <option value="NT2AA">NT2AA</option> | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
291 <help><![CDATA[ |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
292 Note that -st CODON is always necessary when using codon models and you also need to specify a genetic code like this if differed from the standard genetic code. |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
293 <br/><i>-st</i> NT2AA tells IQ-TREE to translate protein-coding DNA into AA sequences and then subsequent analysis will work on the AA sequences. |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
294 You can also use a genetic code like -st NT2AA5 for the Invertebrate Mitochondrial Code (see genetic code table). |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
295 ]]> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
296 </help> |
| 0 | 297 </param> |
| 298 <param argument="-seed" type="integer" optional="true" label="Specify a random number seed to reproduce a previous run (leave blank to randomize)"/> | |
| 299 <param argument="-keep_ident" type="boolean" truevalue="-keep-ident" falsevalue="" checked="false" label="Keep identical sequences in the alignment" help="By default: IQ-TREE will remove them during the analysis and add them in the end."/> | |
| 300 <param argument="-safe" type="boolean" truevalue="-safe" falsevalue="" checked="false" label="Turn on safe numerical mode to avoid numerical underflow for large data sets with many sequences (typically in the order of thousands)" help="This mode is automatically turned on when having more than 2000 sequences."/> | |
| 301 </section> | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
302 |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
303 <section name="initial_tree" expanded="True" title="Input and initial trees"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
304 <conditional name="mode"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
305 <param name="tree" type="select" label="Initial tree search"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
306 <option value="default" selected="true">parsimony using PLL library (default) </option> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
307 <option value="fast">fasttree emulation (option "-fast")</option> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
308 <option value="pars">fast parsimony, using IQTREE internal algorithm (option "-t PARS")</option> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
309 <option value="bionj">bioNJ tree (option "-t BIONJ")</option> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
310 <option value="fixed">fixed tree, no optimisation (option "-te"). Used in e.g. ancestral reconstruction</option> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
311 <option value="file">initial trees from file (option "-t"). Also used for consensus and for tree distance calculations</option> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
312 <help><![CDATA[ |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
313 The options for the initial tree estimation are parsimony (default), fast parsimony, fasttree emulation, bioNJ, fixed (no optimisation), and from file. <br/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
314 The default parsimony estimation is through the PLL library, which is linear on the number of <i>sites</i> and not <i>patterns</i> (so slower than other algorithms). |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
315 Notice that this library is also used to populate the "initial parsimony trees" set ("<i>-ninit</i>" below), so keep that in mind. <br/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
316 |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
317 <br/>According to the documentation (see http://www.iqtree.org/doc/Command-Reference#tree-search-parameters), the option "-fast" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
318 "Turn on the fast tree search mode, where IQ-TREE will just construct two starting trees: maximum parsimony and BIONJ, which are then optimized |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
319 by nearest neighbor interchange (NNI)." |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
320 If you chose fast parsimony, fasttree, or bioNJ, it is assumed that the initial set of trees is two ("-ninit 2") since doing otherwise would defeat the purpose |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
321 of avoiding the default PLL.<br/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
322 <br/>You can also select a tree file as input to IQ_TREE. In most cases this tree file is used as |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
323 initial state for optimisation (or as a fixed topology). But in some cases it is used for consensus tree calculation or Robinson-Foulds distances. |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
324 There are a few reasons why you should provide trees to IQ-TREE: |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
325 <ul> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
326 <li>Through a <b>constrained tree search</b>, where you give a (potentially small, multifurcating) tree which should be respected when searching for optimal trees. |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
327 This is given through option <i>"-g"</i>. You use this option when you "know" or want to test a specific grouping of taxa. (Section "Tree search parameters")</li> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
328 <li>Input trees for <b>tree fit comparison</b>, as in section "Tree topology tests" below. IQ-TREE may still do model optimisation, initial tree search, etc. since |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
329 it compares the user-provided trees assuming a common model. </li> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
330 <li>Suggest the <b>initial trees</b> (option "<i>-t</i>"). From these, IQ-TREE will start exploring trees (to speed up the tree search in case you already have |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
331 good candidate trees). This same command (i.e. same tree file) can also be used for calculating consensus trees or distances between trees. The other options will |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
332 inform IQ-TREE how it should use this file (here, the default is to use it as starting tree). </li> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
333 <li>By <b>fixing the tree</b> (option "<i>-te</i>"). This means that IQ-TREE will not try any optimisation with other topologies, and can be used e.g. in ancestral reconstruction. </li> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
334 </ul> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
335 ]]></help> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
336 </param> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
337 <when value="default"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
338 <param argument="-ninit" type="integer" value="100" optional="true" label="Specify number of initial parsimony trees (default=100)" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
339 help="This value must be strictly higher than the number of best trees to maintain (i.e. higher than 5 under default settings)."/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
340 <param argument="-sprrad" type="integer" value="6" optional="true" label="Specify SPR radius for the initial parsimony tree search"/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
341 </when> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
342 <when value="file"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
343 <param argument="-t" type="data" format="nhx" optional="false" label="Specify a file containing the starting tree for tree search" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
344 help="also used by consensus or distance calculations, if you set the appropriate variables in 'Tree Parameters' below"/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
345 </when> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
346 <when value="fixed"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
347 <param argument="-te" type="data" format="nhx" optional="false" label="file with fixed user tree" help="No tree search is performed and IQ-TREE |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
348 computes the log-likelihood of the fixed user tree."/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
349 </when> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
350 <when value="fast"/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
351 <when value="pars"/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
352 <when value="bionj"/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
353 </conditional> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
354 </section> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
355 |
| 0 | 356 <section name="likelihood_mapping" expanded="False" title="Likelihood mapping analysis"> |
| 357 <param argument="-lmap" type="integer" optional="true" label="Specify the number of quartets to be randomly drawn" help="If you specify -lmap ALL, all unique quartets will be drawn, instead."/> | |
| 358 <param argument="-lmclust" type="data" format="txt" optional="true" label="Specify a NEXUS file containing taxon clusters (see below for example) for quartet mapping analysis."/> | |
| 359 <param argument="-wql" type="boolean" truevalue="-wql" falsevalue="" checked="false" label="Write quartet log-likelihoods into .lmap.quartetlh file (typically not needed)."/> | |
| 360 </section> | |
| 361 <section name="modelling_parameters" title="Modelling Parameters"> | |
| 362 <section name="automatic_model" expanded="False" title="Automatic model selection"> | |
| 363 <conditional name="cond_model" > | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
364 <param name="opt_custommodel" type="boolean" checked="false" label="Use Custom Model"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
365 <help><![CDATA[ |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
366 With this option, you can restrict the analysis to a fixed model instead of letting iq-tree select the best one (slow model selection). |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
367 Suggestions include "HKY+G" (good compromise speed-performance) or "GTR+G+I" (parameter-rich model). |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
368 <br>See http://www.iqtree.org/doc/Substitution-Models for a full list]]> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
369 </help></param> |
| 0 | 370 <when value="true"> |
| 371 <param argument="-m" type="text" label="Model"> | |
| 372 <sanitizer> | |
| 1 | 373 <valid initial="string.ascii_uppercase"> |
| 374 <add value="+" /> | |
| 375 </valid> | |
| 0 | 376 </sanitizer> |
| 377 </param> | |
| 378 </when> | |
| 379 <when value="false"> | |
| 380 <param argument="-m" type="select" label="Perform standard model selection like jModelTest (for DNA) and ProtTest (for protein)" > | |
| 381 <help><![CDATA[ | |
| 382 <b>Note:</b> Only <code>TEST</code> and custom models that do not rely on <code>MF</code> can be used in conjunction with bootstrap parameters (-b)<br/><br/> | |
| 383 IQ-TREE also works for codon, binary and morphogical data. | |
| 384 <table> | |
| 385 <tr> | |
| 386 <td><i>TESTONLY</i></td> | |
| 387 <td>Perform standard model selection like jModelTest (for DNA) and ProtTest (for protein). Moreover, IQ-TREE also works for codon, binary and morphogical data.</td> | |
| 388 </tr> | |
| 389 <tr> | |
| 390 <td><i>TEST</i></td> | |
| 391 <td>Like -m TESTONLY but immediately followed by tree reconstruction using the best-fit model found. So this performs both model selection and tree inference within a single run.</td> | |
| 392 </tr> | |
| 393 <tr> | |
| 394 <td><i>TESTNEWONLY</i></td> | |
| 395 <td>Perform an extended model selection that additionally includes FreeRate model compared with -m TESTONLY. Recommended as replacement for -m TESTONLY. Note that LG4X is a FreeRate model, but by default is not included because it is also a protein mixture model. To include it, use -madd option (see table below).</td> | |
| 396 </tr> | |
| 397 <tr> | |
| 398 <td><i>TESTNEW</i></td> | |
| 399 <td>Like -m MF but immediately followed by tree reconstruction using the best-fit model found.</td> | |
| 400 </tr> | |
| 401 <tr> | |
| 402 <td><i>TESTMERGEONLY</i></td> | |
| 403 <td>Select best-fit partitioning scheme like PartitionFinder.</td> | |
| 404 </tr> | |
| 405 <tr> | |
| 406 <td><i>TESTMERGE</i></td> | |
| 407 <td>Like -m TESTMERGEONLY but immediately followed by tree reconstruction using the best partitioning scheme found.</td> | |
| 408 </tr> | |
| 409 <tr> | |
| 410 <td><i>TESTNEWMERGEONLY</i></td> | |
| 411 <td>Like -m TESTMERGEONLY but additionally includes FreeRate model.</td> | |
| 412 </tr> | |
| 413 <tr> | |
| 414 <td><i>TESTNEWMERGE</i></td> | |
| 415 <td>Like -m MF+MERGE but immediately followed by tree reconstruction using the best partitioning scheme found.</td> | |
| 416 </tr> | |
| 417 </table> | |
| 418 ]]> | |
| 419 </help> | |
| 420 <option value="" selected="true">(None)</option> | |
| 421 <option value="TESTONLY">TESTONLY</option> | |
| 422 <option value="TEST">TEST</option> | |
| 423 <option value="TESTNEWONLY">TESTNEWONLY</option> | |
| 424 <option value="TESTNEW">TESTNEW</option> | |
| 425 <option value="TESTMERGEONLY">TESTMERGEONLY</option> | |
| 426 <option value="TESTMERGE">TESTMERGE</option> | |
| 427 <option value="TESTNEWMERGEONLY">TESTNEWMERGEONLY</option> | |
| 428 <option value="TESTNEWMERGE">TESTNEWMERGE</option> | |
| 429 </param> | |
| 430 </when> | |
| 431 </conditional> | |
| 432 <param argument="-rcluster" type="integer" optional="true" label="Specify the percentage for the relaxed clustering algorithm (Lanfear et al., 2014)" help="This is similar to --rcluster-percent option of PartitionFinder. For example, with -rcluster 10 only the top 10% partition schemes are considered to save computations." /> | |
| 433 <param argument="-mset" type="text" optional="true" label="Specify the name of a program (raxml, phyml or mrbayes) to restrict to only those models supported by the specified program" help="Alternatively, one can specify a comma-separated list of base models. For example, -mset WAG,LG,JTT will restrict model selection to WAG, LG, and JTT instead of all 18 AA models to save computations."> | |
| 434 <expand macro="sanitize_query" | |
| 435 validinitial="string.ascii_uppercase,string.punctuation" /> | |
| 436 </param> | |
| 437 <param argument="-msub" type="select" label="Specify either nuclear, mitochondrial, chloroplast or viral to restrict to those AA models designed for specified source." help=""> | |
| 438 <option value="nuclear">nuclear</option> | |
| 439 <option value="mitochondrial">mitochondrial</option> | |
| 440 <option value="chloroplast">chloroplast</option> | |
| 441 <option value="viral">viral</option> | |
| 442 </param> | |
| 443 <param argument="-mfreq" type="text" optional="true" label="Specify a comma-separated list of frequency types for model selection" > | |
| 444 <expand macro="sanitize_query" | |
| 445 validinitial="string.ascii_uppercase,string.digits,x,string.punctuation" /> | |
| 446 </param> | |
| 447 <param argument="-mrate" type="text" optional="true" label="Specify a comma-separated list of rate heterogeneity types for model selection"> | |
| 448 <expand macro="sanitize_query" | |
| 449 validinitial="string.ascii_uppercase,string.punctuation" /> | |
| 450 </param> | |
| 451 <param argument="-cmin" type="integer" value="2" optional="true" label="Specify minimum number of categories for FreeRate model"/> | |
| 452 <param argument="-cmax" type="integer" value="10" optional="true" label="Specify maximum number of categories for FreeRate model"/> | |
| 453 <param argument="-merit" type="select" label="Specify either AIC, AICc or BIC for the optimality criterion to apply for new procedure" help=""> | |
| 454 <option value="AIC">AIC</option> | |
| 455 <option value="AICc">AICc</option> | |
| 456 <option value="BIC">BIC</option> | |
| 457 </param> | |
| 458 <param argument="-mtree" type="boolean" truevalue="-mtree" falsevalue="" checked="false" label="Turn on full tree search for each model considered, to obtain more accurate result"/> | |
| 459 <param argument="-madd" type="text" optional="true" label="Specify a comma-separated list of mixture models to additionally consider for model selection" help="For example, -madd LG4M,LG4X to additionally include these two protein mixture models."> | |
| 460 <expand macro="sanitize_query" | |
| 461 validinitial="string.ascii_uppercase,string.digits,string.punctuation" /> | |
| 462 </param> | |
| 463 <param argument="-mdef" type="data" format="txt" optional="true" label="Specify a NEXUS model file to define new models."/> | |
| 464 </section> | |
| 465 <section name="specifying_substitution" expanded="False" title="Specifying substitution models"> | |
| 466 <param argument="-mwopt" type="boolean" truevalue="-mwopt" falsevalue="" checked="false" label="Turn on optimizing weights of mixture models" help="Note that for models like LG+C20+F+G this mode is automatically turned on, but not for LG+C20+G."/> | |
| 467 </section> | |
| 468 <section name="rate_heterogeneity" expanded="False" title="Rate heterogeneity"> | |
| 469 <param argument="-a" type="text" optional="true" label="Specify the Gamma shape parameter (default: estimate)"> | |
| 470 <expand macro="sanitize_query" /> | |
| 471 </param> | |
| 472 <param argument="-gmedian" type="boolean" truevalue="-gmedian" falsevalue="" checked="false" label="Perform the median approximation for Gamma rate heterogeneity instead of the default mean approximation (Yang, 1994)"/> | |
| 473 <param argument="-i" type="text" optional="true" label="Specify the proportion of invariable sites (default: estimate)"> | |
| 474 <expand macro="sanitize_query" /> | |
| 475 </param> | |
| 476 <param argument="--opt_gamma_inv" type="boolean" truevalue="--opt-gamma-inv" falsevalue="" checked="false" label="Perform more thorough estimation for +I+G model parameters"/> | |
| 477 <param argument="-wsr" type="boolean" truevalue="-wsr" falsevalue="" checked="false" label="Write per-site rates to .rate file"/> | |
| 478 </section> | |
| 479 <section name="partition_model" expanded="False" title="Partition model options"> | |
| 480 <param argument="-q" type="data" format="txt" optional="true" label="Specify partition file for edge-equal partition model" help="That means, all partitions share the same set of branch lengths (like -q option of RAxML)."/> | |
| 481 <param argument="-spp" type="boolean" truevalue="-spp" falsevalue="" checked="false" label="Like -q but allowing partitions to have different evolutionary speeds (edge-proportional partition model)."/> | |
| 482 <param argument="-sp" type="data" format="txt" optional="true" label="Specify partition file for edge-unlinked partition model" help="That means, each partition has its own set of branch lengths (like -M option of RAxML). This is the most parameter-rich partition model to accomodate heterotachy."/> | |
| 483 </section> | |
| 484 <section name="site_specific_frequency" expanded="False" title="Site-specific frequency model options"> | |
| 485 <param argument="-ft" type="data" format="nhx" optional="true" label="Specify a guide tree to infer site frequency profiles."/> | |
| 486 <param argument="-fs" type="float" optional="true" label="Specify a site frequency file, e.g" help="the .sitefreq file obtained from -ft run. This will save memory used for the first phase of the analysis."/> | |
| 487 <param argument="-fmax" type="boolean" truevalue="-fmax" falsevalue="" checked="false" label="Switch to posterior maximum mode for obtaining site-specific profiles" help="Default: posterior mean."/> | |
| 488 </section> | |
| 489 </section> | |
| 490 <section name="tree_parameters" title="Tree Parameters"> | |
| 491 <section name="tree_search" expanded="False" title="Tree search parameters"> | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
492 <param argument="-nstop" type="integer" value="100" optional="true" label="Specify number of unsuccessful iterations to stop |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
493 (-nstop and -pers are the most relevant when increasing exploration space, according to IQ_TREE reference manual)"/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
494 <param argument="-pers" type="float" value="0.5" optional="true" label="Specify perturbation strength (between 0 and 1) for randomized NNI |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
495 (-nstop and -pers are the most relevant when increasing exploration space, according to IQ_TREE reference manual)"/> |
| 2 | 496 <param argument="-n" type="integer" optional="true" label="Skip subsequent tree search (n = 0), useful when you only want to assess the phylogenetic information of the alignment. Otherwise specify number of iterations to stop (this option overrides -nstop criterion)." /> |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
497 <param argument="-nbest" type="integer" value="5" optional="true" label="Specify number of trees in the candidate set to maintain during ML tree search (default=5)" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
498 help="This value must be strictly smaller than the number of initial pasimony trees (i.e. less than 100 for default settings)"/> |
| 2 | 499 <param argument="-ntop" type="integer" value="20" optional="true" label="Specify number of top initial parsimony trees to optimize with ML nearest neighbor interchange (NNI) search to initialize the candidate set"/> |
| 500 <param argument="-allnni" type="boolean" truevalue="-allnni" falsevalue="" checked="false" label="Turn on more thorough and slower NNI search"/> | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
501 <param argument="-djc" type="boolean" truevalue="-djc" falsevalue="" checked="false" label="Avoid computing ML pairwise distances and BIONJ tree (when generating tree set)" |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
502 help="This is independent on how the initial tree/trees were found, and is added to set of initial trees (since it was estimated using likelihood parameters) "/> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
503 <param argument="-g" type="data" format="txt" optional="true" label="Specify a topological constraint tree file in NEWICK format" help="The constraint tree can be a multifurcating tree and need not to include all taxa. (see https://github.com/Cibiv/IQ-TREE/wiki/Advanced-Tutorial#constrained-tree-search)"/> |
| 0 | 504 </section> |
| 505 <section name="single_branch" expanded="False" title="Single branch tests"> | |
| 506 <param argument="-alrt" type="integer" optional="true" label="Specify number of replicates (>=1000) to perform SH-like approximate likelihood ratio test (SH-aLRT) (Guindon et al., 2010)" help="If number of replicates is set to 0 (-alrt 0), then the parametric aLRT test (Anisimova and Gascuel 2006) is performed, instead of SH-aLRT."/> | |
| 507 <param argument="-abayes" type="boolean" truevalue="-abayes" falsevalue="" checked="false" label="Perform approximate Bayes test (Anisimova et al., 2011)."/> | |
| 508 <param argument="-lbp" type="integer" optional="true" label="Specify number of replicates (>=1000) to perform fast local bootstrap probability method (Adachi and Hasegawa, 1996)."/> | |
| 509 </section> | |
| 510 <section name="tree_topology" expanded="False" title="Tree topology tests"> | |
|
3
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
511 <param argument="-z" type="data" format="txt" optional="true" label="Specify a file containing a set of trees"> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
512 <help><![CDATA[ |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
513 IQ-TREE will compute and compare the log-likelihoods of all trees in this file. It may still search for an optimal tree, in order to find best model and parameters |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
514 (unless you fixed those as well). It is strongly suggested to also perform some tests below (like RELL, AU, etc.). See more info at |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
515 https://github.com/Cibiv/IQ-TREE/wiki/Advanced-Tutorial#tree-topology-tests |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
516 ]]> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
517 </help> |
|
a396ccf6079c
"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
leomrtns
parents:
2
diff
changeset
|
518 </param> |
| 0 | 519 <param argument="-zb" type="integer" optional="true" label="Specify the number of RELL (Kishino et al., 1990) replicates (>=1000) to perform several tree topology tests for all trees passed via -z" help="The tests include bootstrap proportion (BP), KH test (Kishino and Hasegawa, 1989), SH test (Shimodaira and Hasegawa, 1999) and expected likelihood weights (ELW) (Strimmer and Rambaut, 2002)."/> |
| 520 <param argument="-zw" type="boolean" truevalue="-zw" falsevalue="" checked="false" label="Used together with -zb to additionally perform the weighted-KH and weighted-SH tests."/> | |
| 521 <param argument="-au" type="boolean" truevalue="-au" falsevalue="" checked="false" label="Used together with -zb to additionally perform the approximately unbiased (AU) test (Shimodaira, 2002)" help="Note that you have to specify the number of replicates for the AU test via -zb."/> | |
| 522 </section> | |
| 523 <section name="constructing_consensus" expanded="False" title="Constructing consensus tree"> | |
| 524 <param argument="-con" type="boolean" truevalue="-con" falsevalue="" checked="false" label="Compute consensus tree of the trees passed via -t" help="Resulting consensus tree is written to .contree file."/> | |
| 525 <param argument="-net" type="boolean" truevalue="-net" falsevalue="" checked="false" label="Compute consensus network of the trees passed via -t" help="Resulting consensus network is written to .nex file."/> | |
| 526 <param argument="-minsup" type="float" value="0.0" optional="true" label="Specify a minimum threshold (between 0 and 1) to keep branches in the consensus tree"/> | |
| 527 <param argument="-bi" type="integer" optional="true" label="Specify a burn-in, which is the number of beginning trees passed via -t to discard before consensus construction" help="This is useful e.g. when summarizing trees from MrBayes analysis."/> | |
| 528 <param argument="-sup" type="data" format="nhx" optional="true" label="Specify an input “target” tree file" help="That means, support values are first extracted from the trees passed via -t, and then mapped onto the target tree. Resulting tree with assigned support values is written to .suptree file. This option is useful to map and compare support values from different approaches onto a single tree."/> | |
| 529 <param argument="-suptag" type="text" optional="true" label="Specify name of a node in -sup target tree" help="The corresponding node of .suptree will then be assigned with IDs of trees where this node appears. Special option -suptag ALL will assign such IDs for all nodes of the target tree."/> | |
| 530 </section> | |
| 531 <section name="computing_robinson_foulds" expanded="False" title="Computing Robinson-Foulds distance"> | |
| 532 <param argument="-rf" type="data" format="nhx" optional="true" label="Specify a second set of trees" help="IQ-TREE computes all pairwise RF distances between two tree sets passed via -t and -rf"/> | |
| 533 <param argument="-rf_all" type="boolean" truevalue="-rf_all" falsevalue="" checked="false" label="Compute all-to-all RF distances between all trees passed via -t"/> | |
| 534 <param argument="-rf_adj" type="boolean" truevalue="-rf_adj" falsevalue="" checked="false" label="Compute RF distances between adjacent trees passed via -t"/> | |
| 535 </section> | |
| 536 <section name="generating_random" expanded="False" title="Generating random trees"> | |
| 537 <param argument="-r" type="integer" optional="true" label="Specify number of taxa" help="IQ-TREE will create a random tree under Yule-Harding model with specified number of taxa"/> | |
| 538 <param argument="-ru" type="boolean" truevalue="-ru" falsevalue="" checked="false" label="Like -r, but a random tree is created under uniform model."/> | |
| 539 <param argument="-rcat" type="boolean" truevalue="-rcat" falsevalue="" checked="false" label="Like -r, but a random caterpillar tree is created."/> | |
| 540 <param argument="-rbal" type="boolean" truevalue="-rbal" falsevalue="" checked="false" label="Like -r, but a random balanced tree is created."/> | |
| 541 <param argument="-rcsg" type="boolean" truevalue="-rcsg" falsevalue="" checked="false" label="Like -r, bur a random circular split network is created."/> | |
| 542 <param argument="-rlen" type="text" optional="true" label="Specify three numbers: minimum, mean and maximum branch lengths of the random tree"> | |
| 543 <sanitizer> | |
| 544 <valid initial="string.digits,string.whitespace" /> | |
| 545 </sanitizer> | |
| 546 </param> | |
| 547 </section> | |
| 548 </section> | |
| 549 <section name="bootstrap_parameters" title="Bootstrap Parameters"> | |
| 550 <section name="ultrafast_bootstrap" expanded="False" title="Ultrafast bootstrap parameters"> | |
| 551 <param argument="-bb" type="integer" optional="true" label="Specify number of bootstrap replicates (>=1000)."/> | |
| 552 <param argument="-wbt" type="boolean" truevalue="-wbt" falsevalue="" checked="false" label="Turn on writing bootstrap trees to .ufboot file"/> | |
| 553 <param argument="-wbtl" type="boolean" truevalue="-wbtl" falsevalue="" checked="false" label="Like -wbt but bootstrap trees written with branch lengths"/> | |
| 554 <param argument="-nm" type="integer" value="1000" optional="true" label="Specify maximum number of iterations to stop"/> | |
| 555 <param argument="-bcor" type="float" value="0.99" optional="true" label="Specify minimum correlation coefficient for UFBoot convergence criterion"/> | |
| 556 <param argument="-nstep" type="integer" value="100" optional="true" label="Specify iteration interval checking for UFBoot convergence"/> | |
| 557 <param argument="-beps" type="float" value="0.5" optional="true" label="Specify a small epsilon to break tie in RELL evaluation for bootstrap trees"/> | |
| 558 <param argument="-bspec" type="text" optional="true" label="Specify the resampling strategies for partitioned analysis" help="By default, IQ-TREE resamples alignment sites within partitions. With -bspec GENE IQ-TREE will resample partitions. With -bspec GENESITE IQ-TREE will resample partitions and then resample sites within resampled partitions (Gadagkar et al., 2005)."> | |
| 559 <sanitizer> | |
| 560 <valid initial="string.ascii_uppercase" /> | |
| 561 </sanitizer> | |
| 562 </param> | |
| 563 <param argument="-bnni" type="boolean" truevalue="-bnni" falsevalue="" checked="false" label="Perform an additional step to further optimize UFBoot trees by nearest neighbor interchange (NNI) based directly on bootstrap alignments" help="This option is recommended in the presence of severe model violations. It increases computing time by 2-fold but reduces the risk of overestimating branch supports due to severe model violations. Introduced in IQ-TREE 1.6."/> | |
| 564 </section> | |
| 565 <section name="nonparametric_bootstrap" expanded="False" title="Nonparametric bootstrap"> | |
| 566 <param argument="-b" type="integer" optional="true" label="Specify number of bootstrap replicates (recommended >=100)" help="This will perform both bootstrap and analysis on original alignment and provide a consensus tree."/> | |
| 567 <param argument="-bc" type="boolean" truevalue="-bc" falsevalue="" checked="false" label="Like -b but omit analysis on original alignment."/> | |
| 568 <param argument="-bo" type="boolean" truevalue="-bo" falsevalue="" checked="false" label="Like -b but only perform bootstrap analysis (no analysis on original alignment and no consensus tree)."/> | |
| 569 </section> | |
| 570 </section> | |
| 571 <section name="miscellaneous_options" expanded="False" title="Miscellaneous options"> | |
| 572 <param argument="-fconst" type="text" optional="true" label="Specify a list of comma-separated integer numbers" help="The number of entries should be equal to the number of states in the model (e.g. 4 for DNA and 20 for protein). IQ-TREE will then add a number of constant sites accordingly. For example, -fconst 10,20,15,40 will add 10 constant sites of all A, 20 constant sites of all C, 15 constant sites of all G and 40 constant sites of all T into the alignment."> | |
| 573 <expand macro="sanitize_query" | |
| 574 validinitial="string.digits,string.punctuation" /> | |
| 575 </param> | |
| 576 </section> | |
| 577 </inputs> | |
| 578 <outputs> | |
| 579 <data name='bionj' format='nhx' from_work_dir="*.bionj" label="${tool.name} on ${on_string}: BIONJ Tree" /> | |
| 580 <data name='treefile' format='nhx' from_work_dir="*.treefile" label="${tool.name} on ${on_string}: MaxLikelihood Tree" /> | |
| 581 <data name='contree' format='nhx' from_work_dir="*.contree" label="${tool.name} on ${on_string}: Consensus Tree" /> | |
| 582 <data name='mldist' format='mldist' from_work_dir="*.mldist" label="${tool.name} on ${on_string}: MaxLikelihood Distance Matrix"/> | |
| 583 <data name='splits.nex' format='nex' from_work_dir="*.splits.nex" label="${tool.name} on ${on_string}: Occurence Frequencies in Bootstrap Trees" /> | |
| 584 <data name='iqtree' format='iqtree' from_work_dir="*.iqtree" label="${tool.name} on ${on_string}: Report and Final Tree" /> | |
| 585 </outputs> | |
| 586 <tests> | |
| 587 <test> | |
| 588 <param name="seed" value="1257" /> | |
| 589 <param name="st" value="AA" /> | |
| 590 <param name="s" value="example.phy" /> | |
| 591 <param name="m" value="TESTNEW" /> | |
| 592 <param name="msub" value="nuclear" /> | |
| 593 <param name="madd" value="LG4M,LG4X" /> | |
| 594 <param name="merit" value="AICc" /> | |
| 595 <param name="bb" value="1000" /> | |
| 596 <output name='bionj'> | |
| 597 <assert_contents> | |
| 598 <has_text_matching expression=".*Human.*Whale.*" /> | |
| 599 </assert_contents> | |
| 600 </output> | |
| 601 <output name='iqtree'> | |
| 602 <assert_contents> | |
| 603 <has_text_matching expression="VT\+F\+R3(\s+((-|\d|\.)+))+" /> | |
| 604 </assert_contents> | |
| 605 </output> | |
| 606 <output name='mldist'> | |
| 607 <assert_contents> | |
| 608 <has_line_matching expression="^Frog(\s+((\d|\.)+))+\s+$" /> | |
| 609 <has_line_matching expression="^Whale(\s+((\d|\.)+))+\s+$" /> | |
| 610 <has_line_matching expression="^Cow(\s+((\d|\.)+))+\s+$" /> | |
| 611 </assert_contents> | |
| 612 </output> | |
| 613 <output name='treefile'> | |
| 614 <assert_contents> | |
| 615 <has_line_matching expression="\(LngfishAu:(\d|\..)+,\(LngfishSA:(\d.)+,.*" /> | |
| 616 </assert_contents> | |
| 617 </output> | |
| 618 <output name='contree' > | |
| 619 <assert_contents> | |
| 620 <has_line_matching | |
| 621 expression="\(LngfishAu:(\d|\..)+,\(LngfishSA:(\d.)+,.*\)\d+:(\d|\.)+,.*" /> | |
| 622 </assert_contents> | |
| 623 </output> | |
| 624 <output name='splits.nex'> | |
| 625 <assert_contents> | |
| 626 <has_line line="BEGIN Splits;" /> | |
| 627 <has_line line="END; [Splits]" /> | |
| 628 </assert_contents> | |
| 629 </output> | |
| 630 </test> | |
| 631 <test> | |
| 632 <!-- bootstrap sans model --> | |
| 633 <param name="seed" value="1257" /> | |
| 634 <param name="st" value="AA" /> | |
| 635 <param name="s" value="example.phy" /> | |
| 636 <!-- <param name="m" value="TESTONLY" /> --> | |
| 637 <param name="msub" value="nuclear" /> | |
| 638 <param name="cmin" value="2" /> | |
| 639 <param name="cmax" value="10" /> | |
| 640 <param name="merit" value="AICc" /> | |
| 641 <param name="madd" value="LG4M,LG4X" /> | |
| 642 <param name="ninit" value="100" /> | |
| 643 <param name="ntop" value="20" /> | |
| 644 <param name="nbest" value="5" /> | |
| 645 <param name="nstop" value="100" /> | |
| 646 <param name="sprrad" value="6" /> | |
| 647 <param name="pers" value="0.5" /> | |
| 648 <param name="minsup" value="0.0" /> | |
| 649 <param name="nm" value="1000" /> | |
| 650 <param name="bcor" value="0.99" /> | |
| 651 <param name="nstep" value="100" /> | |
| 652 <param name="beps" value="0.5" /> | |
| 653 <param name="b" value="2" /> | |
| 654 <output name='iqtree'> | |
| 655 <assert_contents> | |
| 656 <has_text_matching expression="VT\+F\+I(\s+((-|\d|\.)+))+" /> | |
| 657 </assert_contents> | |
| 658 </output> | |
| 659 <output name='treefile'> | |
| 660 <assert_contents> | |
| 661 <has_line_matching expression="\(LngfishAu:(\d|\..)+,\(LngfishSA:(\d.)+,.*" /> | |
| 662 </assert_contents> | |
| 663 </output> | |
| 664 </test> | |
| 665 <test> | |
| 666 <!-- model sans bootstrap --> | |
| 667 <param name="seed" value="1257" /> | |
| 668 <param name="st" value="AA" /> | |
| 669 <param name="s" value="example.phy" /> | |
| 670 <param name="m" value="TESTONLY" /> | |
| 671 <param name="msub" value="nuclear" /> | |
| 672 <param name="cmin" value="2" /> | |
| 673 <param name="cmax" value="10" /> | |
| 674 <param name="merit" value="AICc" /> | |
| 675 <param name="madd" value="LG4M,LG4X" /> | |
| 676 <param name="ninit" value="100" /> | |
| 677 <param name="ntop" value="20" /> | |
| 678 <param name="nbest" value="5" /> | |
| 679 <param name="nstop" value="100" /> | |
| 680 <param name="sprrad" value="6" /> | |
| 681 <param name="pers" value="0.5" /> | |
| 682 <param name="minsup" value="0.0" /> | |
| 683 <param name="nm" value="1000" /> | |
| 684 <param name="bcor" value="0.99" /> | |
| 685 <param name="nstep" value="100" /> | |
| 686 <param name="beps" value="0.5" /> | |
| 687 <!-- <param name="b" value="100" /> --> | |
| 688 <output name='iqtree'> | |
| 689 <assert_contents> | |
| 690 <has_text_matching expression="VT\+F\+I(\s+((-|\d|\.)+))+" /> | |
| 691 </assert_contents> | |
| 692 </output> | |
| 693 <output name='treefile'> | |
| 694 <assert_contents> | |
| 695 <has_line_matching expression="\(LngfishAu:(\d|\..)+,\(LngfishSA:(\d.)+,.*" /> | |
| 696 </assert_contents> | |
| 697 </output> | |
| 698 </test> | |
| 699 <test> | |
| 700 <param name='s' value='example.phy' /> | |
| 701 <output name='iqtree'> | |
| 702 <assert_contents> | |
| 703 <has_text_matching expression=".*Human.*Whale.*" /> | |
| 704 </assert_contents> | |
| 705 </output> | |
| 706 <output name='treefile'> | |
| 707 <assert_contents> | |
| 708 <has_text_matching expression="\(LngfishAu:(\d|\..)+,\(LngfishSA:(\d.)+,.*" /> | |
| 709 </assert_contents> | |
| 710 </output> | |
| 711 </test> | |
| 712 </tests> | |
| 713 <help><![CDATA[ | |
| 714 IQ-TREE | |
| 715 ======= | |
| 716 | |
| 717 The full documentation can be found here_. | |
| 718 | |
| 719 .. _here : http://www.iqtree.org/doc/ | |
| 720 | |
| 721 General Tutorial | |
| 722 =================== | |
| 723 | |
| 724 Input data | |
| 725 ---------- | |
| 726 | |
| 727 IQ-TREE takes as input a *multiple sequence alignment* and will reconstruct an evolutionary tree that is best explained by the input data. The input alignment can be in various common formats. For example the PHYLIP_ format which may look like: | |
| 728 | |
| 729 .. _PHYLIP: http://evolution.genetics.washington.edu/phylip/doc/sequence.html | |
| 730 | |
| 731 :: | |
| 732 | |
| 733 7 28 | |
| 734 Frog AAATTTGGTCCTGTGATTCAGCAGTGAT | |
| 735 Turtle CTTCCACACCCCAGGACTCAGCAGTGAT | |
| 736 Bird CTACCACACCCCAGGACTCAGCAGTAAT | |
| 737 Human CTACCACACCCCAGGAAACAGCAGTGAT | |
| 738 Cow CTACCACACCCCAGGAAACAGCAGTGAC | |
| 739 Whale CTACCACGCCCCAGGACACAGCAGTGAT | |
| 740 Mouse CTACCACACCCCAGGACTCAGCAGTGAT | |
| 741 | |
| 742 This tiny alignment contains 7 DNA sequences from several animals with the sequence length of 28 nucleotides. IQ-TREE also supports other file formats such as FASTA, NEXUS, CLUSTALW. The FASTA file for the above example may look like this: | |
| 743 | |
| 744 :: | |
| 745 | |
| 746 >Frog | |
| 747 AAATTTGGTCCTGTGATTCAGCAGTGAT | |
| 748 >Turtle | |
| 749 CTTCCACACCCCAGGACTCAGCAGTGAT | |
| 750 >Bird | |
| 751 CTACCACACCCCAGGACTCAGCAGTAAT | |
| 752 >Human | |
| 753 CTACCACACCCCAGGAAACAGCAGTGAT | |
| 754 >Cow | |
| 755 CTACCACACCCCAGGAAACAGCAGTGAC | |
| 756 >Whale | |
| 757 CTACCACGCCCCAGGACACAGCAGTGAT | |
| 758 >Mouse | |
| 759 CTACCACACCCCAGGACTCAGCAGTGAT | |
| 760 | |
| 761 **NOTE**: If you have raw sequences, you need to first apply alignment programs like MAFFT_ or ClustalW_ to align the sequences, before feeding them into IQ-TREE. | |
| 762 | |
| 763 .. _MAFFT: http://mafft.cbrc.jp/alignment/software/ | |
| 764 .. _ClustalW: http://www.clustal.org | |
| 765 | |
| 766 | |
| 767 Running example | |
| 768 --------------------- | |
| 769 | |
| 770 From the download_ there is an example alignment called `example.phy` | |
| 771 in PHYLIP format. This example contains parts of the mitochondrial DNA sequences of several animals (Source: `Phylogenetic Handbook`_) | |
| 772 | |
| 773 .. _`Phylogenetic Handbook` : http://www.kuleuven.be/aidslab/phylogenybook/home.html)). | |
| 774 .. _download: http://www.iqtree.org/#download | |
| 775 | |
| 776 You can now start to reconstruct a maximum-likelihood tree | |
| 777 from this alignment by entering (assuming that you are now in the same folder with `example.phy`): | |
| 778 | |
| 779 ::iqtree -s example.phy -st AA -seed 9999 -m TESTNEW -msub nuclear -madd LG4M,LG4x -merit AICc -bb 2000 | |
| 780 | |
| 781 * **-s** is the option to specify the name of the alignment file | |
| 782 | |
| 783 * **-st** specifies the sequence type as amino-acid | |
| 784 | |
| 785 * **-seed** ensures that the output files remain the same for subsequent runs | |
| 786 | |
| 787 * **-m** sets the modelling parameter for standard model selection | |
| 788 | |
| 789 * **-msub** determines the type sub-modelling | |
| 790 | |
| 791 * **-madd** provides an additional selection mixed models | |
| 792 | |
| 793 * **-merit** asserts the type of optimality criterion | |
| 794 | |
| 795 * **-bb** defines the number of replicates | |
| 796 | |
| 797 | |
| 798 Each of these parameters are available under the relevant sub-sections in the main tool interface. | |
| 799 | |
| 800 At the end of the run IQ-TREE will write several output files including: | |
| 801 | |
| 802 * **example.phy.iqtree**: the main report file that is self-readable. You should look at this file to see the computational results. It also contains a textual representation of the final tree (see below). | |
| 803 | |
| 804 * **example.phy.treefile**: the ML tree in NEWICK format, which can be visualized by any supported tree viewer programs like FigTree or iTOL. | |
| 805 | |
| 806 **NOTE**: Starting with version 1.5.4, with this simple command IQ-TREE will by default perform ModelFinder to find the best-fit substitution model and then infer a phylogenetic tree using the selected model. | |
| 807 | |
| 808 For this example data the resulting maximum-likelihood tree may look like this (extracted from **.iqtree** file): | |
| 809 | |
| 810 **NOTE**: Tree is UNROOTED although outgroup taxon 'LngfishAu' is drawn at root | |
| 811 | |
| 812 :: | |
| 813 | |
| 814 +--------------LngfishAu | |
| 815 | | |
| 816 | +--------------LngfishSA | |
| 817 +--------| | |
| 818 | +--------------LngfishAf | |
| 819 | | |
| 820 | +-------------------Frog | |
| 821 +------| | |
| 822 | +-----------------Turtle | |
| 823 | +-----| | |
| 824 | | | +-----------------------Sphenodon | |
| 825 | | | +--| | |
| 826 | | | | +--------------------------Lizard | |
| 827 | | +---| | |
| 828 | | | +---------------------Crocodile | |
| 829 | | +------| | |
| 830 | | +------------------Bird | |
| 831 +---------| | |
| 832 | +----------------Human | |
| 833 | +--| | |
| 834 | | | +--------Seal | |
| 835 | | +--| | |
| 836 | | | +-------Cow | |
| 837 | | +---| | |
| 838 | | +---------Whale | |
| 839 | +----| | |
| 840 | | | +------Mouse | |
| 841 | | +---------| | |
| 842 | | +--------Rat | |
| 843 +----------| | |
| 844 | +----------------Platypus | |
| 845 +---| | |
| 846 +-------------Opossum | |
| 847 | |
| 848 | |
| 849 This makes sense as the mammals (**Human** to **Opossum**) form a clade, whereas the reptiles **Turtle** to **Crocodile**) and **Bird** form a separate sister clade. Here the tree is drawn at the *outgroup* Lungfish which is more accient than other species in this example. However, please note that IQ-TREE always produces an **unrooted tree** as it knows nothing about this biological background; IQ-TREE simply draws the tree this way as **LngfishAu** is the first sequence occuring in the alignment. | |
| 850 | |
| 851 Choosing the right substitution model | |
| 852 ------------------------------------- | |
| 853 | |
| 854 IQ-TREE will choose the best model for you automatically if specify any of the TEST models, but valid custom models can also be specified that conform to those found in Models_ page. | |
| 855 | |
| 856 .. _Models: http://www.iqtree.org/doc/Substitution-Models#binary-and-morphological-models | |
| 857 | |
| 858 | |
| 859 ------------------------------ | |
| 860 | |
| 861 Advanced Parameter Selection | |
| 862 ============================ | |
| 863 | |
| 864 Using codon models | |
| 865 ------------------ | |
| 866 | |
| 867 IQ-TREE supports a number of codon models. You need to input a protein-coding DNA alignment and specify codon data by option **-st CODON** (Otherwise, IQ-TREE applies DNA model because it detects that your alignment has DNA sequences): | |
| 868 | |
| 869 :: | |
| 870 | |
| 871 iqtree -s coding_gene.phy -st CODON | |
| 872 | |
| 873 If your alignment length is not divisible by 3, IQ-TREE will stop with an error message. IQ-TREE will group sites 1,2,3 into codon site 1; sites 4,5,6 to codon site 2; etc. Moreover, any codon, which has at least one gap/unknown/ambiguous nucleotide, will be treated as unknown codon character. | |
| 874 | |
| 875 Note that the above command assumes the standard genetic code. If your sequences follow 'The Invertebrate Mitochondrial Code', then run: | |
| 876 | |
| 877 :: | |
| 878 | |
| 879 iqtree -s coding_gene.phy -st CODON5 | |
| 880 | |
| 881 Note that ModelFinder works for codon alignments. IQ-TREE version >= 1.5.4 will automatically invokes ModelFinder to find the best-fit codon model. | |
| 882 | |
| 883 | |
| 884 Assessing branch supports with ultrafast bootstrap approximation | |
| 885 ---------------------------------------------------------------- | |
| 886 | |
| 887 To overcome the computational burden required by the nonparametric bootstrap, IQ-TREE introduces an ultrafast bootstrap approximation (UFBoot) ([Minh et al., 2013]) that is orders of magnitude faster than the standard procedure and provides relatively unbiased branch support values. Citation for UFBoot: | |
| 888 | |
| 889 B.Q. Minh, M.A.T. Nguyen, and A. von Haeseler (2013) Ultrafast approximation for phylogenetic bootstrap. _Mol. Biol. Evol., 30:1188-1195. | |
| 890 | |
| 891 :: | |
| 892 | |
| 893 iqtree -s example.phy -m TIM2+I+G -bb 1000 | |
| 894 | |
| 895 **-bb** specifies the number of bootstrap replicates where 1000 is the minimum number recommended. The section **MAXIMUM LIKELIHOOD TREE** in **example.phy.iqtree** shows a textual representation of the maximum likelihood tree with branch support values in percentage. The NEWICK format of the tree is printed to the file **example.phy.treefile**. In addition, IQ-TREE writes the following files: | |
| 896 | |
| 897 * **example.phy.contree**: the consensus tree with assigned branch supports where branch lengths are optimized on the original alignment. | |
| 898 | |
| 899 * **example.phy.splits**: support values in percentage for all splits (bipartitions), computed as the occurence frequencies in the bootstrap trees. This file is in "star-dot" format. | |
| 900 | |
| 901 * **example.phy.splits.nex**: has the same information as **example.phy.splits** but in NEXUS format, which can be viewed with the program SplitsTree_ to explore the conflicting signals in the data. So it is more informative than consensus tree, e.g. you can see how highly supported the second best conflicting split is, which had no chance to enter the consensus tree. | |
| 902 | |
| 903 .. _SplitsTree: http://www.splitstree.org | |
| 904 | |
| 905 | |
| 906 Reducing impact of severe model violations with UFBoot | |
| 907 ------------------------------------------------------ | |
| 908 | |
| 909 Starting with IQ-TREE version 1.6 we provide a new option **-bnni** to reduce the risk of overestimating branch supports with UFBoot due to severe model violations. With this option UFBoot will further optimize each bootstrap tree using a hill-climbing nearest neighbor interchange (NNI) search based directly on the corresponding bootstrap alignment. | |
| 910 | |
| 911 Thus, if severe model violations are present in the data set at hand, users are advised to append **-bnni** to the regular UFBoot command: | |
| 912 | |
| 913 iqtree -s example.phy -m TIM2+I+G -bb 1000 -bnni | |
| 914 | |
| 915 For more details see: | |
| 916 | |
| 917 D.T. Hoang, O. Chernomor, A. von Haeseler, B.Q. Minh, L.S. Vinh (2017) UFBoot2: Improving the ultrafast bootstrap approximation. | |
| 918 | |
| 919 | |
| 920 Assessing branch supports with standard nonparametric bootstrap | |
| 921 ---------------------------------------------------------------- | |
| 922 | |
| 923 The standard nonparametric bootstrap is invoked by the **-b** option: | |
| 924 | |
| 925 iqtree -s example.phy -m TIM2+I+G -b 100 | |
| 926 | |
| 927 **-b** specifies the number of bootstrap replicates where 100 | |
| 928 is the minimum recommended number. The output files are similar to those produced by the UFBoot procedure. | |
| 929 | |
| 930 | |
| 931 | |
| 932 Assessing branch supports with single branch tests | |
| 933 -------------------------------------------------- | |
| 934 | |
| 935 IQ-TREE provides an implementation of the SH-like approximate likelihood ratio test ([Guindon et al., 2010]). To perform this test, run: | |
| 936 | |
| 937 :: | |
| 938 | |
| 939 iqtree -s example.phy -m TIM2+I+G -alrt 1000 | |
| 940 | |
| 941 **-alrt** specifies the number of bootstrap replicates for SH-aLRT where 1000 is the minimum number recommended. | |
| 942 | |
| 943 IQ-TREE also supports other tests such as the aBayes test (Anisimova et al., 2011) and the local bootstrap test (Adachi and Hasegawa, 1996). | |
| 944 | |
| 945 You can also perform both SH-aLRT and the ultrafast bootstrap within one single run: | |
| 946 | |
| 947 :: | |
| 948 | |
| 949 iqtree -s example.phy -m TIM2+I+G -alrt 1000 -bb 1000 | |
| 950 | |
| 951 The branches of the resulting **.treefile** will be assigned with both SH-aLRT and UFBoot support values, which are readable by any tree viewer program like FigTree, Dendroscope or ETE. You can also look at the textual tree figure in **.iqtree** file: | |
| 952 | |
| 953 **NOTE**: Tree is UNROOTED although outgroup taxon 'LngfishAu' is drawn at root | |
| 954 Numbers in parentheses are SH-aLRT support (%) / ultrafast bootstrap support (%) | |
| 955 | |
| 956 :: | |
| 957 | |
| 958 +-------------LngfishAu | |
| 959 | | |
| 960 | +--------------LngfishSA | |
| 961 +-------| (100/100) | |
| 962 | +------------LngfishAf | |
| 963 | | |
| 964 | +--------------------Frog | |
| 965 +------| (99.8/100) | |
| 966 | +-----------------Turtle | |
| 967 | +--| (85/72) | |
| 968 | | | +------------------------Crocodile | |
| 969 | | +----| (96.5/97) | |
| 970 | | +------------------Bird | |
| 971 | +--| (39/51) | |
| 972 | | +---------------------------Sphenodon | |
| 973 | +-----| (98.2/99) | |
| 974 | | +-------------------------------Lizard | |
| 975 +---------| (100/100) | |
| 976 | +--------------Human | |
| 977 | +--| (92.3/93) | |
| 978 | | | +------Seal | |
| 979 | | +--| (68.3/75) | |
| 980 | | | +-----Cow | |
| 981 | | +--| (99.7/100) | |
| 982 | | +-------Whale | |
| 983 | +----| (99.1/100) | |
| 984 | | | +---Mouse | |
| 985 | | +---------| (100/100) | |
| 986 | | +------Rat | |
| 987 +-----------| (100/100) | |
| 988 | +--------------Platypus | |
| 989 +--| (93/98) | |
| 990 +-----------Opossum | |
| 991 | |
| 992 | |
| 993 From this figure, the branching patterns within reptiles are poorly supported (e.g. **Sphenodon** with SH-aLRT: 39%, UFBoot: 51% and **Turtle** with SH-aLRT: 85%, UFBoot: 72%) as well as the phylogenetic position of **Seal** within mammals (SH-aLRT: 68.3%, UFBoot: 75%). Other branches appear to be well supported. | |
| 994 ]]></help> | |
| 995 <citations> | |
| 996 <citation type="doi">10.1093/molbev/msu300</citation> | |
| 997 <citation type="doi">10.1093/molbev/mst024</citation> | |
| 998 </citations> | |
| 999 </tool> |