
seed used = 139989957

K12                                  TTG GAT GTG AAC GTT GAT CAG TTC GAT ACT GAA GCT TTC CGT ACT GAC AAA CTG GAA CTG ACC AGC GGC AAC ATC GCT GAC CAT AAC GGT AAC GTA GTA TCT GGT GTG TTC GAT ATC CAT AGC AGC GAT TAC GTT CTG AAC GCT GAT CTG GTG AAC GAC CGT ACC TGG GAT ACT TCC AAG TCT AAC TAC GGT TAC GGT ATT GTT GCT ATG AAC TCT GAT GGT CAC CTG ACT ATC AAC GGT AAC GGC GAC GTA GAC AAC --- --- GGT ACT GAA CTG GAT AAC AGC TCT GTA GAC AAT --- --- --- --- GTT GTT GCT GCA ACC GGT AAC TAC AAA GTT CGT ATC GAC AAC GCA ACT GGC GCT GGC GCT ATC GCT GAT TAC AAA GAT AAA GAA ATT ATC TAC GTA AAC GAC GTC AAC AGC AAC GCG ACC TTC TCT GCT GCT AAC AAA GCT GAC CTG GGT GCA TAC ACC TAT CAG GCT GAA CAG CGC GGT AAC ACC GTT GTT CTG CAA CAG ATG GAG CTG ACC GAC TAC GCT AAC ATG GCG CTG AGC ATC CCG TCT GCG AAC ACC AAT ATC TGG AAC CTG GAA CAA GAC ACC GTT GGT ACT CGT CTG ACC AAC TCT CGT CAT GGC CTG GCT GAT AAC GGC GGC GCA TGG GTA AGC TAC TTC GGT GGT AAC TTC AAC GGC GAC AAC GGC ACC ATC AAC TAT GAT CAG GAT GTT AAC GGC ATC ATG GTC GGT GTT GAT ACC AAA ATT GAC GGT AAC AAC GCT AAG TGG ATC GTC GGT GCG GCT GCA GGC TTC GCT AAA GGT GAC ATG AAT GAC CGT TCT GGT CAG GTG GAT CAA GAC AGC CAG ACT GCC TAC ATC TAC TCT TCT GCT CAC TTC GCG AAC AAC GTC TTT GTT GAT GGT AGC TTG AGC TAC TCT CAC TTC AAC AAC GAC CTG TCT GCA ACC ATG AGC AAC GGT ACT TAC GTT GAC GGT AGC ACC AAC TCC GAC GCT TGG GGC TTC GGT TTG AAA GCC GGT TAC GAC TTC AAA CTG GGT GAT GCT GGT TAC GTG ACT CCT TAC GGC AGC GTT TCT GGT CTG TTC CAG TCT GGT GAT GAC TAC CAG CTG AGC AAC GAC ATG AAA GTT GAC GGT CAG TCT TAC GAC AGC ATG CGT TAT GAA CTG GGT GTA GAT GCA GGT TAT ACC TTC ACC TAC AGC GAA GAT CAG GCT CTG ACT CCG TAC TTC AAA CTG GCT TAC GTC TAC GAC GAC TCT AAC AAC GAT AAC GAT GTG AAC GGC GAT TCC ATC GAT AAC GGT ACT GAA GGG TCT GCG GTA CGT GTT GGT CTG GGT ACT CAG TTT AGC TTC ACC AAG AAC TTC AGC GCC TAT ACC GAT GCT AAC TAC CTC GGT GGT GGT GAC GTA GAT CAA GAC TGG TCC GCG AAC GTG GGT GTT AAA TAT ACC TGG 
CFT073                               CTG GAT GTG AAC GTT GAT CAG TTT GAT ACT GAA GCT TTC CGT ACT GAC AAA CTG GAA CTG ACT AGC GGC AAC ATC GCT GAC CAC AAC GGC AAT GTA GTT TCT GGT GTG TTC GAT ATC AAC AGC AGC GAT TAT GTT CTG AAC GCT GAC CTG GTG AAC GAC CGT ACT TGG GAT ACC GCT CAG GCT AAC TAC GGC TAC GGT GTT GTT GCT ATG AAC TCT GAT GGT CAT CTG ACT ATC AAC GGT AAC GGC GAC GTA GAC AAC --- --- GGT ACT GAA CTG GAT AAC AGT TCT GTA GAC AAC --- --- --- --- GTT GTT GCG GCA ACC GGT AAC TAC AAA GTT CGT ATC GAC AAC GCA ACT GGC GCT GGC GCA ATC GCT GAT TAC AAA GAT AAA GAA ATT ATC TAT GTA AAC GAC GTC AAC AGC AAC GCG ACC TTC TCT GCT GCT AAC AAA GCT GAC CTG GGT GCA TAC ACC TAT CAG GCT GAG CAG CGT GGT AAC ACC GTT GTT CTG CAA CAG ATG GAG CTG ACC GAC TAC GCT AAC ATG GCG CTG AGC ATC CCG TCT GCG AAC ACC AAT ATC TGG AAC CTG GAA CAA GAC ACC GTT GGT ACT CGT CTG GCC AAC TCT CGT CAT GGT CTG GCT GAT AAC GGC GGC GCA TGG GTA AGC TAC TTC GGT GGT AAC TTC AAC GGC GAC AAC GGC ACC ATC AAC TAT GAT CAG GAT GTT AAC GGC ATC ATG GTC GGT GTT GAT ACC AAA ATT GAC GGT AAC AAC GCT AAG TGG ATC GTC GGT GCG GCT GCA GGC TTC GCT AAA GGT GAT ATG AAT GAC CGT TCT GGT CAG GTG GAT CAA GAC AGC CAG ACT GCC TAC ATC TAC TCT TCT GCT CAC TTC GCG AAC AAC GTC TTT GTT GAC GGT AGC TTG AGC TAC TCT CAC TTC AAC AAC GAC CTG TCT GCA ACC ATG AGC AAC GGT ACT TAC GTT GAC GGC AGC ACC AAC TCC GAC GCT TGG GGC TTT GGT TTG AAA GCC GGT TAC GAC TTC AAA CTG GGT GAT GCA GGT TAC GTG ACT CCT TAC GGC AGC GTT TCT GGT CTG TTC CAG TCT GGT GAT GAT TAC CAG CTG AGC AAC GAC ATG AAA GTT GAT GGT CAG TCT TAC GAC AGC ATG CGT TAT GAA CTG GGT GTA GAT GCA GGT TAT ACC TTC ACC TAC AGC GAA GAC CAG GCT CTG ACT CCG TAC TTC AAA CTG GCT TAC GTC TAC GAC GAC TCT AAC AAC GAT AAC GAT GTG AAC GGT GAT TCC ATC GAT AAC GGT ACT GAA GGG TCT GCG GTA CGT GTT GGT CTG GGT ACT CAG TTC AGC TTC ACC AAG AAC TTC AGC GCC TAT ACC GAT GCT AAC TAC CTC GGT GGT GGT GAC GTA GAT CAA GAC TGG TCC GCG AAT GTG GGT GTT AAA TAT ACC TGG 
O157                                 TTG GAT ATT AAC GTT GAT CAG TTC GAT GCT AAC CCG TTC CAG GCC GAT ACC CTG AAT CTG ACC AGT --- --- ACC ACT GAT ACT AAC GGC AAC ATT CAC GCT GGT GTA TTC GAT ATC CAT AGC AGT GAT TAC GTA ATG GAT ACC GAT CTG GTC AAC GAT CGT ACC AAC GAT ACT ACC AAG TCA AAC TAC GGT TAT GGC TTA ATC GCA ATG AAC TCT GAT GGT CAC CTG ACT ATT AAC GGT AAC GGC GAT AAC GAC AAC ACT GCT TCT ATC GAA GCT GGT CAG AAC GAA GTT GAT AAC AAC GGT GAC CAT GTT GCA GCC GCG ACC GGT AAC TAC AAA GTT CGT ATC GAC AAC GCT ACT GGT GCT GGT TCT ATC GCT GAC TAC AAC GGC AAC GAG CTG ATC TAC GTC AAC GAC AAA AAC AGC AAC GCG ACC TTC TCT GCT GCT AAC AAA GCT GAC CTG GGT GCA TAC ACC TAT CAG GCT GAA CAG CGC GGT AAC ACC GTT GTT CTG CAA CAG ATG GAG TTG ACC GAC TAC GCT AAC ATG GCG CTG AGC ATC CCA TCT GCG AAC ACC AAT ATC TGG AAC CTG GAA CAA GAC ACC GTT GGT ACT CGT TTG ACC AAC TCT CGT CAT GGC CTG GCT GAT AAC GGC GGC GCA TGG GTA AGC TAC TTC GGT GGT AAC TTC AAC GGC GAC AAC GGC ACC ATC AAC TAT GAT CAG GAT GTT AAC GGC ATC ATG GTC GGT GTT GAT ACC AAA ATT GAC GGT AAC AAC GCT AAG TGG ATC GTC GGT GCG GCT GCA GGC TTC GCT AAA GGT GAC ATG AAT GAC CGT TCT GGT CAG GTG GAT CAA GAC AGC CAG ACT GCC TAC ATC TAC TCT TCT GCT CAC TTC GCG AAC AAC GTC TTT GTT GAT GGT AGC TTG AGC TAC TCT CAC TTC AAC AAC GAC CTG TCT GCA ACC ATG AGC AAC GGT ACT TAC GTT GAC GGT AGC ACC AAC TCC GAC GCT TGG GGC TTC GGT TTG AAA GCC GGT TAC GAC TTC AAA CTG GGT GAT GCT GGT TAC GTG ACT CCT TAC GGC AGC ATT TCT GGT CTG TTC CAG TCT GGT GAT GAC TAC CAG CTG AGC AAC GAC ATG AAA GTT GAC GGT CAG TCT TAC GAC AGC ATG CGT TAT GAA CTG GGT GTA GAT GCA GGT TAT ACC TTC ACC TAC AGC GAA GAT CAG GCT CTG ACT CCG TAC TTC AAA CTG GCT TAC GTC TAC GAC GAC TCT AAC AAC GAT AAC GAT GTG AAC GGT GAT TCC ATC GAT AAC GGT ACT GAA GGG TCT GCG GTA CGT GTT GGT CTG GGT ACT CAG TTC AGC TTC ACC AAG AAC TTC AGC GCC TAT ACC GAT GCT AAC TAC CTC GGT GGT GGT GAC GTA GAT CAA GAC TGG TCC GCG AAC GTG GGT GTT AAA TAT ACC TGG 


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

