Supplemental results for BASEML

seqf:  trio2.phy
treef: trio2.phy_phyml_tree.txt


TREE #  1

Ancestral reconstruction by BASEML.

(((((((((((((((((H010: 0.031418, H033: 0.031418): 0.045635, (H027: 0.050987, (H017: 0.050987, H014: 0.050987): 0.000000): 0.026065): 0.000000, H011: 0.077052): 0.000000, H003: 0.077053): 0.023554, H018: 0.100607): 0.000000, (((H022: 0.031920, (H013: 0.014845, H019: 0.014845): 0.017075): 0.019267, H034: 0.051187): 0.020614, H029: 0.071801): 0.028806): 0.013321, H020: 0.113928): 0.000000, H001: 0.113928): 0.000000, H015: 0.113928): 0.015048, H028: 0.128977): 0.044636, H030: 0.173613): 0.022224, ((H006: 0.195836, ((((((H031: 0.064366, H005: 0.064366): 0.000000, (H032: 0.027219, H016: 0.027219): 0.037147): 0.000000, H008: 0.064366): 0.000000, H024: 0.064366): 0.028012, H004: 0.092378): 0.038256, H025: 0.130634): 0.065201): 0.000000, OUTG: 0.195836): 0.000000): 0.024127, H035: 0.219963): 0.000000, H007: 0.219963): 0.032073, (H026: 0.073275, H023: 0.073275): 0.178761): 0.043019, H012: 0.295055): 0.000000, (H021: 0.117306, H002: 0.117306): 0.177750, H009: 0.295056);

(((((((((((((((((13, 19), (6, (16, 34))), 32), 22), 35), (((33, (24, 1)), 5), 21)), 4), 12), 10), 14), 28), ((11, ((((((29, 31), (25, 8)), 20), 15), 3), 26)), 36)), 23), 18), (2, 7)), 30), (27, 9), 17);

  37..38   38..39   39..40   40..41   41..42   42..43   43..44   44..45   45..46   46..47   47..48   48..49   49..50   50..51   51..52   52..53   53..13   53..19   52..54   54..6    54..55   55..16   55..34   51..32   50..22   49..35   48..56   56..57   57..58   58..33   58..59   59..24   59..1    57..5    56..21   47..4    46..12   45..10   44..14   43..28   42..60   60..61   61..11   61..62   62..63   63..64   64..65   65..66   66..67   67..29   67..31   66..68   68..25   68..8    65..20   64..15   63..3    62..26   60..36   41..23   40..18   39..69   69..2    69..7    38..30   37..70   70..27   70..9    37..17 

tree with node labels for Rod Page's TreeView
(((((((((((((((((13_H010: 0.031418, 19_H033: 0.031418) 53 : 0.045635, (6_H027: 0.050987, (16_H017: 0.050987, 34_H014: 0.050987) 55 : 0.000000) 54 : 0.026065) 52 : 0.000000, 32_H011: 0.077052) 51 : 0.000000, 22_H003: 0.077053) 50 : 0.023554, 35_H018: 0.100607) 49 : 0.000000, (((33_H022: 0.031920, (24_H013: 0.014845, 1_H019: 0.014845) 59 : 0.017075) 58 : 0.019267, 5_H034: 0.051187) 57 : 0.020614, 21_H029: 0.071801) 56 : 0.028806) 48 : 0.013321, 4_H020: 0.113928) 47 : 0.000000, 12_H001: 0.113928) 46 : 0.000000, 10_H015: 0.113928) 45 : 0.015048, 14_H028: 0.128977) 44 : 0.044636, 28_H030: 0.173613) 43 : 0.022224, ((11_H006: 0.195836, ((((((29_H031: 0.064366, 31_H005: 0.064366) 67 : 0.000000, (25_H032: 0.027219, 8_H016: 0.027219) 68 : 0.037147) 66 : 0.000000, 20_H008: 0.064366) 65 : 0.000000, 15_H024: 0.064366) 64 : 0.028012, 3_H004: 0.092378) 63 : 0.038256, 26_H025: 0.130634) 62 : 0.065201) 61 : 0.000000, 36_OUTG: 0.195836) 60 : 0.000000) 42 : 0.024127, 23_H035: 0.219963) 41 : 0.000000, 18_H007: 0.219963) 40 : 0.032073, (2_H026: 0.073275, 7_H023: 0.073275) 69 : 0.178761) 39 : 0.043019, 30_H012: 0.295055) 38 : 0.000000, (27_H021: 0.117306, 9_H002: 0.117306) 70 : 0.177750, 17_H009: 0.295056) 37 ;

