
seed used = 140031181

K12                                  ATG AGA CGA ACA TTT ACA GCA GAG GAA AAA GCC TCT GTT TTT GAA CTA TGG AAG AAC GGA ACA GGC TTC AGT GAA ATA GCG AAT ATC CTG GGT TCA AAA CCC GGA ACG ATC TTC ACT ATG TTA AGG GAT ACT GGC GGC ATA AAA CCC CAT GAG CGT AAG CGG GCT GTA GCT CAC CTG ACA CTG TCT GAG CGC GAG GAG ATA CGA GCT GGT TTG TCA GCC AAA ATG AGC ATT CGT GCG ATA GCT ACT GCG CTG AAT CGC AGT CCT TCG ACG ATC TCA CGT GAA GTT CAG CGT AAT CGG GGC AGA CGC TAT TAC AAA GCT GTT GAT GCT AAT AAC CGA GCC AAC AGA ATG GCG AAA AGG CCA AAA CCG TGC TTA CTG GAT CAA AAT TTA CCA TTG CGA AAG CTT GTT CTG GAA AAG CTG GAG ATG AAA TGG TCT CCA GAG CAA ATA TCA GGA TGG TTA AGG CGA ACA AAA CCA CGT CAA AAA ACG CTG CGA ATA TCA CCT GAG ACA ATT TAT AAA ACG CTG TAC TTT CGT AGC CGT GAA GCG CTA CAC CAC CTG AAT ATA CAG CAT CTG CGA CGG TCG CAT AGC CTT CGC CAT GGC AGG CGT CAT ACC CGC AAA GGC GAA AGA GGT ACG ATT AAC ATA GTG AAC GGA ACA CCA ATT CAC GAA CGT TCC CGA AAT ATC GAT AAC AGA CGC TCT CTA GGG CAT TGG GAG GGC GAT TTA GTC TCA GGT ACA AAA AAC TCT CAT ATA GCC ACA CTT GTA GAC CGA AAA TCA CGT TAT ACG ATC ATC CTT AGA CTC AGG GGC AAA GAT TCT GTC TCA GTA AAT CAG GCT CTT ACC GAC AAA TTC CTG AGT TTA CCG TCA GAA CTC AGA AAA TCA CTG ACA TGG GAC AGA GGA ATG GAA CTG GCC AGA CAT CTA GAA TTT ACT GTC AGC ACC GGC GTT AAA GTT TAC TTC TGC GAT CCT CAG AGT CCT TGG CAG CGG GGA ACA AAT GAG AAC ACA AAT GGG CTA ATT CGG CAG TAC TTT CCT AAA AAG ACA TGT CTT GCC CAA TAT ACT CAA CAT GAA CTA GAT CTG GTT GCT GCT CAG CTA AAC AAC AGA CCG AGA AAG ACA CTG AAG TTC AAA ACA CCG AAA GAG ATA ATT GAA AGG GGT GTT GCA TTG ACA GAT 
CFT073                               ATG AGA CGA ACA TTT ACA GCA GAG GAA AAA GCC TCT GTT TTT GAA CTA TGG AAG AAC GGA ACA GGC TTC AGT GAA ATA GCG AAT ATC CTG GGT TCA AAA CCC GGA ACG ATC TTC ACT ATG TTA AGG GAT ACT GGC GGC ATA AAA CCC CAT GAG CGT AAG CGG GCT GTA GCT CAC CTG ACA CTG TCT GAG CGC GAG GAG ATA CGA GCT GGT TTG TCA GCC AAA ATG AGC ATT CGT GCG ATA GCT ACT GCG CTG AAT CGC AGT CCT TCG ACG ATC TCA CGT GAA GTT CAG CGT AAT CGG GGC AGA CGC TAT TAC AAA GCT GTT GAT GCT AAT AAC CGA GCC AAC AGA ATG GCG AAA AGG CCA AAA CCG TGC TTA CTG GAT CAA AAT TTG CCA TTG CGA AAG CTT GTT CTG GAA AAG CTG GAG ATG AAA TGG TCT TCA GAG CAA ATA TCA GGA TGG TTA AGG CGA ACA AAA CCA CGT CAA AAA ACG CTG CGA ATA TCA CCT GAG ACA ATT TAT AAA ACG CTG TAC TTT CGT AGC CGT GAA GCG CTA CAC CAC CTG AAT ATA CAG CAT CTG CGA CGG TCG CAT AGC CTT CGC TAT GGC AGG CGT CAT ACC CGC AAA GGC GAA AGA GGT ACG ATT AAC ATA GTG AAC GGA ACA CCA ATT CAC GAA CGT TCC CGA AAT ATC GAT AAC AGA CGC TCT CTG GGG CAT TGG GAG GGC GAT TTA GTC TCA GGT ACA AAA AAC TCT CAT ATA GCC ACA CTT GTA GAC CGA AAA TCA CGT TAT ACG ATC ATC CTT AGA CTC AGG GGC AAA GAT TCT GTC TCA GTA AAT CAG GCT CTT ACC GAC AAA TTC CTG AGT TTA CCG TCA GAA CTC AGA AAA TCA CTG ACA TGG GAC AGA GGA ATG GAA CTG GCC AGA CAT CTA GAA TTT ACT GTC AGC ACC GGC GTT AAA GTT TAC TTC TGC GAT CCT CAG AGT CCT TGG CAG CGG GGA ACA AAT GAG AAC ACA AAT GGG CTA ATT CGG CAG TAC TTT CCT AAA AAG ACA TGT CTT GCC CAA TAT ACT CAA CAT GAA CTA GAT CTG GTT GCT GCT CAG CTA AAC AAC AGA CCG AGA AAG ACA CTG AAG TTC AAA ACA CCG AAA GAG ATA ATT GAA AGG GGT GTT GCA TTG ACA GAT 
EDL933                               CAA CGC AAC AGT TAT GTG AAA ACA TGG GGT TGC GGA GGT TTG GAA CTA TGG AAG AAC GGA ACA GGC TTC AGT GAA ATA GCG AAT ATC CTG GGT TCA AAA CCC GGA ACG ATC TTC ACT ATG TTA AGG GAT ACT GGC GGC ATA AAA CCC CAT GAG CGT AAG CGG GCT GTA GCT CAC CTG ACA CTG TCT GAG CGC GAG GAG ATA CGA GCT GGT TTG TCA GCC AAA ATG AGC ATT CGT GCG ATA GCT ACT GCG CTG AAT CGC AGT CCT TCG ACG ATC TCA CGT GAA GTT CAG CGT AAT CGG GGC AGA CGC TAT TAC AAA GCT GTT GAT GCT AAT AAC CGA GCC AAC AGA ATG GCG AAA AGG CCA AAA CCG TGC TTA CTG GAT CAA AAT TTA CCA TTG CGA AAG CTT GTT CTG GAA AAG CTG GAG ATG AAA TGG TCT CCA GAG CAA ATA TCA GGA TGG TTA AGG CGA ACA AAA CCA CGT CAA AAA ACG CTG CGA ATA TCA CCT GAG ACA ATT TAT AAA ACG CTG TAC TTT CGT AGC CGT GAA GCG CTA CAC CAC CTG AAT ATA CAG CAT CTG CGA CGG TCG CAT AGC CTT CGC CAT GGC AGG CGT CAT ACC CGC AAA GGC GAA AGA GGT ACG ATT AAC ATA GTG AAC GGA ACA CCA ATT CAC GAA CGT TCC CGA AAT ATC GAT AAC AGA CGC TCT CTA GGG CAT TGG GAG GGC GAT TTA GTC TCA GGT ACA AAA AAC TCT CAT ATA GCC ACA CTT GTA GAC CGA AAA TCA CGT TAT ACG ATC ATC CTC AGA CTC AGG GGC AAA GAT TCT GTC TCA GTA AAT CAG GCT CTT ACC GAC AAA TTC CTG AGT TTA CCG TCA GAA CTC AGA AAA TCA CTG ACA TGG GAC AGA GGA ATG GAA CTG GCC AGA CAT CTA GAA TTT ACT GTC AGC ACC GGC GTT AAA GTT TAC TTC TGC GAT CCT CAG AGT CCT TGG CAG CGG GGA ACA AAT GAG AAC ACA AAT GGG CTA ATT CGG CAG TAC TTT CCT AAA AAG ACA TGT CTT GCC CAA TAT ACT CAA CAT GAA --- GAT CTG GTT GCT GCT CAG CTA AAC AAC AGA CCG AGA AAG ACA CTG AAG TTC AAA ACA CCG AAA GAG ATA ATT GAA AGG GGT GTT GCA TTG ACA GAT 