ns = 3  	ls = 473
# site patterns = 127
    1   21    1   41   13   12    2    1    1    1   14    1    1    5    1
   17    1    2    2    1    1    1    1    1    1    2    1    1    1   32
    1   11    1   15    2    1    1    1    1    3    1    7    1    1    1
    1    1   20    2    1    1    1    1    1    8   14    1    9    1    9
    1   20    1    1    1    1    5    1    1    1    1    1    1    1    1
    1    1    1    1    1    1   16    9    1    2   16    1    1    2    1
    1    1    1    1    7    5    6    1    1    4    1    2    1    2    6
    1    1    1    4    4    1    2    3    4    3    2    1    2    3    1
    1    1    1    1    1    1    1

K12                                  TTG GAT GTG AAC GTT CAG TTC ACT GAA GCT TTC CGT ACT GAC AAA CTG GAA ACC AGC GGC AAC ATC GCT CAT GGT AAC GTA GTA TCT GGT GTG ATC CAT AGC TAC GTT CTG AAC GCT GAT GTG CGT TGG ACT TCC AAG TCT TAC GGT TAC GGT ATT GTT GCT ATG TCT CAC ACT ATC GGC GTA GAC --- --- GGT ACT GAA CTG GAT AAC AGC TCT GTA AAT --- --- --- --- GTT GCT GCA ACC AAA GCA GGC GCT GCT GAT AAA GAT GAA ATT GTA GTC GCG GCA TAT GAA CGC CAA GAG CTG CCG AAT TGG ACC CAT GGC GTA GTC ATT AAG GAC GTG GCC CAC TTT TTG TCC GCT CCT GTT CCG GGC GGG TTT CTC 
CFT073                               C.. ... ... ... ... ... ..T ... ... ... ... ... ... ... ... ... ... ..T ... ... ... ... ... ..C ..C ..T ... ..T ... ... ... ... A.C ... ..T ... ... ... ... ..C ... ... ... ..C G.T C.. G.. ... ..C ... ... G.. ... ... ... ... ..T ... ... ... ... ... --- --- ... ... ... ... ... ... ..T ... ... ..C --- --- --- --- ... ..G ... ... ... ... ... ... ..A ... ... ... ... ... ... ... ... ... ... ..G ..T ... ... ... ... ... ... G.. ... ..T ... ... ... ... ..T ... ... ... ... ... ... ..A ... ... ... ..T ... ..C ... 
O157                                 ... ... A.T ... ... ... ... G.. A.C C.G ... .AG G.C ..T .CC ... A.T ... ..T --- --- .C. A.. AC. ..C ... A.T CAC G.. ... ..A ... ... ... ... ..A A.. G.T A.C ... ..C ... AAC ... A.. ... ..A ... ... ..T ..C T.A A.C ..A ... ... ... ... ..T ... AAC ... ACT GCT TC. .TC ... GCT .G. C.G .A. GAA ..T ..C AAC GGT GAC CAT .CA ..C ..G ... ... ..T ..T ... T.. ..C ..C .GC ..G C.G ..C AAA ... ... ... ... ... ... ... T.. ..A ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... A.. ... ..T ... ..C ... 



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

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