Nodes 37 to 70 are ancestral

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


Prob of best character at each node, listed by site 

Site   Freq   Data: 

   1      1   GGGCGCCGCCCCCCGCCCGGGCCGGGGCGCGCGGCC: C(0.994) C(0.994) C(0.996) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(0.999) C(0.999) C(1.000) C(1.000) C(1.000) C(0.957) C(0.999) C(0.999) G(0.948) G(0.998) G(1.000) G(1.000) C(1.000) C(1.000) G(0.924) G(0.995) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) C(0.848) C(0.847) 
   2      1   CCCGCCCCCGGGCGCCCCCCCCCCCCCCCCCCCCCC: C(0.999) C(0.999) C(1.000) C(1.000) C(1.000) C(0.995) C(0.938) G(0.745) G(0.774) G(0.774) G(0.774) C(0.996) C(0.996) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(0.995) C(0.995) C(0.999) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(0.999) C(0.998) 
   3      1   GCGGGCCGCGGGGGGCCGGGGGGGGGCGGCGGGCGG: C(0.937) C(0.937) C(0.760) G(0.973) G(0.973) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) C(0.995) C(0.995) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) C(0.992) C(0.993) 
   4      1   GCGGGCCCCGGGCGGCCCCGGCCGCGCGGCGCGCGG: C(0.994) C(0.994) C(0.981) C(0.926) C(0.926) G(0.997) G(0.999) G(1.000) G(1.000) G(1.000) G(1.000) G(0.998) G(0.998) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(0.997) G(0.997) G(0.999) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) C(0.981) C(0.999) C(0.997) 
   5      1   CGCGGGGCCGGGGGCGGGGCGCGCCGCGCGCGCGGG: G(0.979) G(0.979) G(0.997) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(0.995) G(0.944) C(0.974) C(1.000) G(1.000) G(1.000) G(0.826) C(0.948) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) G(0.999) C(0.930) 
   6      1   GCCCCGCGGCCCCCGCGCCGCCCGGCGCGGGCCCCC: G(0.755) G(0.755) C(0.908) C(0.997) C(0.997) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(0.999) C(0.999) C(1.000) C(0.998) C(0.964) G(0.988) C(1.000) C(1.000) C(0.976) C(0.869) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) C(0.996) G(0.980) 
   7      1   CGGCCCGGGCGCCCGCGGCGCCGCCGGGGGGCCCCG: G(0.999) G(0.999) G(1.000) G(1.000) G(1.000) G(1.000) G(0.923) C(0.962) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(0.964) G(0.999) G(0.998) 
   8      1   CGGCCCGGGCGCCGGCCGCGCCGCGGCGCGGCCCCG: G(0.952) G(0.952) G(0.994) G(1.000) G(1.000) G(1.000) G(0.976) G(0.705) C(0.999) C(0.999) C(0.999) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(0.999) G(0.816) 
   9      1   CGGCCCGGGGGCCCCCGGCGCCGCGGGCGGGCCCCG: G(0.999) G(0.999) G(1.000) G(1.000) G(1.000) G(0.998) C(0.737) C(0.988) C(0.999) C(0.999) C(0.999) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) G(0.998) G(0.998) G(0.999) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(0.999) G(0.998) 
  10      1   CCCGGGCCCCCCGCCCCCGCCCGGCCCCCCGCGCCC: C(0.999) C(0.999) C(0.999) C(0.999) C(0.999) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) G(0.979) C(0.999) C(0.999) C(0.920) G(0.894) G(0.932) G(0.918) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(1.000) C(0.999) C(0.998) 
  11      1   GGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGGC: G(0.999) G(0.999) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(1.000) G(0.999) G(0.998) 

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

