RAxml-8.2.X使用说明

This is RAxML version 8.2.12 released by Alexandros Stamatakis on May 2018.

简单示例：raxmlHPC-PTHREADS -s 比对文件 -n 输出文件名 -m 替代模型（核苷酸或氨基酸替代模型+变异速率模型） -T 线程数 -p 12345 -f a -N 快速自展值 -x 12345

-s sequenceFileName -n outputFileName -m substitutionModel

[-a weightFileName] [-A secondaryStructureSubstModel]

[-b bootstrapRandomNumberSeed] [-B wcCriterionThreshold]

[-c numberOfCategories] [-C] [-d] [-D]

[-e likelihoodEpsilon] [-E excludeFileName]

[-f a|A|b|B|c|C|d|D|e|E|F|g|G|h|H|i|I|j|J|k|m|n|N|o|p|P|q|r|R|s|S|t|T|u|v|V|w|W|x|y] [-F]

[-g groupingFileName] [-G placementThreshold] [-h] [-H]

[-i initialRearrangementSetting] [-I autoFC|autoMR|autoMRE|autoMRE_IGN]

[-L MR|MRE|T_<PERCENT>] [-M]

[-o outGroupName1[,outGroupName2[,...]]][-O]

[-p parsimonyRandomSeed] [-P proteinModel]

[-q multipleModelFileName] [-r binaryConstraintTree]

[-R binaryModelParamFile] [-S secondaryStructureFile] [-t userStartingTree]

[-T numberOfThreads] [-u] [-U] [-v] [-V] [-w outputDirectory] [-W slidingWindowSize]

[-x rapidBootstrapRandomNumberSeed] [-X] [-y] [-Y quartetGroupingFileName|ancestralSequenceCandidatesFileName]

[--mesquite][--silent][--no-seq-check][--no-bfgs]

[--asc-corr=stamatakis|felsenstein|lewis]

[--flag-check][--auto-prot=ml|bic|aic|aicc]

[--epa-keep-placements=number][--epa-accumulated-threshold=threshold]

[--epa-prob-threshold=threshold]

[--JC69][--K80][--HKY85]

[--bootstop-perms=number]

[--quartets-without-replacement]

[---without-replacement]

[--print-identical-sequences]

参数说明：

-a Specify a column weight file name to assign individual weights to each column of 设定每个位点的权重

the alignment. Those weights must be integers separated by any type and number

of whitespaces whithin a separate file, see file "example_weights" for an example.

-A Specify one of the secondary structure substitution models implemented in RAxML. 指定一个二级结构替换模型

The same nomenclature as in the PHASE manual is used, available models:

S6A, S6B, S6C, S6D, S6E, S7A, S7B, S7C, S7D, S7E, S7F, S16, S16A, S16B

DEFAULT: 16-state GTR model (S16)

-b Specify an integer number (random seed) and turn on bootstrapping 设定bootstrap起始随机数，默认是关闭的。-b 12345

DEFAULT: OFF

-B specify a floating point number between 0.0 and 1.0 that will be used as cutoff threshold 用于切断基于MR的自举停止标准，对大多数一致树的收敛有积极作用

for the MR-based bootstopping criteria. The recommended setting is 0.03. 默认值是0.03。可调节范围是0.0~1.0。

DEFAULT: 0.03 (recommended empirically determined setting) 一般情况下，我们用默认值就行，默认值0.03是基于经验得出来的，在实际中不需要输入-B。

-c Specify number of distinct rate catgories for RAxML when model of rate heterogeneity 当使用CAT变异速率模型时去指定distinct rate categories，默认值是25

is set to CAT 注意，如果你使用的是Gamma变异速率模型，不管你将-c 设置为多少值都是没有作用的，因为Gamma只有4个变异速率。

Individual per-site rates are categorized into numberOfCategories rate

categories to accelerate computations.

DEFAULT: 25

-C Enable verbose output for the "-L" and "-f i" options. This will produce more, as well as more verbose output files 激活“-L”和“-f i”选项的冗长输出，默认是关闭。

DEFAULT: OFF

-d start ML optimization from random starting tree 完全随机的搜索进化树，而不是从maximum parsimony tree 开始，在100 至200 个分类单元间，该选项可能会生成拓扑结构完全不同的局部最大似然树。默认值是关闭。

DEFAULT: OFF

-D ML search convergence criterion. This will break off ML searches if the relative ML法搜寻收敛的标准。在非常非常大的数据集或许可以采用这个，