#1: K12            
position  1:    T:0.16702    C:0.11991    A:0.29550    G:0.41756
position  2:    T:0.22484    C:0.20557    A:0.38544    G:0.18415
position  3:    T:0.34261    C:0.41542    A:0.08779    G:0.15418

#2: CFT073         
position  1:    T:0.16060    C:0.12206    A:0.29122    G:0.42612
position  2:    T:0.22484    C:0.20557    A:0.38544    G:0.18415
position  3:    T:0.35118    C:0.40257    A:0.08779    G:0.15846

#3: O157           
position  1:    T:0.16773    C:0.11890    A:0.31210    G:0.40127
position  2:    T:0.21231    C:0.21656    A:0.39278    G:0.17834
position  3:    T:0.34607    C:0.42251    A:0.07856    G:0.15287

Sums of codon usage counts
------------------------------------------------------------------------------
Phe F TTT       6 | Ser S TCT      49 | Tyr Y TAT      21 | Cys C TGT       0
      TTC      48 |       TCC      10 |       TAC      66 |       TGC       0
Leu L TTA       1 |       TCA       1 | *** * TAA       0 | *** * TGA       0
      TTG      10 |       TCG       0 |       TAG       0 | Trp W TGG      20
------------------------------------------------------------------------------
Leu L CTT       0 | Pro P CCT       3 | His H CAT       8 | Arg R CGT      24
      CTC       3 |       CCC       0 |       CAC      10 |       CGC       2
      CTA       0 |       CCA       1 | Gln Q CAA      12 |       CGA       0
      CTG      61 |       CCG       6 |       CAG      39 |       CGG       0