Branch 1:   37..38



Branch 2:   38..39

	   6 G 0.755 -> C 0.908


Branch 3:   39..40

	   3 C 0.760 -> G 0.973


Branch 4:   40..41



Branch 5:   41..42

	   4 C 0.926 -> G 0.997


Branch 6:   42..43

	   9 G 0.998 -> C 0.737


Branch 7:   43..44

	   2 C 0.938 -> G 0.745
	   7 G 0.923 -> C 0.962


Branch 8:   44..45

	   8 G 0.705 -> C 0.999


Branch 9:   45..46



Branch 10:   46..47



Branch 11:   47..48

	   2 G 0.774 -> C 0.996


Branch 12:   48..49



Branch 13:   49..50

	   4 G 0.998 -> C 1.000


Branch 14:   50..51



Branch 15:   51..52



Branch 16:   52..53

	  10 C 1.000 -> G 0.979


Branch 17:   53..13 (H010) 



Branch 18:   53..19 (H033) 

	   1 C 0.957 -> G


Branch 19:   52..54

	   3 G 1.000 -> C 0.995


Branch 20:   54..6  (H027) 

	   6 C 0.999 -> G
	  10 C 0.999 -> G


Branch 21:   54..55



Branch 22:   55..16 (H017) 



Branch 23:   55..34 (H014) 

	   1 C 0.999 -> G


Branch 24:   51..32 (H011) 



Branch 25:   50..22 (H003) 

	   5 G 1.000 -> C


Branch 26:   49..35 (H018) 



Branch 27:   48..56

	   1 C 0.999 -> G 0.948


Branch 28:   56..57

	  10 C 0.920 -> G 0.894


Branch 29:   57..58

	   5 G 0.944 -> C 0.974


Branch 30:   58..33 (H022) 



Branch 31:   58..59

	   6 C 0.964 -> G 0.988


Branch 32:   59..24 (H013) 



Branch 33:   59..1  (H019) 

	  10 G 0.918 -> C


Branch 34:   57..5  (H034) 



Branch 35:   56..21 (H029) 



Branch 36:   47..4  (H020) 

	  10 C 1.000 -> G


Branch 37:   46..12 (H001) 



Branch 38:   45..10 (H015) 

	   9 C 0.999 -> G


Branch 39:   44..14 (H028) 



Branch 40:   43..28 (H030) 



Branch 41:   42..60



Branch 42:   60..61



Branch 43:   61..11 (H006) 

	   2 C 0.995 -> G


Branch 44:   61..62

	   1 C 1.000 -> G 0.924


Branch 45:   62..63

	   5 G 0.826 -> C 0.948


Branch 46:   63..64

	   6 C 0.869 -> G 1.000


Branch 47:   64..65



Branch 48:   65..66



Branch 49:   66..67



Branch 50:   67..29 (H031) 

	   8 G 1.000 -> C


Branch 51:   67..31 (H005) 

	  10 C 1.000 -> G


Branch 52:   66..68

	   4 G 1.000 -> C 0.981


Branch 53:   68..25 (H032) 

	   7 G 0.964 -> C


Branch 54:   68..8  (H016) 



Branch 55:   65..20 (H008) 



Branch 56:   64..15 (H024) 

	   9 G 1.000 -> C


Branch 57:   63..3  (H004) 



Branch 58:   62..26 (H025) 



Branch 59:   60..36 (OUTG) 

	  11 G 1.000 -> C


Branch 60:   41..23 (H035) 

	  10 C 0.999 -> G


Branch 61:   40..18 (H007) 



Branch 62:   39..69



Branch 63:   69..2  (H026) 

	   1 C 0.848 -> G


Branch 64:   69..7  (H023) 



Branch 65:   38..30 (H012) 



Branch 66:   37..70

	   5 G 0.979 -> C 0.930


Branch 67:   70..27 (H021) 

	   1 C 0.847 -> G
	   8 G 0.816 -> C


Branch 68:   70..9  (H002) 



Branch 69:   37..17 (H009) 

	   8 G 0.952 -> C




List of extant and reconstructed sequences