Robinson-Foulds distance between the trees obtained from two consecutive lazy SPR cycles 但是可能找到的是不好的树，所以这个选项默认在关闭的。

is smaller or equal to 1%. Usage recommended for very large datasets in terms of taxa.

On trees with more than 500 taxa this will yield execution time improvements of approximately 50%

While yielding only slightly worse trees.

DEFAULT: OFF

-e set model optimization precision in log likelihood units for final

optimization of tree topology

DEFAULT: 0.1 for models not using proportion of invariant sites estimate

0.001 for models using proportion of invariant sites estimate

-E specify an exclude file name, that contains a specification of alignment positions you wish to exclude.

Format is similar to Nexus, the file shall contain entries like "100-200 300-400", to exclude a

single column write, e.g., "100-100", if you use a mixed model, an appropriately adapted model file

will be written.

-f select algorithm:

"-f a": rapid Bootstrap analysis and search for best-scoring ML tree in one program run

"-f A": compute marginal ancestral states on a ROOTED reference tree provided with "-t"

"-f b": draw bipartition information on a tree provided with "-t" based on multiple trees

(e.g., from a bootstrap) in a file specified by "-z"

"-f B": optimize br-len scaler and other model parameters (GTR, alpha, etc.) on a tree provided with "-t".

The tree needs to contain branch lengths. The branch lengths will not be optimized, just scaled by a single common value.

"-f c": check if the alignment can be properly read by RAxML

"-f C": ancestral sequence test for Jiajie, users will also need to provide a list of taxon names via -Y separated by whitespaces

"-f d": new rapid hill-climbing

DEFAULT: ON

"-f D": rapid hill-climbing with RELL bootstraps

"-f e": optimize model+branch lengths for given input tree under GAMMA/GAMMAI only

"-f E": execute very fast experimental tree search, at present only for testing

"-f F": execute fast experimental tree search, at present only for testing

"-f g": compute per site log Likelihoods for one or more trees passed via

"-z" and write them to a file that can be read by CONSEL

The model parameters will be estimated on the first tree only!

"-f G": compute per site log Likelihoods for one or more trees passed via

"-z" and write them to a file that can be read by CONSEL.

The model parameters will be re-estimated for each tree

"-f h": compute log likelihood test (SH-test) between best tree passed via "-t"

and a bunch of other trees passed via "-z"

The model parameters will be estimated on the first tree only!

"-f H": compute log likelihood test (SH-test) between best tree passed via "-t"

and a bunch of other trees passed via "-z"

The model parameters will be re-estimated for each tree

"-f i": calculate IC and TC scores (Salichos and Rokas 2013) on a tree provided with "-t" based on multiple trees

(e.g., from a bootstrap) in a file specified by "-z"

"-f I": a simple tree rooting algorithm for unrooted trees.

It roots the tree by rooting it at the branch that best balances the subtree lengths

(sum over branches in the subtrees) of the left and right subtree.

A branch with an optimal balance does not always exist!

You need to specify the tree you want to root via "-t".

"-f j": generate a bunch of bootstrapped alignment files from an original alignemnt file.

You need to specify a seed with "-b" and the number of replicates with "-#"

"-f J": Compute SH-like support values on a given tree passed via "-t".

"-f k": Fix long branch lengths in partitioned data sets with missing data using the

branch length stealing algorithm.

This option only works in conjunction with "-t", "-M", and "-q".

It will print out a tree with shorter branch lengths, but having the same likelihood score.

"-f m": compare bipartitions between two bunches of trees passed via "-t" and "-z"

respectively. This will return the Pearson correlation between all bipartitions found

in the two tree files. A file called RAxML_bipartitionFrequencies.outpuFileName

will be printed that contains the pair-wise bipartition frequencies of the two sets

"-f n": compute the log likelihood score of all trees contained in a tree file provided by

"-z" under GAMMA or GAMMA+P-Invar

The model parameters will be estimated on the first tree only!

"-f N": compute the log likelihood score of all trees contained in a tree file provided by

"-z" under GAMMA or GAMMA+P-Invar

The model parameters will be re-estimated for each tree

"-f o": old and slower rapid hill-climbing without heuristic cutoff

"-f p": perform pure stepwise MP addition of new sequences to an incomplete starting tree and exit

"-f P": perform a phylogenetic placement of sub trees specified in a file passed via "-z" into a given reference tree

in which these subtrees are contained that is passed via "-t" using the evolutionary placement algorithm.

"-f q": fast quartet calculator