CODONML (in paml 3.14, January 2004)    H:\users\mt269\runpaml\NEWdna80\b4284.phy   Model: One dN/dS ratio 
Codon frequencies: F3x4
Site-class models:  PositiveSelection

ns = 3  	ls = 383
# site patterns = 81
    1    1    1    1    1    1    1    1    1    1    1    1    1    1   12
    5    6    7   10    6   14    9    5    4   10    5   10    6   16    5
   10   20    2    6    5    5    5    6    9    8   10   10    5    9    3
    4    5    6    8    3    5    4    6    5    2    7    7   10    4    4
    4    4    2    5    1    5    1    3    1    3    1    1    1    2    3
    3    1    2    1    1    1

K12                                  ATG AGA CGA ACA TTT ACA GCA GAG GAA AAA GCC TCT GTT TTT GAA CTA TGG AAG AAC GGA ACA GGC TTC AGT ATA GCG AAT ATC CTG GGT TCA AAA CCC ACG ACT ATG TTA AGG GAT CAT GAG CGT CGG GCT GTA CAC TCT CGC CGA TTG GCC AGC ATT CCT TCG GTT CAG AGA TAT TAC CCA CCG TGC CAA TTA CTT CCA TTT CAT ACC GTG TCC CTA GGG GTC GAC CTT CTC TGT CTA GCA 
CFT073                               ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..G ... T.. ... T.. ... ... ... ..G ... ... ... ... ... ... ... ... 
EDL933                               CAA C.C AAC .GT .A. GTG AA. ACA TGG GGT TG. GGA .G. ..G ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..C ... ... --- ... 



