1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
|
seed used = 861385637
Before deleting alignment gaps
4 258
Gorilla_gorilla ATG GCC AGG TAC AGA TGC TGT CGC AGC CAG AGC CGC AGC AGA TGT TAC CGG CAG AGC CGG AGC AGG TGT TAC CGG CAG AGA CAA AGC CAG AGC CGG AGC AGA TGC TAC CGG CAG AGC CAA AGC CGG AGC AGG TGT TAC CGG CAG AGA CAA AGA AGT CGC AGA CGT AGG CGG AGG AGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGC TGC CGC CGC AGG TAC AGA CTG AGA CGT AGA AGA CCC TAT CAT ATT GTA TCT
Hylobates_lar ATG GCC AGG TAC AGA TGC TGC CGC AGC CAG AGC CGG AGC AGA TGT TAC CGC CAG AGC CGG AGC AGA TGT TAC CGC CAG AGG CAA AGC CAG AGT CGG AGC AGA TGT TAC CGC CAG AGC CAG AGC CGG AGC AGA TGT TAC CGC CAG AGA CAA AGA AGT CGG AGA CGA AGG AGG CGG AGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGT --- CGC CGC AGG TAC AGG CTG AGA CGT AGA AGC TGT TAC CAC ATT GTA TCT
Pan_troglodytes ATG GCC AGG TAC AGA TGC TGT CGC AGC CAG AGC CGG AGC AGA TGT TAC CGG CAG AGA CGG AGC AGG TGT TAC CGG CAA AGG CAA AGC CAA AGT CGG AGC AGA TGT TAC CGG CAG AGC CAG AGA CGG AGC AGG TGT TAC CGG CAA AGA CAA AGA AGT CGC AGA CGA AGG CGA CGG AGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGC TGC CGC CGC AGG TAC AGA CTG AGA CGT AAA AGA TGT TAC CAT ATT GTA TCT
Papio_cynocephalus ATG GCC AGG TAC AGA TGC TGC CGC AGC CAG AGC CGA AGC AGA TGC TAT CGC CAG AGC CGG AGC AGA TGT AAC CGC CAG AGA CAG AGC CAA AGC CGG AGA AGC TGC TAT CGC CAG AGC CAA AGC CGG AGC AGA TGT TAC CGC CAG AGA CAG AGA AGT CGT AGA CGA AGG AGG CGA CGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGC TTC CGC CGC AGG TAC AGG CTG AGG CGT AGG AGG CCC TAT CAC ATC GTG TCT
After deleting gaps. 255 sites
4 255
Gorilla_gorilla ATG GCC AGG TAC AGA TGC TGT CGC AGC CAG AGC CGC AGC AGA TGT TAC CGG CAG AGC CGG AGC AGG TGT TAC CGG CAG AGA CAA AGC CAG AGC CGG AGC AGA TGC TAC CGG CAG AGC CAA AGC CGG AGC AGG TGT TAC CGG CAG AGA CAA AGA AGT CGC AGA CGT AGG CGG AGG AGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGC CGC CGC AGG TAC AGA CTG AGA CGT AGA AGA CCC TAT CAT ATT GTA TCT
Hylobates_lar ATG GCC AGG TAC AGA TGC TGC CGC AGC CAG AGC CGG AGC AGA TGT TAC CGC CAG AGC CGG AGC AGA TGT TAC CGC CAG AGG CAA AGC CAG AGT CGG AGC AGA TGT TAC CGC CAG AGC CAG AGC CGG AGC AGA TGT TAC CGC CAG AGA CAA AGA AGT CGG AGA CGA AGG AGG CGG AGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGT CGC CGC AGG TAC AGG CTG AGA CGT AGA AGC TGT TAC CAC ATT GTA TCT
Pan_troglodytes ATG GCC AGG TAC AGA TGC TGT CGC AGC CAG AGC CGG AGC AGA TGT TAC CGG CAG AGA CGG AGC AGG TGT TAC CGG CAA AGG CAA AGC CAA AGT CGG AGC AGA TGT TAC CGG CAG AGC CAG AGA CGG AGC AGG TGT TAC CGG CAA AGA CAA AGA AGT CGC AGA CGA AGG CGA CGG AGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGC CGC CGC AGG TAC AGA CTG AGA CGT AAA AGA TGT TAC CAT ATT GTA TCT
Papio_cynocephalus ATG GCC AGG TAC AGA TGC TGC CGC AGC CAG AGC CGA AGC AGA TGC TAT CGC CAG AGC CGG AGC AGA TGT AAC CGC CAG AGA CAG AGC CAA AGC CGG AGA AGC TGC TAT CGC CAG AGC CAA AGC CGG AGC AGA TGT TAC CGC CAG AGA CAG AGA AGT CGT AGA CGA AGG AGG CGA CGC TGC CAG ACA CGG AGG AGA GCC ATG AGG TGC CGC CGC AGG TAC AGG CTG AGG CGT AGG AGG CCC TAT CAC ATC GTG TCT
Printing out site pattern counts
4 144 P
Gorilla_gorilla ACA AGA AGA AGA AGA AGA AGA AGA AGC AGC AGC AGC AGC AGG AGG AGG AGT ATG ATT CAA CAA CAG CAG CAG CAT CCC CGC CGC CGC CGG CGG CGG CGT CGT CTG GCC GTA TAC TAC TAC TAT TCT TGC TGC TGC TGT TGT TGT
Hylobates_lar ... ... ... ... ... ..C ..G ..G ... ... ... ... ..T ..A ... C.. ... ... ... ... ..G ... ... ... ..C TGT ... ..G ..G A.. ..C ... ..A ... ... ... ... ... ... ... ..C ... ... ..T ..T ..C ... ...
Pan_troglodytes ... .A. ... ... ... ... ... ..G ..A ... ... ... ..T ... ... C.. ... ... ... ... ..G ..A ..A ... ... TGT ... ... ..G ..A ... ... ..A ... ... ... ... ... ... ... ..C ... ... ... ..T ... ... ...
Papio_cynocephalus ... ..G ... ..C ..G ..G ..G ... ... ..A ... C.. ... ..A ... C.A ... ... ..C ..G ... ..A ... ... ..C ... ... ..T ..A A.. ..C ... ..A ... ... ... ..G A.. ... ..T ... ... ... ... ... ..C ..C ...
1 1 6 1 1 1 1 1 2 1 7 1 1 2 5
1 1 2 1 2 1 1 2 4 1 1 3 1 1 1
4 4 1 1 1 2 1 1 3 2 1 1 2 1 1
1 1 2
CODONML (in paml version 4.4, January 2010) algn
Model: One dN/dS ratio for branches
Codon frequency model: F3x4
Site-class models: NearlyNeutral
ns = 4 ls = 85
Codon usage in sequences
--------------------------------------------------------------------------------------
Phe TTT 0 0 0 0 | Ser TCT 1 1 1 1 | Tyr TAT 1 0 0 3 | Cys TGT 4 6 6 2
TTC 0 0 0 0 | TCC 0 0 0 0 | TAC 6 7 7 3 | TGC 4 3 3 6
Leu TTA 0 0 0 0 | TCA 0 0 0 0 | *** TAA 0 0 0 0 | *** TGA 0 0 0 0
TTG 0 0 0 0 | TCG 0 0 0 0 | TAG 0 0 0 0 | Trp TGG 0 0 0 0
--------------------------------------------------------------------------------------
Leu CTT 0 0 0 0 | Pro CCT 0 0 0 0 | His CAT 1 0 1 0 | Arg CGT 2 1 1 2
CTC 0 0 0 0 | CCC 1 0 0 1 | CAC 0 1 0 1 | CGC 5 7 4 8
CTA 0 0 0 0 | CCA 0 0 0 0 | Gln CAA 3 2 5 2 | CGA 0 1 2 3
CTG 1 1 1 1 | CCG 0 0 0 0 | CAG 7 8 5 8 | CGG 9 7 10 4
--------------------------------------------------------------------------------------
Ile ATT 1 1 1 0 | Thr ACT 0 0 0 0 | Asn AAT 0 0 0 0 | Ser AGT 1 2 2 1
ATC 0 0 0 1 | ACC 0 0 0 0 | AAC 0 0 0 1 | AGC 12 12 9 11
ATA 0 0 0 0 | ACA 1 1 1 1 | Lys AAA 0 0 1 0 | Arg AGA 12 11 12 10
Met ATG 2 2 2 2 | ACG 0 0 0 0 | AAG 0 0 0 0 | AGG 8 8 8 10
--------------------------------------------------------------------------------------
Val GTT 0 0 0 0 | Ala GCT 0 0 0 0 | Asp GAT 0 0 0 0 | Gly GGT 0 0 0 0
GTC 0 0 0 0 | GCC 2 2 2 2 | GAC 0 0 0 0 | GGC 0 0 0 0
GTA 1 1 1 0 | GCA 0 0 0 0 | Glu GAA 0 0 0 0 | GGA 0 0 0 0
GTG 0 0 0 1 | GCG 0 0 0 0 | GAG 0 0 0 0 | GGG 0 0 0 0
--------------------------------------------------------------------------------------
Codon position x base (3x4) table for each sequence.
#1: Gorilla_gorilla
position 1: T:0.18824 C:0.34118 A:0.43529 G:0.03529
position 2: T:0.05882 C:0.05882 A:0.21176 G:0.67059
position 3: T:0.12941 C:0.35294 A:0.20000 G:0.31765
Average T:0.12549 C:0.25098 A:0.28235 G:0.34118
#2: Hylobates_lar
position 1: T:0.20000 C:0.32941 A:0.43529 G:0.03529
position 2: T:0.05882 C:0.04706 A:0.21176 G:0.68235
position 3: T:0.12941 C:0.37647 A:0.18824 G:0.30588
Average T:0.12941 C:0.25098 A:0.27843 G:0.34118
#3: Pan_troglodytes
position 1: T:0.20000 C:0.34118 A:0.42353 G:0.03529
position 2: T:0.05882 C:0.04706 A:0.22353 G:0.67059
position 3: T:0.14118 C:0.29412 A:0.25882 G:0.30588
Average T:0.13333 C:0.22745 A:0.30196 G:0.33725
#4: Papio_cynocephalus
position 1: T:0.17647 C:0.35294 A:0.43529 G:0.03529
position 2: T:0.05882 C:0.05882 A:0.21176 G:0.67059
position 3: T:0.10588 C:0.40000 A:0.18824 G:0.30588
Average T:0.11373 C:0.27059 A:0.27843 G:0.33725
Sums of codon usage counts
------------------------------------------------------------------------------
Phe F TTT 0 | Ser S TCT 4 | Tyr Y TAT 4 | Cys C TGT 18
TTC 0 | TCC 0 | TAC 23 | TGC 16
Leu L TTA 0 | TCA 0 | *** * TAA 0 | *** * TGA 0
TTG 0 | TCG 0 | TAG 0 | Trp W TGG 0
------------------------------------------------------------------------------
Leu L CTT 0 | Pro P CCT 0 | His H CAT 2 | Arg R CGT 6
CTC 0 | CCC 2 | CAC 2 | CGC 24
CTA 0 | CCA 0 | Gln Q CAA 12 | CGA 6
CTG 4 | CCG 0 | CAG 28 | CGG 30
------------------------------------------------------------------------------
Ile I ATT 3 | Thr T ACT 0 | Asn N AAT 0 | Ser S AGT 6
ATC 1 | ACC 0 | AAC 1 | AGC 44
ATA 0 | ACA 4 | Lys K AAA 1 | Arg R AGA 45
Met M ATG 8 | ACG 0 | AAG 0 | AGG 34
------------------------------------------------------------------------------
Val V GTT 0 | Ala A GCT 0 | Asp D GAT 0 | Gly G GGT 0
GTC 0 | GCC 8 | GAC 0 | GGC 0
GTA 3 | GCA 0 | Glu E GAA 0 | GGA 0
GTG 1 | GCG 0 | GAG 0 | GGG 0
------------------------------------------------------------------------------
Codon position x base (3x4) table, overall
position 1: T:0.19118 C:0.34118 A:0.43235 G:0.03529
position 2: T:0.05882 C:0.05294 A:0.21471 G:0.67353
position 3: T:0.12647 C:0.35588 A:0.20882 G:0.30882
Average T:0.12549 C:0.25000 A:0.28529 G:0.33922
Codon frequencies under model, for use in evolver (TTT TTC TTA TTG ... GGG):
0.00149417 0.00420453 0.00246712 0.00364856
0.00134475 0.00378408 0.00222041 0.00328370
0.00545372 0.01534653 0.00000000 0.00000000
0.01710826 0.04814185 0.00000000 0.04177598
0.00266652 0.00750347 0.00440286 0.00651127
0.00239987 0.00675312 0.00396257 0.00586014
0.00973280 0.02738765 0.01607044 0.02376614
0.03053166 0.08591468 0.05041275 0.07455406
0.00337913 0.00950870 0.00557949 0.00825135
0.00304121 0.00855783 0.00502154 0.00742622
0.01233381 0.03470676 0.02036513 0.03011744
0.03869099 0.10887464 0.06388512 0.09447799
0.00027585 0.00077622 0.00045547 0.00067358
0.00024826 0.00069860 0.00040992 0.00060622
0.00100684 0.00283321 0.00166246 0.00245857
0.00315845 0.00888773 0.00521511 0.00771249
Nei & Gojobori 1986. dN/dS (dN, dS)
(Note: This matrix is not used in later ML. analysis.
Use runmode = -2 for ML pairwise comparison.)
Gorilla_gorilla
Hylobates_lar 0.0292 (0.0152 0.5209)
Pan_troglodytes 0.0952 (0.0257 0.2703) 0.0576 (0.0204 0.3547)
Papio_cynocephalus 0.0287 (0.0204 0.7121) 0.0740 (0.0360 0.4864) 0.0458 (0.0469 1.0261)
TREE # 1: ((2, (1, 3)), 4); MP score: 53
This is a rooted tree. Please check!
lnL(ntime: 6 np: 9): -511.919210 +0.000000
5..6 6..2 6..7 7..1 7..3 5..4
0.223828 0.032940 0.166566 0.156609 0.096748 0.230835 3.542367 0.938133 0.060763
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).
tree length = 0.90753
((2: 0.032940, (1: 0.156609, 3: 0.096748): 0.166566): 0.223828, 4: 0.230835);
((Hylobates_lar: 0.032940, (Gorilla_gorilla: 0.156609, Pan_troglodytes: 0.096748): 0.166566): 0.223828, Papio_cynocephalus: 0.230835);
Detailed output identifying parameters
kappa (ts/tv) = 3.54237
dN/dS (w) for site classes (K=2)
p: 0.93813 0.06187
w: 0.06076 1.00000
dN & dS for each branch
branch t N S dN/dS dN dS N*dN S*dS
5..6 0.224 159.3 95.7 0.1189 0.0197 0.1659 3.1 15.9
6..2 0.033 159.3 95.7 0.1189 0.0029 0.0244 0.5 2.3
6..7 0.167 159.3 95.7 0.1189 0.0147 0.1235 2.3 11.8
7..1 0.157 159.3 95.7 0.1189 0.0138 0.1161 2.2 11.1
7..3 0.097 159.3 95.7 0.1189 0.0085 0.0717 1.4 6.9
5..4 0.231 159.3 95.7 0.1189 0.0203 0.1711 3.2 16.4
Time used: 0:03
|