"-f r": compute pairwise Robinson-Foulds (RF) distances between all pairs of trees in a tree file passed via "-z"

if the trees have node labales represented as integer support values the program will also compute two flavors of

the weighted Robinson-Foulds (WRF) distance

"-f R": compute all pairwise Robinson-Foulds (RF) distances between a large reference tree passed via "-t"

and many smaller trees (that must have a subset of the taxa of the large tree) passed via "-z".

This option is intended for checking the plausibility of very large phylogenies that can not be inspected

visually any more.

"-f s": split up a multi-gene partitioned alignment into the respective subalignments

"-f S": compute site-specific placement bias using a leave one out test inspired by the evolutionary placement algorithm

"-f t": do randomized tree searches on one fixed starting tree

"-f T": do final thorough optimization of ML tree from rapid bootstrap search in stand-alone mode

"-f u": execute morphological weight calibration using maximum likelihood, this will return a weight vector.

you need to provide a morphological alignment and a reference tree via "-t"

"-f v": classify a bunch of environmental sequences into a reference tree using thorough read insertions

you will need to start RAxML with a non-comprehensive reference tree and an alignment containing all sequences (reference + query)

"-f V": classify a bunch of environmental sequences into a reference tree using thorough read insertions

you will need to start RAxML with a non-comprehensive reference tree and an alignment containing all sequences (reference + query)

WARNING: this is a test implementation for more efficient handling of multi-gene/whole-genome datasets!

"-f w": compute ELW test on a bunch of trees passed via "-z"

The model parameters will be estimated on the first tree only!

"-f W": compute ELW test on a bunch of trees passed via "-z"

The model parameters will be re-estimated for each tree

"-f x": compute pair-wise ML distances, ML model parameters will be estimated on an MP

starting tree or a user-defined tree passed via "-t", only allowed for GAMMA-based

models of rate heterogeneity

"-f y": classify a bunch of environmental sequences into a reference tree using parsimony

you will need to start RAxML with a non-comprehensive reference tree and an alignment containing all sequences (reference + query)

DEFAULT for "-f": new rapid hill climbing

-F enable ML tree searches under CAT model for very large trees without switching to

GAMMA in the end (saves memory).

This option can also be used with the GAMMA models in order to avoid the thorough optimization

of the best-scoring ML tree in the end.

DEFAULT: OFF

-g specify the file name of a multifurcating constraint tree

this tree does not need to be comprehensive, i.e. must not contain all taxa

-G enable the ML-based evolutionary placement algorithm heuristics

by specifiyng a threshold value (fraction of insertion branches to be evaluated

using slow insertions under ML).

-h Display this help message.

-H Disable pattern compression.

DEFAULT: ON

-i Initial rearrangement setting for the subsequent application of topological

changes phase

-I a posteriori bootstopping analysis. Use:

"-I autoFC" for the frequency-based criterion

"-I autoMR" for the majority-rule consensus tree criterion

"-I autoMRE" for the extended majority-rule consensus tree criterion

"-I autoMRE_IGN" for metrics similar to MRE, but include bipartitions under the threshold whether they are compatible

or not. This emulates MRE but is faster to compute.

You also need to pass a tree file containg several bootstrap replicates via "-z"

-j Specifies that intermediate tree files shall be written to file during the standard ML and BS tree searches.

DEFAULT: OFF

-J Compute majority rule consensus tree with "-J MR" or extended majority rule consensus tree with "-J MRE"

or strict consensus tree with "-J STRICT". For a custom consensus threshold >= 50%, specify T_<NUM>, where 100 >= NUM >= 50.

Options "-J STRICT_DROP" and "-J MR_DROP" will execute an algorithm that identifies dropsets which contain

rogue taxa as proposed by Pattengale et al. in the paper "Uncovering hidden phylogenetic consensus".

You will also need to provide a tree file containing several UNROOTED trees via "-z"

-k Specifies that bootstrapped trees should be printed with branch lengths.

The bootstraps will run a bit longer, because model parameters will be optimized

at the end of each run under GAMMA or GAMMA+P-Invar respectively.

DEFAULT: OFF

-K Specify one of the multi-state substitution models (max 32 states) implemented in RAxML.

Available models are: ORDERED, MK, GTR

DEFAULT: GTR model

-L Compute consensus trees labelled by IC supports and the overall TC value as proposed in Salichos and Rokas 2013.

Compute a majority rule consensus tree with "-L MR" or an extended majority rule consensus tree with "-L MRE".

For a custom consensus threshold >= 50%, specify "-L T_<NUM>", where 100 >= NUM >= 50.

