1#NEXUS
2
3BEGIN TAXA;
4	DIMENSIONS NTAX = 13;
5	TAXLABELS
6		'U68496' 'U68497' 'U68498' 'U68499' 'U68500' 'U68501' 'U68502' 'U68503' 'U68504' 'U68505' 'U68506' 'U68507' 'U68508' ;
7END;
8
9BEGIN CHARACTERS;
10	DIMENSIONS NCHAR = 273;
11	FORMAT
12		DATATYPE = DNA
13		GAP=-
14		MISSING=?
15	;
16
17MATRIX
18	'U68496' GTAGTAATTAGATCTGAAAACTTCTCGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAGAAGTATACATTTTGGACCAGGGAAAGCATTTTATGCAGGAGAAATAATAGGAGATATAAGACAAGCATATTGTACCCTTAATGGAACAGAATGGAATAACACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTATTAAAACAATAGTTTTTAATCAATCC
19	'U68497' ATAGTAATTAGATCTGAAAACTTCTCGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAGAAGTATACATTTCGGACCAGGGAAAGCATTTTATGCAGGAGAAATAATAGGAGATATAAGACAAGCATATTGTACTCTTAATGGAGCAGAATGGAATAACACTGTAAAACAGGTAGCTGCAAAATTAAGAGAACAATTTAATAAAACAATAATCTTTAATCAATCC
20	'U68498' GTAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTGTAAGACCCGGCAACAATACAAGAAGAAGTATACATATAGGACCAGGGAGAGCATATTATACAGGAGAAGTAATAGGAGATATAAGACAAGCACATTGTAACCTTAGTAGAACAGACTGGAATAAAACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTAATACAACAATAGTCTTTAATCAATCC
21	'U68499' ATAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAAGTGTGAAAGACCCAACAACAATACAAGAAAAAGTGTACATATAGGACCAGGGAGAGCATATTATACAGGAGAAATAATAGGAGATATAAGACAAGCACATTGTAACCTTAGTGGAACAGAATGGAGGGAAACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTAATAAAACAATAGTCTTTAATCAATCC
22	'U68500' ATAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACATCTAAATGAATCTGTAGAAATTATTTGTGAAAGACCCAACAACAATACAAGAAAAAGTGTACATATGGGACCAGGGAGAGCATATTACACAGGAGAAATAATAGGAGATATAAGACAAGCACATTGTAACATTAGTAGAACAAATTGGACGGAAACTTTAAAACAGGTAGCTGAAAAATTAAGAGAACAATTTAATAAAACAATAGTCTTTAATCAATCC
23	'U68501' GTAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAGACCATAATAGTACAGCTAAATAAACCTGTAAAAATTAATTGTACAAGACCCAACAACAATGCAAAAATAAGAATACATATAGGACCAGGGAGACCATTTTATACAGCAGGAGAAATAGGAAATATAAGACAAGCACATTGTAACCTTAGTAGAACAGACTGGAATAACACTTTAAAACTGGTAGCTGAAAAATTAAGAGAACAATTTAATAAAACAATAGTCTTTAATCAATCC
24	'U68502' GTAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAGACCATAATAGTACAGCTAAATAACTCTGTAGCAATTAAGTGTGAAAGACCCAACAACAATACAAGAAAAAGTATACCTATAGGACCAGGGAGAGCCTTTTATACAACAGGAGACATAGGAGATATAAGACAAGCACATTGTAACCTTAGTAGAAAAGACTGGAATGACACTTTAAGACAGGTAGTTGGAAAGTTAAGAGAACAATTTGGAAGAACAATAATCTTTAATCAATCC
25	'U68503' ATAGTAATTAGATCTGAAAACTTCACGAACAATGCTAAAACCATAATAGTACAGCTAAAGGAACCTGTAGACATTACTTGTGAAAGACCCAGCAACAATACAAGAAAAAGTATACATATAGGACCAGGAAAAGCATTTTATGCAACAGGAGAAATAGGAGATATAAGACGAGCACATTGTAACCTTAATAGAACAGCATGGAATAAAACTTTAAAACAGGTAGTTGAAAAATTAAGAGAACAATTTAAGAAAACAATAACCTTTAACCAATCC
26	'U68504' ATAGTAATTAGATCTGAAAACTTCACAGACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAAAAGTATACATATAGGACCAGGGAGAGCTTTTTATACAACAGGTGAAATAGGAGATTTAAGACAAGCACATTGTAACCTTAGTAGAACAGCATGGAATGAAACTTTAAGACAAGTAGCTCAAAAATTAAAAGAACAATTTAATAGAACAATAGTTTTTAATCAATCC
27	'U68505' ATAGTAATTAGATCTGAAAACTTCACAGACAATGCTAAAACCATAATAGTACAGCTAAACAAATCTGTAGAAATTAATTGTATAAGACCCAACAACAATACAAGAAAAAGTATACATATAGGACCAGGGAGAGCATTTTATACAACAGGAGACATAGGAGATATAAGACAAGCATATTGTAACCTTAGTAGAACAGCATGGAATGAAACTTTAAGACAAGTAGCTCAAAAATTAAAAGAACAATTTAATAGAACAATAGTCTTTAATCAATCC
28	'U68506' ATAGTAATTAGATCTGAAAACTTCACGGACAATGCTAAAACCATAATAGTACAGCTAAATAAATCTGTAGAAATTAATTGTACAAGACCCAACAACAATACAAGAAAAAGTATAAATTTTGGACCAGGGAGAGCATTTTATACAACAGGTGAAATAGGAAATTTAAGACAAGCACATTGTAACCTTAGTAGAACAGCATGGAATGAAACTTTAAGACAAGTAGCTAAAAAACTAAAAGAACAATTTAATAGAACAATAGTTTTTAATCAATCC
29	'U68507' GTAGTAATTCGATCTGAAAACTTCACGGACAATGCTAAAACCATAATAGTACAGCTAAACAAATCTGTAGAAATTACTTGTGTAAGACCCAACAACAATACAAGAAAAAGTATAAATATAAGACCAGGGAGAGCATTTTATACAACAGGAGAAATAGGAGATATAAGACAAGCACATTGTAACCTTAGTAGAACAGCATGGAATGAAGCTTTAAGACAAGTAGCTAAAAAATTAAAAGAACAATTTAATAGAACAATAGTCTTTAATCAATCC
30	'U68508' ATAGTAATCAGATCTGAAAACTTCTCGGACAATGCTAAAACCATAATAGTACAGCTAAACAACACTGTAAACATTACTTGTGAAAGACCCAACAACAATACAAGAAAAAGGATACATATAGGACCAGGGAGAGCAGTTTATACAACAGGACAAATAGGAGATATAAGAAAAGCACATTGTAACCTTAGTAGAACAAATTGGACTGAAACTTTAAGACAAGTAGCTGAAAAATTAAAAGAACAATTTAATAAAACAATAATCTTTAATAATTCC;
31END;
32
33BEGIN TREES;
34	TREE tree = (U68496: 0.023746, U68497: 0.079178, ((U68498: 0.050407, (U68499: 0.028429, U68500: 0.099890): 0.083460): 0.023647, (U68501: 0.131681, (U68502: 0.191150, (U68503: 0.193336, ((((U68504: 0.018984, U68506: 0.070370): 0.055603, U68505: 0.030301): 0.036688, U68507: 0.061120): 0.034773, U68508: 0.198346): 0.050591): 0.064013): 0.033695): 0.048101): 0.152813);
35END;
36
37BEGIN HYPHY;
38
39
40global P2=0.7104121522639678;
41P2:<1;
42global W_1=0.05768869986246565;
43W_1:<1;
44global W_2=3.641844577696067;
45W_2:>1;
46global P1=0.3739141342952348;
47P1:<1;
48global kappa=0.3597183898000834;
49
50c.weights={1,3};
51c.weights[0][0]:=P1;
52c.weights[0][1]:=(1-P1)*P2;
53c.weights[0][2]:=(1-P1)*(1-P2);
54
55
56c.points={1,3};
57c.points[0][0]:=W_1;
58c.points[0][1]:=1;
59c.points[0][2]:=W_2;
60
61
62category c=(3,c.weights,MEAN,,c.points,0,1,);
63
64modelMatrix={61,61};
65modelMatrix[0][1]:=kappa*c*t;
66modelMatrix[0][2]:=t;
67modelMatrix[0][3]:=kappa*c*t;
68modelMatrix[0][4]:=kappa*c*t;
69modelMatrix[0][8]:=c*t;
70modelMatrix[0][12]:=kappa*c*t;
71modelMatrix[0][16]:=kappa*c*t;
72modelMatrix[0][32]:=c*t;
73modelMatrix[1][0]:=kappa*c*t;
74modelMatrix[1][2]:=kappa*c*t;
75modelMatrix[1][3]:=t;
76modelMatrix[1][5]:=kappa*c*t;
77modelMatrix[1][9]:=c*t;
78modelMatrix[1][13]:=kappa*c*t;
79modelMatrix[1][17]:=kappa*c*t;
80modelMatrix[1][33]:=c*t;
81modelMatrix[1][48]:=kappa*c*t;
82modelMatrix[2][0]:=t;
83modelMatrix[2][1]:=kappa*c*t;
84modelMatrix[2][3]:=kappa*c*t;
85modelMatrix[2][6]:=kappa*c*t;
86modelMatrix[2][10]:=c*t;
87modelMatrix[2][14]:=kappa*c*t;
88modelMatrix[2][18]:=kappa*c*t;
89modelMatrix[2][34]:=c*t;
90modelMatrix[3][0]:=kappa*c*t;
91modelMatrix[3][1]:=t;
92modelMatrix[3][2]:=kappa*c*t;
93modelMatrix[3][7]:=kappa*c*t;
94modelMatrix[3][11]:=c*t;
95modelMatrix[3][15]:=kappa*c*t;
96modelMatrix[3][19]:=kappa*c*t;
97modelMatrix[3][35]:=c*t;
98modelMatrix[3][49]:=kappa*c*t;
99modelMatrix[4][0]:=kappa*c*t;
100modelMatrix[4][5]:=kappa*t;
101modelMatrix[4][6]:=t;
102modelMatrix[4][7]:=kappa*t;
103modelMatrix[4][8]:=kappa*c*t;
104modelMatrix[4][12]:=c*t;
105modelMatrix[4][20]:=kappa*c*t;
106modelMatrix[4][36]:=c*t;
107modelMatrix[4][50]:=kappa*c*t;
108modelMatrix[5][1]:=kappa*c*t;
109modelMatrix[5][4]:=kappa*t;
110modelMatrix[5][6]:=kappa*t;
111modelMatrix[5][7]:=t;
112modelMatrix[5][9]:=kappa*c*t;
113modelMatrix[5][13]:=c*t;
114modelMatrix[5][21]:=kappa*c*t;
115modelMatrix[5][37]:=c*t;
116modelMatrix[5][51]:=kappa*c*t;
117modelMatrix[6][2]:=kappa*c*t;
118modelMatrix[6][4]:=t;
119modelMatrix[6][5]:=kappa*t;
120modelMatrix[6][7]:=kappa*t;
121modelMatrix[6][10]:=kappa*c*t;
122modelMatrix[6][14]:=c*t;
123modelMatrix[6][22]:=kappa*c*t;
124modelMatrix[6][38]:=c*t;
125modelMatrix[6][52]:=kappa*c*t;
126modelMatrix[7][3]:=kappa*c*t;
127modelMatrix[7][4]:=kappa*t;
128modelMatrix[7][5]:=t;
129modelMatrix[7][6]:=kappa*t;
130modelMatrix[7][11]:=kappa*c*t;
131modelMatrix[7][15]:=c*t;
132modelMatrix[7][23]:=kappa*c*t;
133modelMatrix[7][39]:=c*t;
134modelMatrix[7][53]:=kappa*c*t;
135modelMatrix[8][0]:=c*t;
136modelMatrix[8][4]:=kappa*c*t;
137modelMatrix[8][9]:=kappa*c*t;
138modelMatrix[8][10]:=t;
139modelMatrix[8][11]:=kappa*c*t;
140modelMatrix[8][12]:=kappa*c*t;
141modelMatrix[8][24]:=kappa*t;
142modelMatrix[8][40]:=c*t;
143modelMatrix[9][1]:=c*t;
144modelMatrix[9][5]:=kappa*c*t;
145modelMatrix[9][8]:=kappa*c*t;
146modelMatrix[9][10]:=kappa*c*t;
147modelMatrix[9][11]:=t;
148modelMatrix[9][13]:=kappa*c*t;
149modelMatrix[9][25]:=kappa*c*t;
150modelMatrix[9][41]:=c*t;
151modelMatrix[9][54]:=kappa*c*t;
152modelMatrix[10][2]:=c*t;
153modelMatrix[10][6]:=kappa*c*t;
154modelMatrix[10][8]:=t;
155modelMatrix[10][9]:=kappa*c*t;
156modelMatrix[10][11]:=kappa*c*t;
157modelMatrix[10][14]:=kappa*c*t;
158modelMatrix[10][26]:=kappa*t;
159modelMatrix[10][42]:=c*t;
160modelMatrix[10][55]:=kappa*c*t;
161modelMatrix[11][3]:=c*t;
162modelMatrix[11][7]:=kappa*c*t;
163modelMatrix[11][8]:=kappa*c*t;
164modelMatrix[11][9]:=t;
165modelMatrix[11][10]:=kappa*c*t;
166modelMatrix[11][15]:=kappa*c*t;
167modelMatrix[11][27]:=kappa*c*t;
168modelMatrix[11][43]:=c*t;
169modelMatrix[11][56]:=kappa*c*t;
170modelMatrix[12][0]:=kappa*c*t;
171modelMatrix[12][4]:=c*t;
172modelMatrix[12][8]:=kappa*c*t;
173modelMatrix[12][13]:=kappa*t;
174modelMatrix[12][14]:=c*t;
175modelMatrix[12][15]:=kappa*t;
176modelMatrix[12][28]:=kappa*c*t;
177modelMatrix[12][44]:=c*t;
178modelMatrix[12][57]:=kappa*c*t;
179modelMatrix[13][1]:=kappa*c*t;
180modelMatrix[13][5]:=c*t;
181modelMatrix[13][9]:=kappa*c*t;
182modelMatrix[13][12]:=kappa*t;
183modelMatrix[13][14]:=kappa*c*t;
184modelMatrix[13][15]:=t;
185modelMatrix[13][29]:=kappa*c*t;
186modelMatrix[13][45]:=c*t;
187modelMatrix[13][58]:=kappa*c*t;
188modelMatrix[14][2]:=kappa*c*t;
189modelMatrix[14][6]:=c*t;
190modelMatrix[14][10]:=kappa*c*t;
191modelMatrix[14][12]:=c*t;
192modelMatrix[14][13]:=kappa*c*t;
193modelMatrix[14][15]:=kappa*c*t;
194modelMatrix[14][30]:=kappa*c*t;
195modelMatrix[14][46]:=c*t;
196modelMatrix[14][59]:=kappa*c*t;
197modelMatrix[15][3]:=kappa*c*t;
198modelMatrix[15][7]:=c*t;
199modelMatrix[15][11]:=kappa*c*t;
200modelMatrix[15][12]:=kappa*t;
201modelMatrix[15][13]:=t;
202modelMatrix[15][14]:=kappa*c*t;
203modelMatrix[15][31]:=kappa*c*t;
204modelMatrix[15][47]:=c*t;
205modelMatrix[15][60]:=kappa*c*t;
206modelMatrix[16][0]:=kappa*c*t;
207modelMatrix[16][17]:=kappa*c*t;
208modelMatrix[16][18]:=t;
209modelMatrix[16][19]:=kappa*c*t;
210modelMatrix[16][20]:=kappa*c*t;
211modelMatrix[16][24]:=c*t;
212modelMatrix[16][28]:=kappa*c*t;
213modelMatrix[16][32]:=kappa*c*t;
214modelMatrix[17][1]:=kappa*c*t;
215modelMatrix[17][16]:=kappa*c*t;
216modelMatrix[17][18]:=kappa*c*t;
217modelMatrix[17][19]:=t;
218modelMatrix[17][21]:=kappa*c*t;
219modelMatrix[17][25]:=c*t;
220modelMatrix[17][29]:=kappa*c*t;
221modelMatrix[17][33]:=kappa*c*t;
222modelMatrix[17][48]:=c*t;
223modelMatrix[18][2]:=kappa*c*t;
224modelMatrix[18][16]:=t;
225modelMatrix[18][17]:=kappa*c*t;
226modelMatrix[18][19]:=kappa*c*t;
227modelMatrix[18][22]:=kappa*c*t;
228modelMatrix[18][26]:=c*t;
229modelMatrix[18][30]:=kappa*c*t;
230modelMatrix[18][34]:=kappa*c*t;
231modelMatrix[19][3]:=kappa*c*t;
232modelMatrix[19][16]:=kappa*c*t;
233modelMatrix[19][17]:=t;
234modelMatrix[19][18]:=kappa*c*t;
235modelMatrix[19][23]:=kappa*c*t;
236modelMatrix[19][27]:=c*t;
237modelMatrix[19][31]:=kappa*c*t;
238modelMatrix[19][35]:=kappa*c*t;
239modelMatrix[19][49]:=c*t;
240modelMatrix[20][4]:=kappa*c*t;
241modelMatrix[20][16]:=kappa*c*t;
242modelMatrix[20][21]:=kappa*t;
243modelMatrix[20][22]:=t;
244modelMatrix[20][23]:=kappa*t;
245modelMatrix[20][24]:=kappa*c*t;
246modelMatrix[20][28]:=c*t;
247modelMatrix[20][36]:=kappa*c*t;
248modelMatrix[20][50]:=c*t;
249modelMatrix[21][5]:=kappa*c*t;
250modelMatrix[21][17]:=kappa*c*t;
251modelMatrix[21][20]:=kappa*t;
252modelMatrix[21][22]:=kappa*t;
253modelMatrix[21][23]:=t;
254modelMatrix[21][25]:=kappa*c*t;
255modelMatrix[21][29]:=c*t;
256modelMatrix[21][37]:=kappa*c*t;
257modelMatrix[21][51]:=c*t;
258modelMatrix[22][6]:=kappa*c*t;
259modelMatrix[22][18]:=kappa*c*t;
260modelMatrix[22][20]:=t;
261modelMatrix[22][21]:=kappa*t;
262modelMatrix[22][23]:=kappa*t;
263modelMatrix[22][26]:=kappa*c*t;
264modelMatrix[22][30]:=c*t;
265modelMatrix[22][38]:=kappa*c*t;
266modelMatrix[22][52]:=c*t;
267modelMatrix[23][7]:=kappa*c*t;
268modelMatrix[23][19]:=kappa*c*t;
269modelMatrix[23][20]:=kappa*t;
270modelMatrix[23][21]:=t;
271modelMatrix[23][22]:=kappa*t;
272modelMatrix[23][27]:=kappa*c*t;
273modelMatrix[23][31]:=c*t;
274modelMatrix[23][39]:=kappa*c*t;
275modelMatrix[23][53]:=c*t;
276modelMatrix[24][8]:=kappa*t;
277modelMatrix[24][16]:=c*t;
278modelMatrix[24][20]:=kappa*c*t;
279modelMatrix[24][25]:=kappa*t;
280modelMatrix[24][26]:=t;
281modelMatrix[24][27]:=kappa*t;
282modelMatrix[24][28]:=kappa*c*t;
283modelMatrix[24][40]:=kappa*c*t;
284modelMatrix[25][9]:=kappa*c*t;
285modelMatrix[25][17]:=c*t;
286modelMatrix[25][21]:=kappa*c*t;
287modelMatrix[25][24]:=kappa*t;
288modelMatrix[25][26]:=kappa*t;
289modelMatrix[25][27]:=t;
290modelMatrix[25][29]:=kappa*c*t;
291modelMatrix[25][41]:=kappa*c*t;
292modelMatrix[25][54]:=c*t;
293modelMatrix[26][10]:=kappa*t;
294modelMatrix[26][18]:=c*t;
295modelMatrix[26][22]:=kappa*c*t;
296modelMatrix[26][24]:=t;
297modelMatrix[26][25]:=kappa*t;
298modelMatrix[26][27]:=kappa*t;
299modelMatrix[26][30]:=kappa*c*t;
300modelMatrix[26][42]:=kappa*c*t;
301modelMatrix[26][55]:=c*t;
302modelMatrix[27][11]:=kappa*c*t;
303modelMatrix[27][19]:=c*t;
304modelMatrix[27][23]:=kappa*c*t;
305modelMatrix[27][24]:=kappa*t;
306modelMatrix[27][25]:=t;
307modelMatrix[27][26]:=kappa*t;
308modelMatrix[27][31]:=kappa*c*t;
309modelMatrix[27][43]:=kappa*c*t;
310modelMatrix[27][56]:=c*t;
311modelMatrix[28][12]:=kappa*c*t;
312modelMatrix[28][16]:=kappa*c*t;
313modelMatrix[28][20]:=c*t;
314modelMatrix[28][24]:=kappa*c*t;
315modelMatrix[28][29]:=kappa*t;
316modelMatrix[28][30]:=t;
317modelMatrix[28][31]:=kappa*t;
318modelMatrix[28][44]:=kappa*c*t;
319modelMatrix[28][57]:=t;
320modelMatrix[29][13]:=kappa*c*t;
321modelMatrix[29][17]:=kappa*c*t;
322modelMatrix[29][21]:=c*t;
323modelMatrix[29][25]:=kappa*c*t;
324modelMatrix[29][28]:=kappa*t;
325modelMatrix[29][30]:=kappa*t;
326modelMatrix[29][31]:=t;
327modelMatrix[29][45]:=kappa*c*t;
328modelMatrix[29][58]:=c*t;
329modelMatrix[30][14]:=kappa*c*t;
330modelMatrix[30][18]:=kappa*c*t;
331modelMatrix[30][22]:=c*t;
332modelMatrix[30][26]:=kappa*c*t;
333modelMatrix[30][28]:=t;
334modelMatrix[30][29]:=kappa*t;
335modelMatrix[30][31]:=kappa*t;
336modelMatrix[30][46]:=kappa*c*t;
337modelMatrix[30][59]:=t;
338modelMatrix[31][15]:=kappa*c*t;
339modelMatrix[31][19]:=kappa*c*t;
340modelMatrix[31][23]:=c*t;
341modelMatrix[31][27]:=kappa*c*t;
342modelMatrix[31][28]:=kappa*t;
343modelMatrix[31][29]:=t;
344modelMatrix[31][30]:=kappa*t;
345modelMatrix[31][47]:=kappa*c*t;
346modelMatrix[31][60]:=c*t;
347modelMatrix[32][0]:=c*t;
348modelMatrix[32][16]:=kappa*c*t;
349modelMatrix[32][33]:=kappa*c*t;
350modelMatrix[32][34]:=t;
351modelMatrix[32][35]:=kappa*c*t;
352modelMatrix[32][36]:=kappa*c*t;
353modelMatrix[32][40]:=c*t;
354modelMatrix[32][44]:=kappa*c*t;
355modelMatrix[33][1]:=c*t;
356modelMatrix[33][17]:=kappa*c*t;
357modelMatrix[33][32]:=kappa*c*t;
358modelMatrix[33][34]:=kappa*c*t;
359modelMatrix[33][35]:=t;
360modelMatrix[33][37]:=kappa*c*t;
361modelMatrix[33][41]:=c*t;
362modelMatrix[33][45]:=kappa*c*t;
363modelMatrix[33][48]:=kappa*c*t;
364modelMatrix[34][2]:=c*t;
365modelMatrix[34][18]:=kappa*c*t;
366modelMatrix[34][32]:=t;
367modelMatrix[34][33]:=kappa*c*t;
368modelMatrix[34][35]:=kappa*c*t;
369modelMatrix[34][38]:=kappa*c*t;
370modelMatrix[34][42]:=c*t;
371modelMatrix[34][46]:=kappa*c*t;
372modelMatrix[35][3]:=c*t;
373modelMatrix[35][19]:=kappa*c*t;
374modelMatrix[35][32]:=kappa*c*t;
375modelMatrix[35][33]:=t;
376modelMatrix[35][34]:=kappa*c*t;
377modelMatrix[35][39]:=kappa*c*t;
378modelMatrix[35][43]:=c*t;
379modelMatrix[35][47]:=kappa*c*t;
380modelMatrix[35][49]:=kappa*c*t;
381modelMatrix[36][4]:=c*t;
382modelMatrix[36][20]:=kappa*c*t;
383modelMatrix[36][32]:=kappa*c*t;
384modelMatrix[36][37]:=kappa*t;
385modelMatrix[36][38]:=t;
386modelMatrix[36][39]:=kappa*t;
387modelMatrix[36][40]:=kappa*c*t;
388modelMatrix[36][44]:=c*t;
389modelMatrix[36][50]:=kappa*c*t;
390modelMatrix[37][5]:=c*t;
391modelMatrix[37][21]:=kappa*c*t;
392modelMatrix[37][33]:=kappa*c*t;
393modelMatrix[37][36]:=kappa*t;
394modelMatrix[37][38]:=kappa*t;
395modelMatrix[37][39]:=t;
396modelMatrix[37][41]:=kappa*c*t;
397modelMatrix[37][45]:=c*t;
398modelMatrix[37][51]:=kappa*c*t;
399modelMatrix[38][6]:=c*t;
400modelMatrix[38][22]:=kappa*c*t;
401modelMatrix[38][34]:=kappa*c*t;
402modelMatrix[38][36]:=t;
403modelMatrix[38][37]:=kappa*t;
404modelMatrix[38][39]:=kappa*t;
405modelMatrix[38][42]:=kappa*c*t;
406modelMatrix[38][46]:=c*t;
407modelMatrix[38][52]:=kappa*c*t;
408modelMatrix[39][7]:=c*t;
409modelMatrix[39][23]:=kappa*c*t;
410modelMatrix[39][35]:=kappa*c*t;
411modelMatrix[39][36]:=kappa*t;
412modelMatrix[39][37]:=t;
413modelMatrix[39][38]:=kappa*t;
414modelMatrix[39][43]:=kappa*c*t;
415modelMatrix[39][47]:=c*t;
416modelMatrix[39][53]:=kappa*c*t;
417modelMatrix[40][8]:=c*t;
418modelMatrix[40][24]:=kappa*c*t;
419modelMatrix[40][32]:=c*t;
420modelMatrix[40][36]:=kappa*c*t;
421modelMatrix[40][41]:=kappa*t;
422modelMatrix[40][42]:=t;
423modelMatrix[40][43]:=kappa*t;
424modelMatrix[40][44]:=kappa*c*t;
425modelMatrix[41][9]:=c*t;
426modelMatrix[41][25]:=kappa*c*t;
427modelMatrix[41][33]:=c*t;
428modelMatrix[41][37]:=kappa*c*t;
429modelMatrix[41][40]:=kappa*t;
430modelMatrix[41][42]:=kappa*t;
431modelMatrix[41][43]:=t;
432modelMatrix[41][45]:=kappa*c*t;
433modelMatrix[41][54]:=kappa*c*t;
434modelMatrix[42][10]:=c*t;
435modelMatrix[42][26]:=kappa*c*t;
436modelMatrix[42][34]:=c*t;
437modelMatrix[42][38]:=kappa*c*t;
438modelMatrix[42][40]:=t;
439modelMatrix[42][41]:=kappa*t;
440modelMatrix[42][43]:=kappa*t;
441modelMatrix[42][46]:=kappa*c*t;
442modelMatrix[42][55]:=kappa*c*t;
443modelMatrix[43][11]:=c*t;
444modelMatrix[43][27]:=kappa*c*t;
445modelMatrix[43][35]:=c*t;
446modelMatrix[43][39]:=kappa*c*t;
447modelMatrix[43][40]:=kappa*t;
448modelMatrix[43][41]:=t;
449modelMatrix[43][42]:=kappa*t;
450modelMatrix[43][47]:=kappa*c*t;
451modelMatrix[43][56]:=kappa*c*t;
452modelMatrix[44][12]:=c*t;
453modelMatrix[44][28]:=kappa*c*t;
454modelMatrix[44][32]:=kappa*c*t;
455modelMatrix[44][36]:=c*t;
456modelMatrix[44][40]:=kappa*c*t;
457modelMatrix[44][45]:=kappa*t;
458modelMatrix[44][46]:=t;
459modelMatrix[44][47]:=kappa*t;
460modelMatrix[44][57]:=kappa*c*t;
461modelMatrix[45][13]:=c*t;
462modelMatrix[45][29]:=kappa*c*t;
463modelMatrix[45][33]:=kappa*c*t;
464modelMatrix[45][37]:=c*t;
465modelMatrix[45][41]:=kappa*c*t;
466modelMatrix[45][44]:=kappa*t;
467modelMatrix[45][46]:=kappa*t;
468modelMatrix[45][47]:=t;
469modelMatrix[45][58]:=kappa*c*t;
470modelMatrix[46][14]:=c*t;
471modelMatrix[46][30]:=kappa*c*t;
472modelMatrix[46][34]:=kappa*c*t;
473modelMatrix[46][38]:=c*t;
474modelMatrix[46][42]:=kappa*c*t;
475modelMatrix[46][44]:=t;
476modelMatrix[46][45]:=kappa*t;
477modelMatrix[46][47]:=kappa*t;
478modelMatrix[46][59]:=kappa*c*t;
479modelMatrix[47][15]:=c*t;
480modelMatrix[47][31]:=kappa*c*t;
481modelMatrix[47][35]:=kappa*c*t;
482modelMatrix[47][39]:=c*t;
483modelMatrix[47][43]:=kappa*c*t;
484modelMatrix[47][44]:=kappa*t;
485modelMatrix[47][45]:=t;
486modelMatrix[47][46]:=kappa*t;
487modelMatrix[47][60]:=kappa*c*t;
488modelMatrix[48][1]:=kappa*c*t;
489modelMatrix[48][17]:=c*t;
490modelMatrix[48][33]:=kappa*c*t;
491modelMatrix[48][49]:=t;
492modelMatrix[48][51]:=kappa*c*t;
493modelMatrix[48][54]:=c*t;
494modelMatrix[48][58]:=kappa*c*t;
495modelMatrix[49][3]:=kappa*c*t;
496modelMatrix[49][19]:=c*t;
497modelMatrix[49][35]:=kappa*c*t;
498modelMatrix[49][48]:=t;
499modelMatrix[49][53]:=kappa*c*t;
500modelMatrix[49][56]:=c*t;
501modelMatrix[49][60]:=kappa*c*t;
502modelMatrix[50][4]:=kappa*c*t;
503modelMatrix[50][20]:=c*t;
504modelMatrix[50][36]:=kappa*c*t;
505modelMatrix[50][51]:=kappa*t;
506modelMatrix[50][52]:=t;
507modelMatrix[50][53]:=kappa*t;
508modelMatrix[50][57]:=c*t;
509modelMatrix[51][5]:=kappa*c*t;
510modelMatrix[51][21]:=c*t;
511modelMatrix[51][37]:=kappa*c*t;
512modelMatrix[51][48]:=kappa*c*t;
513modelMatrix[51][50]:=kappa*t;
514modelMatrix[51][52]:=kappa*t;
515modelMatrix[51][53]:=t;
516modelMatrix[51][54]:=kappa*c*t;
517modelMatrix[51][58]:=c*t;
518modelMatrix[52][6]:=kappa*c*t;
519modelMatrix[52][22]:=c*t;
520modelMatrix[52][38]:=kappa*c*t;
521modelMatrix[52][50]:=t;
522modelMatrix[52][51]:=kappa*t;
523modelMatrix[52][53]:=kappa*t;
524modelMatrix[52][55]:=kappa*c*t;
525modelMatrix[52][59]:=c*t;
526modelMatrix[53][7]:=kappa*c*t;
527modelMatrix[53][23]:=c*t;
528modelMatrix[53][39]:=kappa*c*t;
529modelMatrix[53][49]:=kappa*c*t;
530modelMatrix[53][50]:=kappa*t;
531modelMatrix[53][51]:=t;
532modelMatrix[53][52]:=kappa*t;
533modelMatrix[53][56]:=kappa*c*t;
534modelMatrix[53][60]:=c*t;
535modelMatrix[54][9]:=kappa*c*t;
536modelMatrix[54][25]:=c*t;
537modelMatrix[54][41]:=kappa*c*t;
538modelMatrix[54][48]:=c*t;
539modelMatrix[54][51]:=kappa*c*t;
540modelMatrix[54][55]:=kappa*c*t;
541modelMatrix[54][56]:=t;
542modelMatrix[54][58]:=kappa*c*t;
543modelMatrix[55][10]:=kappa*c*t;
544modelMatrix[55][26]:=c*t;
545modelMatrix[55][42]:=kappa*c*t;
546modelMatrix[55][52]:=kappa*c*t;
547modelMatrix[55][54]:=kappa*c*t;
548modelMatrix[55][56]:=kappa*c*t;
549modelMatrix[55][59]:=kappa*c*t;
550modelMatrix[56][11]:=kappa*c*t;
551modelMatrix[56][27]:=c*t;
552modelMatrix[56][43]:=kappa*c*t;
553modelMatrix[56][49]:=c*t;
554modelMatrix[56][53]:=kappa*c*t;
555modelMatrix[56][54]:=t;
556modelMatrix[56][55]:=kappa*c*t;
557modelMatrix[56][60]:=kappa*c*t;
558modelMatrix[57][12]:=kappa*c*t;
559modelMatrix[57][28]:=t;
560modelMatrix[57][44]:=kappa*c*t;
561modelMatrix[57][50]:=c*t;
562modelMatrix[57][58]:=kappa*c*t;
563modelMatrix[57][59]:=t;
564modelMatrix[57][60]:=kappa*c*t;
565modelMatrix[58][13]:=kappa*c*t;
566modelMatrix[58][29]:=c*t;
567modelMatrix[58][45]:=kappa*c*t;
568modelMatrix[58][48]:=kappa*c*t;
569modelMatrix[58][51]:=c*t;
570modelMatrix[58][54]:=kappa*c*t;
571modelMatrix[58][57]:=kappa*c*t;
572modelMatrix[58][59]:=kappa*c*t;
573modelMatrix[58][60]:=t;
574modelMatrix[59][14]:=kappa*c*t;
575modelMatrix[59][30]:=t;
576modelMatrix[59][46]:=kappa*c*t;
577modelMatrix[59][52]:=c*t;
578modelMatrix[59][55]:=kappa*c*t;
579modelMatrix[59][57]:=t;
580modelMatrix[59][58]:=kappa*c*t;
581modelMatrix[59][60]:=kappa*c*t;
582modelMatrix[60][15]:=kappa*c*t;
583modelMatrix[60][31]:=c*t;
584modelMatrix[60][47]:=kappa*c*t;
585modelMatrix[60][49]:=kappa*c*t;
586modelMatrix[60][53]:=c*t;
587modelMatrix[60][56]:=kappa*c*t;
588modelMatrix[60][57]:=kappa*c*t;
589modelMatrix[60][58]:=t;
590modelMatrix[60][59]:=kappa*c*t;
591
592
593vectorOfFrequencies={
594{0.09810421316102973}
595{0.02278249499948084}
596{0.0103838582650695}
597{0.05207427428452762}
598{0.05450234064501652}
599{0.01265694166637824}
600{0.005768810147260831}
601{0.02893015238029313}
602{0.05009576842265348}
603{0.01163361446782}
604{0.005302395709822721}
605{0.02659111878358857}
606{0.0716647798268515}
607{0.01664253180813139}
608{0.007585371640440838}
609{0.03804007270430033}
610{0.02422953796550475}
611{0.005626764740804422}
612{0.002564579847849634}
613{0.01286117655041011}
614{0.01346085442528042}
615{0.003125980411558012}
616{0.001424766582138685}
617{0.00714509808356117}
618{0.0123725300249386}
619{0.00287324156977247}
620{0.001309572688263643}
621{0.006567409302337075}
622{0.01769959156345383}
623{0.004110331690091173}
624{0.001873416484599378}
625{0.009395043863065538}
626{0.04862851325943961}
627{0.01129287748678929}
628{0.005147093820509405}
629{0.02581229139837552}
630{0.02701584069968867}
631{0.006273820825994051}
632{0.002859496566949669}
633{0.0143401618879864}
634{0.02483158123886278}
635{0.005766575737934958}
636{0.002628303227494164}
637{0.01318074454385133}
638{0.03552295649448425}
639{0.00824940695843473}
640{0.003759933783776374}
641{0.0188557873335651}
642{0.006846552901398386}
643{0.01564926377462488}
644{0.01637894174824331}
645{0.003803640500776881}
646{0.001733632064979939}
647{0.008694035430347156}
648{0.003496112119863006}
649{0.00159346606823688}
650{0.007991113416829728}
651{0.02153656595832843}
652{0.005001382615915133}
653{0.002279541736505537}
654{0.0114317316935203}
655};
656Model theModel=(modelMatrix,vectorOfFrequencies);
657TRY_NUMERIC_SEQUENCE_MATCH=0;
658ACCEPT_ROOTED_TREES=0;
659
660UseModel (theModel);
661Tree givenTree=(U68496,U68497,((U68498,(U68499,U68500)Node7)Node5,(U68501,(U68502,(U68503,((((U68504,U68506)Node19,U68505)Node18,U68507)Node17,U68508)Node16)Node14)Node12)Node10)Node4);
662
663givenTree.U68496.t=0.1820765235473329;
664givenTree.U68497.t=0.6321073455214453;
665givenTree.U68498.t=0.3343564501791256;
666givenTree.U68499.t=0.1876640197504523;
667givenTree.U68500.t=0.8005801879569905;
668givenTree.Node7.t=0.750108304945318;
669givenTree.Node5.t=0.189875002788432;
670givenTree.U68501.t=1.055571283712865;
671givenTree.U68502.t=1.747134065548263;
672givenTree.U68503.t=1.779098089271837;
673givenTree.U68504.t=0.1466122835681893;
674givenTree.U68506.t=0.5212546540008767;
675givenTree.Node19.t=0.4171959202637034;
676givenTree.U68505.t=0.2414131900046768;
677givenTree.Node18.t=0.2762126262689943;
678givenTree.U68507.t=0.4784362608407944;
679givenTree.Node17.t=0.006818214720109361;
680givenTree.U68508.t=1.945808887957573;
681givenTree.Node16.t=0.5086054961633685;
682givenTree.Node14.t=0.508023542650578;
683givenTree.Node12.t=0.1816851776949857;
684givenTree.Node10.t=0.394526379685196;
685givenTree.Node4.t=1.247106653589127;SetParameter (DEFER_CONSTRAINT_APPLICATION, 1, 0);
686SetParameter (DEFER_CONSTRAINT_APPLICATION, 0, 0);
687
688DataSet ds = ReadDataFile(USE_NEXUS_FILE_DATA);
689DataSetFilter filteredData = CreateFilter(ds,3,"0-272","0-8,10,9,11,12","TAA,TAG,TGA");
690ASSUME_REVERSIBLE_MODELS=0;
691USE_LAST_RESULTS=1;
692LikelihoodFunction lf = (filteredData,givenTree);
693
694END;