|
|
ednaml |
Please help by correcting and extending the Wiki pages.
% ednaml -categories -catnum 2 -catvals "1.0 3.2" -catprob "0.4 0.6" -lambda 1.5 -deflambda -printdata -progress Phylogenies from nucleic acid maximum likelihood Input (aligned) nucleotide sequence set: dnaml.dat Search for best tree [Y]: Global rearrangements [N]: Randomize input order of species [N]: Use empirical base frequencies [Y]: Outgroup root [N]: Phylip dnaml program output file [ednaml.outfile]: Create a tree file [Y]: Phylip tree output file [ednaml.treefile]: Adding species: Alpha Beta Gamma Delta Epsilon Output written to output file Tree also written onto file |
Go to the input files for this example
Go to the output files for this example
Phylogenies from nucleic acid maximum likelihood
Version: EMBOSS:6.6.0.0
Standard (Mandatory) qualifiers (* if not always prompted):
[-sequence] seqset File containing a sequence alignment
-[no]besttree toggle [Y] Search for best tree
* -lengths boolean [N] Use lengths from user trees
* -global boolean [N] Global rearrangements
* -random toggle [N] Randomize input order of species
* -randseed integer [3] Random number seed (must be odd) (Any
integer value)
* -randtimes integer [3] How many times to randomise (Any integer
value)
-[no]basefrequency toggle [Y] Use empirical base frequencies
* -freqa float [0.25] Frequency for A (Number from 0.000 to
1.000)
* -freqc float [0.25] Frequency for C (Number from 0.000 to
1.000)
* -freqg float [0.25] Frequency for G (Number from 0.000 to
1.000)
* -freqt float [0.25] Frequency for T/U (Number from 0.000
to 1.000)
-categories toggle [N] More then one category of substitution
rates
* -catnum integer [2] Number of categories of substitution
rates (Integer from 2 to 9)
* -catvals string Space separated category values (Any string)
* -catprob string Space separated probability values (Any
string)
* -deflambda toggle [N] Use default mean block length of sites
* -lambda float [1.0] Mean block length of sites having the
same rate (greater than 1) (Number 1.000 or
more)
-og toggle [N] Outgroup root
* -outgnum integer [1] Number of the outgroup (Integer 1 or
more)
[-outfile] outfile [ednaml.outfile] Phylip dnaml program output
file
[-[no]trout] toggle [Y] Create a tree file
* -treefile outfile [ednaml.treefile] Phylip tree output file
Additional (Optional) qualifiers:
-ttratio float [2.0] Transition/transversion ratio (Any
numeric value)
-printdata boolean [N] Print out the data at start of run
-progress boolean [N] Print indications of progress of run
-[no]drawtree boolean [Y] Draw tree
Advanced (Unprompted) qualifiers: (none)
Associated qualifiers:
"-sequence" associated qualifiers
-sbegin1 integer Start of each sequence to be used
-send1 integer End of each sequence to be used
-sreverse1 boolean Reverse (if DNA)
-sask1 boolean Ask for begin/end/reverse
-snucleotide1 boolean Sequence is nucleotide
-sprotein1 boolean Sequence is protein
-slower1 boolean Make lower case
-supper1 boolean Make upper case
-scircular1 boolean Sequence is circular
-squick1 boolean Read id and sequence only
-sformat1 string Input sequence format
-iquery1 string Input query fields or ID list
-ioffset1 integer Input start position offset
-sdbname1 string Database name
-sid1 string Entryname
-ufo1 string UFO features
-fformat1 string Features format
-fopenfile1 string Features file name
"-outfile" associated qualifiers
-odirectory2 string Output directory
"-treefile" associated qualifiers
-odirectory string Output directory
General qualifiers:
-auto boolean Turn off prompts
-stdout boolean Write first file to standard output
-filter boolean Read first file from standard input, write
first file to standard output
-options boolean Prompt for standard and additional values
-debug boolean Write debug output to program.dbg
-verbose boolean Report some/full command line options
-help boolean Report command line options and exit. More
information on associated and general
qualifiers can be found with -help -verbose
-warning boolean Report warnings
-error boolean Report errors
-fatal boolean Report fatal errors
-die boolean Report dying program messages
-version boolean Report version number and exit
|
| Qualifier | Type | Description | Allowed values | Default |
|---|---|---|---|---|
| Standard (Mandatory) qualifiers | ||||
| [-sequence] (Parameter 1) |
seqset | File containing a sequence alignment | Readable set of sequences | Required |
| -[no]besttree | toggle | Search for best tree | Toggle value Yes/No | Yes |
| -lengths | boolean | Use lengths from user trees | Boolean value Yes/No | No |
| -global | boolean | Global rearrangements | Boolean value Yes/No | No |
| -random | toggle | Randomize input order of species | Toggle value Yes/No | No |
| -randseed | integer | Random number seed (must be odd) | Any integer value | 3 |
| -randtimes | integer | How many times to randomise | Any integer value | 3 |
| -[no]basefrequency | toggle | Use empirical base frequencies | Toggle value Yes/No | Yes |
| -freqa | float | Frequency for A | Number from 0.000 to 1.000 | 0.25 |
| -freqc | float | Frequency for C | Number from 0.000 to 1.000 | 0.25 |
| -freqg | float | Frequency for G | Number from 0.000 to 1.000 | 0.25 |
| -freqt | float | Frequency for T/U | Number from 0.000 to 1.000 | 0.25 |
| -categories | toggle | More then one category of substitution rates | Toggle value Yes/No | No |
| -catnum | integer | Number of categories of substitution rates | Integer from 2 to 9 | 2 |
| -catvals | string | Space separated category values | Any string | |
| -catprob | string | Space separated probability values | Any string | |
| -deflambda | toggle | Use default mean block length of sites | Toggle value Yes/No | No |
| -lambda | float | Mean block length of sites having the same rate (greater than 1) | Number 1.000 or more | 1.0 |
| -og | toggle | Outgroup root | Toggle value Yes/No | No |
| -outgnum | integer | Number of the outgroup | Integer 1 or more | 1 |
| [-outfile] (Parameter 2) |
outfile | Phylip dnaml program output file | Output file | ednaml.outfile |
| [-[no]trout] (Parameter 3) |
toggle | Create a tree file | Toggle value Yes/No | Yes |
| -treefile | outfile | Phylip tree output file | Output file | ednaml.treefile |
| Additional (Optional) qualifiers | ||||
| -ttratio | float | Transition/transversion ratio | Any numeric value | 2.0 |
| -printdata | boolean | Print out the data at start of run | Boolean value Yes/No | No |
| -progress | boolean | Print indications of progress of run | Boolean value Yes/No | No |
| -[no]drawtree | boolean | Draw tree | Boolean value Yes/No | Yes |
| Advanced (Unprompted) qualifiers | ||||
| (none) | ||||
| Associated qualifiers | ||||
| "-sequence" associated seqset qualifiers | ||||
| -sbegin1 -sbegin_sequence |
integer | Start of each sequence to be used | Any integer value | 0 |
| -send1 -send_sequence |
integer | End of each sequence to be used | Any integer value | 0 |
| -sreverse1 -sreverse_sequence |
boolean | Reverse (if DNA) | Boolean value Yes/No | N |
| -sask1 -sask_sequence |
boolean | Ask for begin/end/reverse | Boolean value Yes/No | N |
| -snucleotide1 -snucleotide_sequence |
boolean | Sequence is nucleotide | Boolean value Yes/No | N |
| -sprotein1 -sprotein_sequence |
boolean | Sequence is protein | Boolean value Yes/No | N |
| -slower1 -slower_sequence |
boolean | Make lower case | Boolean value Yes/No | N |
| -supper1 -supper_sequence |
boolean | Make upper case | Boolean value Yes/No | N |
| -scircular1 -scircular_sequence |
boolean | Sequence is circular | Boolean value Yes/No | N |
| -squick1 -squick_sequence |
boolean | Read id and sequence only | Boolean value Yes/No | N |
| -sformat1 -sformat_sequence |
string | Input sequence format | Any string | |
| -iquery1 -iquery_sequence |
string | Input query fields or ID list | Any string | |
| -ioffset1 -ioffset_sequence |
integer | Input start position offset | Any integer value | 0 |
| -sdbname1 -sdbname_sequence |
string | Database name | Any string | |
| -sid1 -sid_sequence |
string | Entryname | Any string | |
| -ufo1 -ufo_sequence |
string | UFO features | Any string | |
| -fformat1 -fformat_sequence |
string | Features format | Any string | |
| -fopenfile1 -fopenfile_sequence |
string | Features file name | Any string | |
| "-outfile" associated outfile qualifiers | ||||
| -odirectory2 -odirectory_outfile |
string | Output directory | Any string | |
| "-treefile" associated outfile qualifiers | ||||
| -odirectory | string | Output directory | Any string | |
| General qualifiers | ||||
| -auto | boolean | Turn off prompts | Boolean value Yes/No | N |
| -stdout | boolean | Write first file to standard output | Boolean value Yes/No | N |
| -filter | boolean | Read first file from standard input, write first file to standard output | Boolean value Yes/No | N |
| -options | boolean | Prompt for standard and additional values | Boolean value Yes/No | N |
| -debug | boolean | Write debug output to program.dbg | Boolean value Yes/No | N |
| -verbose | boolean | Report some/full command line options | Boolean value Yes/No | Y |
| -help | boolean | Report command line options and exit. More information on associated and general qualifiers can be found with -help -verbose | Boolean value Yes/No | N |
| -warning | boolean | Report warnings | Boolean value Yes/No | Y |
| -error | boolean | Report errors | Boolean value Yes/No | Y |
| -fatal | boolean | Report fatal errors | Boolean value Yes/No | Y |
| -die | boolean | Report dying program messages | Boolean value Yes/No | Y |
| -version | boolean | Report version number and exit | Boolean value Yes/No | N |
5 13 Alpha AACGTGGCCAAAT Beta AAGGTCGCCAAAC Gamma CATTTCGTCACAA Delta GGTATTTCGGCCT Epsilon GGGATCTCGGCCC |
Nucleic acid sequence Maximum Likelihood method, version 3.57c.650
5 Species, 13 Sites
Name Sequences
---- ---------
Alpha AACGTGGCCAAAT
Beta AAGGTCGCCAAAC
Gamma CATTTCGTCACAA
Delta GGTATTTCGGCCT
Epsilon GGGATCTCGGCCC
Site category Rate of change Probability
1 1.000 0.400
2 3.200 0.600
Empirical Base Frequencies:
A 0.24615
C 0.29231
G 0.24615
T(U) 0.21538
Transition/transversion ratio = 2.000000
(Transition/transversion parameter = 1.523077)
+Epsilon
+-------------------------------------------------3
+--2 +----Delta
! !
! +Beta
!
--1--------------------Gamma
!
+---Alpha
remember: this is an unrooted tree!
Ln Likelihood = -71.83178
Examined 15 trees
Between And Length Approx. Confidence Limits
------- --- ------ ------- ---------- ------
1 2 0.08840 ( zero, 0.48704)
2 3 3.33144 ( zero, 12.20065) **
3 Epsilon 0.00003 ( zero, 0.39969)
3 Delta 0.30872 ( zero, 0.70442) **
2 Beta 0.00003 ( zero, 0.41145)
1 Gamma 1.35995 ( zero, 3.26620) **
1 Alpha 0.24900 ( zero, 0.68195)
* = significantly positive, P < 0.05
** = significantly positive, P < 0.01
Combination of categories that contributes the most to the likelihood:
2122121111 212
|
(((Epsilon:0.00012,Delta:0.30872):3.33144,Beta:0.00012):0.08840, Gamma:1.35995,Alpha:0.24900); |
| Program name | Description |
|---|---|
| distmat | Create a distance matrix from a multiple sequence alignment |
| ednacomp | DNA compatibility algorithm |
| ednadist | Nucleic acid sequence distance matrix program |
| ednainvar | Nucleic acid sequence invariants method |
| ednamlk | Phylogenies from nucleic acid maximum likelihood with clock |
| ednapars | DNA parsimony algorithm |
| ednapenny | Penny algorithm for DNA |
| eprotdist | Protein distance algorithm |
| eprotpars | Protein parsimony algorithm |
| erestml | Restriction site maximum likelihood method |
| eseqboot | Bootstrapped sequences algorithm |
| fdiscboot | Bootstrapped discrete sites algorithm |
| fdnacomp | DNA compatibility algorithm |
| fdnadist | Nucleic acid sequence distance matrix program |
| fdnainvar | Nucleic acid sequence invariants method |
| fdnaml | Estimate nucleotide phylogeny by maximum likelihood |
| fdnamlk | Estimates nucleotide phylogeny by maximum likelihood |
| fdnamove | Interactive DNA parsimony |
| fdnapars | DNA parsimony algorithm |
| fdnapenny | Penny algorithm for DNA |
| fdolmove | Interactive Dollo or polymorphism parsimony |
| ffreqboot | Bootstrapped genetic frequencies algorithm |
| fproml | Protein phylogeny by maximum likelihood |
| fpromlk | Protein phylogeny by maximum likelihood |
| fprotdist | Protein distance algorithm |
| fprotpars | Protein parsimony algorithm |
| frestboot | Bootstrapped restriction sites algorithm |
| frestdist | Calculate distance matrix from restriction sites or fragments |
| frestml | Restriction site maximum likelihood method |
| fseqboot | Bootstrapped sequences algorithm |
| fseqbootall | Bootstrapped sequences algorithm |
This application was modified for inclusion in EMBOSS by Ian Longden (il@sanger.ac.uk) Informatics Division, The Sanger Centre, Wellcome Trust Genome Campus, Hinxton, Cambridge, CB10 1SA, UK.
None