You will of course also need to provide a tree file containing several UNROOTED trees via "-z"!

-m Model of Binary (Morphological), Nucleotide, Multi-State, or Amino Acid Substitution:

BINARY:

"-m BINCAT[X]" : Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under BINGAMMA, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m BINCATI[X]" : Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under BINGAMMAI, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_BINCAT[X]" : Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under BINGAMMA, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

The ASC prefix willl correct the likelihood for ascertainment bias.

"-m BINGAMMA[X]" : GAMMA model of rate heterogeneity (alpha parameter will be estimated).

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_BINGAMMA[X]" : GAMMA model of rate heterogeneity (alpha parameter will be estimated).

The ASC prefix willl correct the likelihood for ascertainment bias.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m BINGAMMAI[X]" : Same as BINGAMMA, but with estimate of proportion of invariable sites.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

NUCLEOTIDES:

"-m GTRCAT[X]" : GTR + Optimization of substitution rates + Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

under GTRGAMMA, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m GTRCATI[X]" : GTR + Optimization of substitution rates + Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

under GTRGAMMAI, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_GTRCAT[X]" : GTR + Optimization of substitution rates + Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

under GTRGAMMA, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

The ASC prefix willl correct the likelihood for ascertainment bias.

"-m GTRGAMMA[X]" : GTR + Optimization of substitution rates + GAMMA model of rate

heterogeneity (alpha parameter will be estimated).

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_GTRGAMMA[X]" : GTR + Optimization of substitution rates + GAMMA model of rate

heterogeneity (alpha parameter will be estimated).

The ASC prefix willl correct the likelihood for ascertainment bias.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m GTRGAMMAI[X]" : Same as GTRGAMMA, but with estimate of proportion of invariable sites.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

MULTI-STATE:

"-m MULTICAT[X]" : Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under MULTIGAMMA, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m MULTICATI[X]" : Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under MULTIGAMMAI, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_MULTICAT[X]" : Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under MULTIGAMMA, depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

The ASC prefix willl correct the likelihood for ascertainment bias.

"-m MULTIGAMMA[X]" : GAMMA model of rate heterogeneity (alpha parameter will be estimated).

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_MULTIGAMMA[X]" : GAMMA model of rate heterogeneity (alpha parameter will be estimated).

The ASC prefix willl correct the likelihood for ascertainment bias.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m MULTIGAMMAI[X]" : Same as MULTIGAMMA, but with estimate of proportion of invariable sites.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

You can use up to 32 distinct character states to encode multi-state regions, they must be used in the following order:

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, A, B, C, D, E, F, G, H, I, J, K, L, M, N, O, P, Q, R, S, T, U, V

i.e., if you have 6 distinct character states you would use 0, 1, 2, 3, 4, 5 to encode these.

The substitution model for the multi-state regions can be selected via the "-K" option

AMINO ACIDS:

"-m PROTCATmatrixName[F|X]" : specified AA matrix + Optimization of substitution rates + Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under PROTGAMMAmatrixName[F|X], depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m PROTCATImatrixName[F|X]" : specified AA matrix + Optimization of substitution rates + Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under PROTGAMMAImatrixName[F|X], depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_PROTCATmatrixName[F|X]" : specified AA matrix + Optimization of substitution rates + Optimization of site-specific

evolutionary rates which are categorized into numberOfCategories distinct

rate categories for greater computational efficiency. Final tree might be evaluated

automatically under PROTGAMMAmatrixName[F|X], depending on the tree search option.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

The ASC prefix willl correct the likelihood for ascertainment bias.

"-m PROTGAMMAmatrixName[F|X]" : specified AA matrix + Optimization of substitution rates + GAMMA model of rate

heterogeneity (alpha parameter will be estimated).

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m ASC_PROTGAMMAmatrixName[F|X]" : specified AA matrix + Optimization of substitution rates + GAMMA model of rate

heterogeneity (alpha parameter will be estimated).

The ASC prefix willl correct the likelihood for ascertainment bias.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

"-m PROTGAMMAImatrixName[F|X]" : Same as PROTGAMMAmatrixName[F|X], but with estimate of proportion of invariable sites.

With the optional "X" appendix you can specify a ML estimate of base frequencies.

Available AA substitution models:

DAYHOFF, DCMUT, JTT, MTREV, WAG, RTREV, CPREV, VT,

BLOSUM62, MTMAM, LG, MTART, MTZOA, PMB, HIVB, HIVW,