------------------------------------------------------------------------------
Ile I ATT      10 | Thr T ACT      40 | Asn N AAT      10 | Ser S AGT       3
      ATC      39 |       ACC      59 |       AAC     143 |       AGC      50
      ATA       0 |       ACA       0 | Lys K AAA      34 | Arg R AGA       0
Met M ATG      25 |       ACG       0 |       AAG       8 |       AGG       0
------------------------------------------------------------------------------
Val V GTT      50 | Ala A GCT      67 | Asp D GAT      84 | Gly G GGT     112
      GTC      16 |       GCC      12 |       GAC      81 |       GGC      42
      GTA      23 |       GCA      23 | Glu E GAA      24 |       GGA       0
      GTG      18 |       GCG      23 |       GAG       5 |       GGG       3
------------------------------------------------------------------------------

(Ambiguity data are not used in the counts.)


Codon position x base (3x4) table, overall

position  1:    T:0.16512    C:0.12028    A:0.29964    G:0.41495
position  2:    T:0.22064    C:0.20925    A:0.38790    G:0.18221
position  3:    T:0.34662    C:0.41352    A:0.08470    G:0.15516

Codon frequencies under model, for use in evolver:
  0.01285876  0.01534074  0.00314208  0.00575608
  0.01219508  0.01454896  0.00297991  0.00545899
  0.02260653  0.02697001  0.00000000  0.00000000
  0.01061885  0.01266848  0.00000000  0.00475341
  0.00936694  0.01117493  0.00228884  0.00419300
  0.00888349  0.01059816  0.00217071  0.00397659
  0.01646769  0.01964626  0.00402393  0.00737157
  0.00773528  0.00922833  0.00189014  0.00346261
  0.02333422  0.02783815  0.00570179  0.01044530
  0.02212987  0.02640134  0.00540750  0.00990618
  0.04102306  0.04894126  0.01002411  0.01836350
  0.01926955  0.02298892  0.00470857  0.00862579
  0.03231318  0.03855022  0.00789583  0.01446463
  0.03064540  0.03656053  0.00748830  0.01371806
  0.05680865  0.06777377  0.01388137  0.02542974
  0.02668443  0.03183502  0.00652043  0.01194498



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.0482 (0.0056 0.1154)
O157                 0.3182 (0.0512 0.1609) 0.2073 (0.0547 0.2639)