Codon usage in sequences
--------------------------------------------------------------------------
Phe TTT  5  5  3 | Ser TCT  6  6  5 | Tyr TAT  4  5  5 | Cys TGT  1  1  1
    TTC  5  5  5 |     TCC  1  1  1 |     TAC  4  4  4 |     TGC  2  2  3
Leu TTA  6  5  6 |     TCA 10 11 10 | *** TAA  0  0  0 | *** TGA  0  0  0
    TTG  3  4  4 |     TCG  2  2  2 |     TAG  0  0  0 | Trp TGG  6  6  7
--------------------------------------------------------------------------
Leu CTT  6  6  5 | Pro CCT  5  5  5 | His CAT  9  8  9 | Arg CGT 10 10 10
    CTC  2  2  3 |     CCC  2  2  2 |     CAC  4  4  4 |     CGC  6  6  7
    CTA  7  6  6 |     CCA  5  4  5 | Gln CAA  5  5  6 |     CGA  9  9  8
    CTG 16 17 16 |     CCG  4  4  4 |     CAG  7  7  7 |     CGG  5  5  5
--------------------------------------------------------------------------
Ile ATT  6  6  6 | Thr ACT  5  5  5 | Asn AAT 10 10 10 | Ser AGT  4  4  5
    ATC  6  6  6 |     ACC  3  3  3 |     AAC 10 10 11 |     AGC  4  4  4
    ATA 10 10 10 |     ACA 16 16 15 | Lys AAA 21 21 21 | Arg AGA 11 11 10
Met ATG  6  6  5 |     ACG  6  6  6 |     AAG  7  7  7 |     AGG  6  6  6
--------------------------------------------------------------------------
Val GTT  8  8  7 | Ala GCT  9  9  9 | Asp GAT  9  9  9 | Gly GGT  5  5  7
    GTC  3  3  3 |     GCC  6  6  5 |     GAC  3  3  3 |     GGC  9  9  9
    GTA  3  3  3 |     GCA  2  2  1 | Glu GAA 13 13 12 |     GGA  6  6  7
    GTG  1  1  2 |     GCG  5  5  5 |     GAG 11 11 10 |     GGG  2  2  2
--------------------------------------------------------------------------

Codon position x base (3x4) table for each sequence.