JTTDCMUT, FLU, STMTREV, DUMMY, DUMMY2, AUTO, LG4M, LG4X,

PROT_FILE, GTR_UNLINKED, GTR

With the optional "F" appendix you can specify if you want to use empirical base frequencies.

AUTOF and AUTOX are not supported any more, if you specify AUTO it will test prot subst. models with and without empirical

base frequencies now!

Please note that for partitioned models you can in addition specify the per-gene AA model in

the partition file (see manual for details). Also note that if you estimate AA GTR parameters on a partitioned

dataset, they will be linked (estimated jointly) across all partitions to avoid over-parametrization

-M Switch on estimation of individual per-partition branch lengths. Only has effect when used in combination with "-q"

Branch lengths for individual partitions will be printed to separate files

A weighted average of the branch lengths is computed by using the respective partition lengths

DEFAULT: OFF

-n Specifies the name of the output file.

-o Specify the name of a single outgroup or a comma-separated list of outgroups, eg "-o Rat"

or "-o Rat,Mouse", in case that multiple outgroups are not monophyletic the first name

in the list will be selected as outgroup, don't leave spaces between taxon names!

-O Disable check for completely undetermined sequence in alignment.

The program will not exit with an error message when "-O" is specified.

DEFAULT: check enabled

-p Specify a random number seed for the parsimony inferences. This allows you to reproduce your results

and will help me debug the program.

-P Specify the file name of a user-defined AA (Protein) substitution model. This file must contain

420 entries, the first 400 being the AA substitution rates (this must be a symmetric matrix) and the

last 20 are the empirical base frequencies

-q Specify the file name which contains the assignment of models to alignment

partitions for multiple models of substitution. For the syntax of this file

please consult the manual.

-r Specify the file name of a binary constraint tree.

this tree does not need to be comprehensive, i.e. must not contain all taxa

-R Specify the file name of a binary model parameter file that has previously been generated

with RAxML using the -f e tree evaluation option. The file name should be:

RAxML_binaryModelParameters.runID

-s Specify the name of the alignment data file in PHYLIP format

-S Specify the name of a secondary structure file. The file can contain "." for

alignment columns that do not form part of a stem and characters "()<>[]{}" to define

stem regions and pseudoknots

-t Specify a user starting tree file name in Newick format

-T PTHREADS VERSION ONLY! Specify the number of threads you want to run.

Make sure to set "-T" to at most the number of CPUs you have on your machine,

otherwise, there will be a huge performance decrease!

-u use the median for the discrete approximation of the GAMMA model of rate heterogeneity

DEFAULT: OFF

-U Try to save memory by using SEV-based implementation for gap columns on large gappy alignments

The technique is described here: http://www.biomedcentral.com/1471-2105/12/470

This will only work for DNA and/or PROTEIN data and only with the SSE3 or AVX-vextorized version of the code.

-v Display version information

-V Disable rate heterogeneity among sites model and use one without rate heterogeneity instead.

Only works if you specify the CAT model of rate heterogeneity.

DEFAULT: use rate heterogeneity

-w FULL (!) path to the directory into which RAxML shall write its output files

DEFAULT: current directory

-W Sliding window size for leave-one-out site-specific placement bias algorithm

only effective when used in combination with "-f S"

DEFAULT: 100 sites

-x Specify an integer number (random seed) and turn on rapid bootstrapping

CAUTION: unlike in version 7.0.4 RAxML will conduct rapid BS replicates under

the model of rate heterogeneity you specified via "-m" and not by default under CAT

-X Same as the "-y" option below, however the parsimony search is more superficial.

RAxML will only do a randomized stepwise addition order parsimony tree reconstruction

without performing any additional SPRs.

This may be helpful for very broad whole-genome datasets, since this can generate topologically

more different starting trees.

DEFAULT: OFF

-y If you want to only compute a parsimony starting tree with RAxML specify "-y",

the program will exit after computation of the starting tree

DEFAULT: OFF

-Y Pass a quartet grouping file name defining four groups from which to draw quartets

The file input format must contain 4 groups in the following form:

(Chicken, Human, Loach), (Cow, Carp), (Mouse, Rat, Seal), (Whale, Frog);

Only works in combination with -f q !

-z Specify the file name of a file containing multiple trees e.g. from a bootstrap

that shall be used to draw bipartition values onto a tree provided with "-t",

It can also be used to compute per site log likelihoods in combination with "-f g"

and to read a bunch of trees for a couple of other options ("-f h", "-f m", "-f n").