TREE #  1:  (1, (2, 3));   MP score: -1
This is a rooted tree.  Please check!
lnL(ntime:  4  np:  9):  -2411.448197     +0.000000
   4..1     4..5     5..2     5..3  
  0.00259  0.00000  0.09859  0.29505  1.56389  0.92228  0.00000  0.04174  3.78788
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.39622

(1: 0.002587, (2: 0.098586, 3: 0.295048): 0.000000);

(K12: 0.002587, (CFT073: 0.098586, O157: 0.295048): 0.000000);

Detailed output identifying parameters

kappa (ts/tv) =  1.56389


dN/dS for site classes (K=3)

p:   0.92228  0.00000  0.07772
w:   0.04174  1.00000  3.78788

dN & dS for each branch

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

   4..1       0.003    322.4   1096.6   0.3329   0.0006   0.0018    0.6    0.6
   4..5       0.000    322.4   1096.6   0.3329   0.0000   0.0000    0.0    0.0
   5..2       0.099    322.4   1096.6   0.3329   0.0226   0.0678   21.9   24.8
   5..3       0.295    322.4   1096.6   0.3329   0.0676   0.2030   65.4   74.1


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

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

     3 V      0.805         3.056
    10 T      0.534         2.043
    11 E      0.985*        3.731
    12 A      0.669         2.549
    14 R      0.969*        3.671
    15 T      0.514         1.969
    17 K      0.970*        3.675
    19 E      0.987*        3.738
    26 A      0.519         1.984
    28 H      0.994**       3.765
    32 V      0.716         2.725
    33 V      0.997**       3.777
    34 S      0.640         2.438
    40 H      0.538         2.057
    46 L      0.769         2.921
    56 W      1.000**       3.788
    59 S      0.993**       3.760
    60 K      0.550         2.101
    61 S      0.734         2.790
    67 I      0.999**       3.786
    84 V      0.997**       3.777
    89 G      0.994**       3.766
    90 T      0.553         2.115
    92 L      0.997**       3.777
    94 N      0.990*        3.750
    96 S      1.000**       3.787
   105 V      0.645         2.458
   123 A      0.644         2.455
   128 K      0.506         1.937
   130 K      0.506         1.937
   132 I      0.901         3.417
   138 V      1.000**       3.787
   352 V      0.508         1.944


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

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

     3 V      0.500         3.306 +- 3.088
    11 E      0.613         4.022 +- 2.838
    14 R      0.760         4.870 +- 2.796
    17 K      0.699         4.512 +- 2.860
    19 E      0.634         4.139 +- 2.831
    28 H      0.808         5.173 +- 2.645
    33 V      0.902         5.705 +- 2.360
    56 W      0.967*        6.046 +- 2.050
    59 S      0.780         5.013 +- 2.707
    67 I      0.949         5.952 +- 2.143
    84 V      0.901         5.702 +- 2.363
    89 G      0.805         5.141 +- 2.641
    92 L      0.923         5.832 +- 2.276
    94 N      0.705         4.563 +- 2.804
    96 S      0.952*        5.960 +- 2.129
   132 I      0.617         3.999 +- 3.027
   138 V      0.906         5.703 +- 2.322



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.999  0.001  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
w2:   0.001  0.017  0.092  0.180  0.212  0.186  0.135  0.089  0.055  0.034

Posterior for p0-p1 (see the ternary graph)

 0.000
 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.128
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.331 0.539

sum of density on p0-p1 =   1.000000