#1: K12            
position  1:    T:0.14360    C:0.26632    A:0.34204    G:0.24804
position  2:    T:0.24282    C:0.22715    A:0.30548    G:0.22454
position  3:    T:0.26632    C:0.18277    A:0.32376    G:0.22715

#2: CFT073         
position  1:    T:0.14883    C:0.26110    A:0.34204    G:0.24804
position  2:    T:0.24282    C:0.22715    A:0.30548    G:0.22454
position  3:    T:0.26632    C:0.18277    A:0.31854    G:0.23238

#3: EDL933         
position  1:    T:0.14660    C:0.26702    A:0.34031    G:0.24607
position  2:    T:0.23560    C:0.21728    A:0.30890    G:0.23822
position  3:    T:0.26440    C:0.19110    A:0.31414    G:0.23037

Sums of codon usage counts
------------------------------------------------------------------------------
Phe F TTT      13 | Ser S TCT      17 | Tyr Y TAT      14 | Cys C TGT       3
      TTC      15 |       TCC       3 |       TAC      12 |       TGC       7
Leu L TTA      17 |       TCA      31 | *** * TAA       0 | *** * TGA       0
      TTG      11 |       TCG       6 |       TAG       0 | Trp W TGG      19
------------------------------------------------------------------------------
Leu L CTT      17 | Pro P CCT      15 | His H CAT      26 | Arg R CGT      30
      CTC       7 |       CCC       6 |       CAC      12 |       CGC      19
      CTA      19 |       CCA      14 | Gln Q CAA      16 |       CGA      26
      CTG      49 |       CCG      12 |       CAG      21 |       CGG      15
------------------------------------------------------------------------------
Ile I ATT      18 | Thr T ACT      15 | Asn N AAT      30 | Ser S AGT      13
      ATC      18 |       ACC       9 |       AAC      31 |       AGC      12
      ATA      30 |       ACA      47 | Lys K AAA      63 | Arg R AGA      32
Met M ATG      17 |       ACG      18 |       AAG      21 |       AGG      18
------------------------------------------------------------------------------
Val V GTT      23 | Ala A GCT      27 | Asp D GAT      27 | Gly G GGT      17
      GTC       9 |       GCC      17 |       GAC       9 |       GGC      27
      GTA       9 |       GCA       5 | Glu E GAA      38 |       GGA      19
      GTG       4 |       GCG      15 |       GAG      32 |       GGG       6
------------------------------------------------------------------------------

(Ambiguity data are not used in the counts.)


Codon position x base (3x4) table, overall

position  1:    T:0.14634    C:0.26481    A:0.34146    G:0.24739
position  2:    T:0.24042    C:0.22387    A:0.30662    G:0.22909
position  3:    T:0.26568    C:0.18554    A:0.31882    G:0.22997

Codon frequencies under model, for use in evolver:
  0.00968961  0.00676684  0.01162753  0.00838707
  0.00902257  0.00630101  0.01082708  0.00780970
  0.01235776  0.00863017  0.00000000  0.00000000
  0.00923321  0.00644811  0.00000000  0.00799203
  0.01753357  0.01224476  0.02104029  0.01517660
  0.01632655  0.01140182  0.01959186  0.01413184
  0.02236166  0.01561650  0.02683399  0.01935567
  0.01670772  0.01166801  0.02004926  0.01446176
  0.02260908  0.01578929  0.02713090  0.01956983
  0.02105266  0.01470235  0.02526319  0.01822263
  0.02883477  0.02013707  0.03460173  0.02495862
  0.02154416  0.01504559  0.02585299  0.01864806
  0.01638005  0.01143918  0.01965606  0.01417814
  0.01525244  0.01065170  0.01830292  0.01320211
  0.02089050  0.01458910  0.02506860  0.01808227
  0.01560852  0.01090038  0.01873023  0.01351033



Nei & Gojobori 1986. dN/dS (dN, dS)
(Pairwise deletion)
(Note: This matrix is not used in later m.l. analysis.
Use runmode = -2 for ML pairwise comparison.)

K12                 
CFT073               0.3345 (0.0023 0.0070)
EDL933               1.0442 (0.0297 0.0284) 0.8999 (0.0321 0.0357)


