Supplemental results for CODEML (seqf: Homogentisate_Aln_Tree.nex_132.nex  treef: OptimizedTreeFile.tre)


TREE #  1

Ancestral reconstruction by AAML.

((HGD_HUM: 0.002246, (B8YB76_B: 0.029948, (HGD_MOU: 0.015790, Q6AYR0: 0.027504): 0.035239): 0.022660): 0.000000, 1EY2_A: 0.000000);

((2, (3, (4, 5))), 1);

   6..7     7..2     7..8     8..3     8..9     9..4     9..5     6..1  

tree with node labels for Rod Page's TreeView
((2_HGD_HUM, (3_B8YB76_B, (4_HGD_MOU, 5_Q6AYR0) 9 ) 8 ) 7 , 1_1EY2_A) 6 ;

Nodes 6 to 9 are ancestral

(1) Marginal reconstruction of ancestral sequences
(eqn. 4 in Yang et al. 1995 Genetics 141:1641-1650).

Prob distribs at nodes, those with p < 0.001 not listed

Prob distribution at node 6, by site

   site  Freq   Data

      1      1   IIIVV: A(0.000) R(0.000) N(0.000) D(0.000) C(0.000) Q(0.000) E(0.000) G(0.000) H(0.000) I(1.000) L(0.000) K(0.000) M(0.000) F(0.000) P(0.000) S(0.000) T(0.000) W(0.000) Y(0.000) V(0.000) 

Prob distribution at node 7, by site

   site  Freq   Data

      1      1   IIIVV: A(0.000) R(0.000) N(0.000) D(0.000) C(0.000) Q(0.000) E(0.000) G(0.000) H(0.000) I(1.000) L(0.000) K(0.000) M(0.000) F(0.000) P(0.000) S(0.000) T(0.000) W(0.000) Y(0.000) V(0.000) 

Prob distribution at node 8, by site

   site  Freq   Data

      1      1   IIIVV: A(0.000) R(0.000) N(0.000) D(0.000) C(0.000) Q(0.000) E(0.000) G(0.000) H(0.000) I(0.991) L(0.000) K(0.000) M(0.000) F(0.000) P(0.000) S(0.000) T(0.000) W(0.000) Y(0.000) V(0.009) 

Prob distribution at node 9, by site

   site  Freq   Data

      1      1   IIIVV: A(0.000) R(0.000) N(0.000) D(0.000) C(0.000) Q(0.000) E(0.000) G(0.000) H(0.000) I(0.005) L(0.000) K(0.000) M(0.000) F(0.000) P(0.000) S(0.000) T(0.000) W(0.000) Y(0.000) V(0.995) 

Prob of best state at each node, listed by site

   site   Freq   Data: 

   1      1   IIIVV: I(1.000) I(1.000) I(0.991) V(0.995) 

Summary of changes along branches.
Check root of tree for directions of change.

Branch 1:    6..7 



Branch 2:    7..2  (HGD_HUM) 



Branch 3:    7..8 



Branch 4:    8..3  (B8YB76_B) 



Branch 5:    8..9 

	   1 I 0.991 -> V 0.995


Branch 6:    9..4  (HGD_MOU) 



Branch 7:    9..5  (Q6AYR0) 



Branch 8:    6..1  (1EY2_A) 





List of extant and reconstructed sequences

     9      1

1EY2_A            I
HGD_HUM           I
B8YB76_B          I
HGD_MOU           V
Q6AYR0            V
node #6           I
node #7           I
node #8           I
node #9           V


Overall accuracy of the 4 ancestral sequences:
  1.00000  1.00000  0.99058  0.99458
for a site.

  1.00000  1.00000  0.99058  0.99458
for the sequence.


Counts of changes at sites

   1  IV (1)


(2) Joint reconstruction of ancestral sequences
(eqn. 2 in Yang et al. 1995 Genetics 141:1641-1650), using the algorithm of Pupko et al. (2000 Mol Biol Evol 17:890-896),
modified to generate sub-optimal reconstructions.

Reconstruction (prob.), listed by pattern (use the observed data to find the right site).

Pattern Freq   Data:

   1       1  IIIVV: IIIV (0.985) IIVV (0.009) IIII (0.005) IIMV (0.000) VVVV (0.000) VVIV (0.000)  (total  1.000)



List of extant and reconstructed sequences

     9      1

1EY2_A            I
HGD_HUM           I
B8YB76_B          I
HGD_MOU           V
Q6AYR0            V
node #6           I
node #7           I
node #8           I
node #9           V