H019              GCGGCGCCCC G
H026              GCCCGCGGGC G
H004              GCGGCCGGGC G
H020              CGGGGCCCCG G
H034              GCGGGCCCCG G
H027              CCCCGGCCCG G
H023              CCCCGCGGGC G
H016              GCGCCGGGGC G
H002              CCCCCGGGGC G
H015              CGGGGCCCGC G
H006              CGGGGCGGGC G
H001              CGGGGCCCCC G
H010              CCGCGCCCCG G
H028              CGGGGCCGCC G
H024              GCGGCGGGCC G
H017              CCCCGCCCCC G
H009              CCCCGGGCGC G
H007              CCGCGCGGGC G
H033              GCGCGCCCCG G
H008              GCGGCGGGGC G
H029              GCGGGCCCCC G
H003              CCGCCCCCCC G
H035              CCGCGCGGGG G
H013              GCGGCGCCCG G
H032              GCGCCGCGGC G
H025              GCGGGCGGGC G
H021              GCCCCGGCGC G
H030              CCGGGCGGCC G
H031              GCGGCGGCGC G
H012              CCCCGGGGGC G
H005              GCGGCGGGGG G
H011              CCGCGCCCCC G
H022              GCGGCCCCCG G
H014              GCCCGCCCCC G
H018              CCGGGCCCCC G
OUTG              CCGGGCGGGC C
node #37          CCCCGGGGGC G
node #38          CCCCGGGGGC G
node #39          CCCCGCGGGC G
node #40          CCGCGCGGGC G
node #41          CCGCGCGGGC G
node #42          CCGGGCGGGC G
node #43          CCGGGCGGCC G
node #44          CGGGGCCGCC G
node #45          CGGGGCCCCC G
node #46          CGGGGCCCCC G
node #47          CGGGGCCCCC G
node #48          CCGGGCCCCC G
node #49          CCGGGCCCCC G
node #50          CCGCGCCCCC G
node #51          CCGCGCCCCC G
node #52          CCGCGCCCCC G
node #53          CCGCGCCCCG G
node #54          CCCCGCCCCC G
node #55          CCCCGCCCCC G
node #56          GCGGGCCCCC G
node #57          GCGGGCCCCG G
node #58          GCGGCCCCCG G
node #59          GCGGCGCCCG G
node #60          CCGGGCGGGC G
node #61          CCGGGCGGGC G
node #62          GCGGGCGGGC G
node #63          GCGGCCGGGC G
node #64          GCGGCGGGGC G
node #65          GCGGCGGGGC G
node #66          GCGGCGGGGC G
node #67          GCGGCGGGGC G
node #68          GCGCCGGGGC G
node #69          CCCCGCGGGC G
node #70          CCCCCGGGGC G


Overall accuracy of the 34 ancestral sequences:
  0.96430  0.96430  0.96673  0.99030  0.99030  0.99897  0.96110  0.94544  0.97923  0.97923  0.97923  0.99936  0.99936  0.99992  0.99992  0.99992  0.99410  0.99918  0.99918  0.98743  0.98479  0.98822  0.99134  0.99897  0.99897  0.97480  0.98281  1.00000  1.00000  1.00000  1.00000  0.99492  0.98441  0.95925
for a site.

  0.64910  0.64910  0.66689  0.89584  0.89584  0.98872  0.62168  0.49918  0.77199  0.77199  0.77199  0.99302  0.99302  0.99913  0.99913  0.99913  0.93600  0.99103  0.99103  0.86669  0.83945  0.87553  0.90582  0.98872  0.98872  0.74263  0.81876  0.99995  0.99995  0.99995  0.99995  0.94489  0.83161  0.61662
for the sequence.


Counts of changes at sites

   1  CG CG CG CG CG CG (6)
   2  CG CG GC (3)
   3  CG GC (2)
   4  CG GC GC (3)
   5  GC GC GC GC (4)
   6  CG GC CG CG (4)
   7  GC GC (2)
   8  GC GC GC GC (4)
   9  CG GC GC (3)
  10  GC CG CG CG CG CG CG (7)
  11  GC (1)