TREE #  1:  (1, (2, 3));   MP score: -1
This is a rooted tree.  Please check!
lnL(ntime:  4  np:  9):  -1717.975953     +0.000000
   4..1     4..5     5..2     5..3  
  0.00000  0.00000  0.01046  0.52726  9.17898  0.95602  0.00000  0.00000 653.05141
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.53772

(1: 0.000004, (2: 0.010461, 3: 0.527255): 0.000000);

(K12: 0.000004, (CFT073: 0.010461, EDL933: 0.527255): 0.000000);

Detailed output identifying parameters

kappa (ts/tv) =  9.17898


dN/dS for site classes (K=3)

p:   0.95602  0.00000  0.04398
w:   0.00000  1.00000653.05141

dN & dS for each branch

 branch           t        S        N    dN/dS       dN       dS   S*dS   N*dN

   4..1       0.000    353.6    795.4  28.7198   0.0000   0.0000    0.0    0.0
   4..5       0.000    353.6    795.4  28.7198   0.0000   0.0000    0.0    0.0
   5..2       0.010    353.6    795.4  28.7198   0.0050   0.0002    0.1    3.9
   5..3       0.527    353.6    795.4  28.7198   0.2500   0.0087    3.1  198.9


Positively selected sites (*: P>95%; **: P>99%)

            Pr(w>1)     post mean +- SE for w

     1 M      1.000**       653.051
     2 R      0.999**       652.385
     3 R      1.000**       653.051
     4 T      1.000**       653.051
     5 F      1.000**       653.051
     6 T      1.000**       653.051
     7 A      1.000**       653.051
     8 E      1.000**       653.051
     9 E      1.000**       653.051
    10 K      1.000**       653.051
    11 A      1.000**       653.051
    12 S      1.000**       653.051
    13 V      1.000**       653.051
    14 F      1.000**       653.051
   145 P      1.000**       653.051
   197 H      1.000**       653.051


Bayes Empirical Bayes (BEB) analysis
Positively selected sites (*: P>95%; **: P>99%)

            Pr(w>1)     post mean +- SE for w

     1 M      0.997**       10.167 +- 0.782
     2 R      0.943         9.653 +- 2.290
     3 R      0.997**       10.169 +- 0.773
     4 T      0.974*        9.952 +- 1.600
     5 F      0.853         8.809 +- 3.377
     6 T      0.995**       10.147 +- 0.894
     7 A      0.985*        10.056 +- 1.274
     8 E      0.986*        10.066 +- 1.239
     9 E      0.987*        10.073 +- 1.210
    10 K      0.998**       10.175 +- 0.735
    11 A      0.987*        10.073 +- 1.212
    12 S      0.991**       10.107 +- 1.078
    13 V      0.852         8.805 +- 3.380
    14 F      0.854         8.822 +- 3.365
   145 P      0.846         8.751 +- 3.433
   197 H      0.846         8.744 +- 3.439



The grid (see ternary graph for p0-p1)

w0:   0.050  0.150  0.250  0.350  0.450  0.550  0.650  0.750  0.850  0.950
w2:   1.500  2.500  3.500  4.500  5.500  6.500  7.500  8.500  9.500 10.500


Posterior on the grid

w0:   0.233  0.181  0.141  0.111  0.088  0.071  0.057  0.047  0.039  0.032
w2:   0.000  0.000  0.000  0.000  0.000  0.001  0.008  0.042  0.191  0.757

Posterior for p0-p1 (see the ternary graph)

 0.000
 0.002 0.000 0.000
 0.012 0.006 0.002 0.000 0.000
 0.007 0.013 0.016 0.007 0.003 0.000 0.000
 0.001 0.004 0.009 0.017 0.021 0.009 0.003 0.000 0.000
 0.000 0.001 0.002 0.006 0.012 0.023 0.027 0.011 0.004 0.000 0.000
 0.000 0.000 0.000 0.001 0.002 0.008 0.017 0.031 0.036 0.014 0.005 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.012 0.024 0.043 0.049 0.018 0.006 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.004 0.017 0.034 0.059 0.067 0.024 0.007 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.003 0.007 0.024 0.049 0.083 0.092 0.030 0.009 0.000 0.000

sum of density on p0-p1 =   1.000000