-#|-N Specify the number of alternative runs on distinct starting trees

In combination with the "-b" option, this will invoke a multiple boostrap analysis

Note that "-N" has been added as an alternative since "-#" sometimes caused problems

with certain MPI job submission systems, since "-#" is often used to start comments.

If you want to use the bootstopping criteria specify "-# autoMR" or "-# autoMRE" or "-# autoMRE_IGN"

for the majority-rule tree based criteria (see -I option) or "-# autoFC" for the frequency-based criterion.

Bootstopping will only work in combination with "-x" or "-b"

DEFAULT: 1 single analysis

--mesquite Print output files that can be parsed by Mesquite.

DEFAULT: Off

--silent Disables printout of warnings related to identical sequences and entirely undetermined sites in the alignment

DEFAULT: Off

--no-seq-check Disables checking the input MSA for identical sequences and entirely undetermined sites.

Enabling this option may save time, in particular for large phylogenomic alignments.

Before using this, make sure to check the alignment using the "-f c" option!

DEFAULT: Off

--no-bfgs Disables automatic usage of BFGS method to optimize GTR rates on unpartitioned DNA datasets

DEFAULT: BFGS on

--asc-corr Allows to specify the type of ascertainment bias correction you wish to use. There are 3

types available:

--asc-corr=lewis: the standard correction by Paul Lewis

--asc-corr=felsenstein: a correction introduced by Joe Felsenstein that allows to explicitely specify

the number of invariable sites (if known) one wants to correct for.

--asc-corr=stamatakis: a correction introduced by myself that allows to explicitely specify

the number of invariable sites for each character (if known) one wants to correct for.

--flag-check When using this option, RAxML will only check if all command line flags specifed are available and then exit

with a message listing all invalid command line flags or with a message stating that all flags are valid.

--auto-prot=ml|bic|aic|aicc When using automatic protein model selection you can chose the criterion for selecting these models.

RAxML will test all available prot subst. models except for LG4M, LG4X and GTR-based models, with and without empirical base frequencies.

You can chose between ML score based selection and the BIC, AIC, and AICc criteria.

DEFAULT: ml

--epa-keep-placements=number specify the number of potential placements you want to keep for each read in the EPA algorithm.

Note that, the actual values printed will also depend on the settings for --epa-prob-threshold=threshold !

DEFAULT: 7

--epa-prob-threshold=threshold specify a percent threshold for including potential placements of a read depending on the

maximum placement weight for this read. If you set this value to 0.01 placements that have a placement weight of 1 per cent of

the maximum placement will still be printed to file if the setting of --epa-keep-placements allows for it

DEFAULT: 0.01

--epa-accumulated-threshold=threshold specify an accumulated likelihood weight threshold for which different placements of read are printed

to file. Placements for a read will be printed until the sum of their placement weights has reached the threshold value.

Note that, this option can neither be used in combination with --epa-prob-threshold nor with --epa-keep-placements!

--JC69 specify that all DNA partitions will evolve under the Jukes-Cantor model, this overrides all other model specifications for DNA partitions.

DEFAULT: Off

--K80 specify that all DNA partitions will evolve under the K80 model, this overrides all other model specifications for DNA partitions.

DEFAULT: Off

--HKY85 specify that all DNA partitions will evolve under the HKY85 model, this overrides all other model specifications for DNA partitions.

DEFAULT: Off

--bootstop-perms=number specify the number of permutations to be conducted for the bootstopping/bootstrap convergence test.

The allowed minimum number is 100!

DEFAULT: 100

--quartets-without-replacement specify that quartets are randomly subsampled, but without replacement.

DEFAULT: random sampling with replacements

--print-identical-sequences specify that RAxML shall automatically generate a .reduced alignment with all

undetermined columns removed, but without removing exactly identical sequences

DEFAULT: identical sequences will also be removed in the .reduced file

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 202,056评论 5赞 474
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 84,842评论 2赞 378
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 148,938评论 0赞 335
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 54,296评论 1赞 272
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 63,292评论 5赞 363
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 48,413评论 1赞 281
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 37,824评论 3赞 393
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 36,493评论 0赞 256
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 40,686评论 1赞 295
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 35,502评论 2赞 318
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 37,553评论 1赞 329
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 33,281评论 4赞 318
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 38,820评论 3赞 305
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 29,873评论 0赞 19
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 31,109评论 1赞 258
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 42,699评论 2赞 348
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 42,257评论 2赞 341

RAxml-8.2.X使用说明

推荐阅读更多精彩内容