
seed used = 139993109

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


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

ns = 5  	ls = 311
# site patterns = 71
    1    1   13    4    4    1   17    3   10    9    7    6    7    5   11
    3    5    6   13    3   11    5    4   10    6    4    1    7    6    1
    3    1    7   11    5    2    6    3    3    2   12    7    1    1    1
    5    3    6    1    5    2    3    2    2    2    3    4    2    5    1
    1    1    2    2    1    1    4    1    1    1    1

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

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

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

#1: K12            
position  1:    T:0.16077    C:0.29904    A:0.19936    G:0.34084
position  2:    T:0.36013    C:0.23473    A:0.24759    G:0.15756
position  3:    T:0.27974    C:0.24759    A:0.15756    G:0.31511

#2: S2457T         
position  1:    T:0.16077    C:0.29904    A:0.20579    G:0.33441
position  2:    T:0.36013    C:0.23473    A:0.24759    G:0.15756
position  3:    T:0.27974    C:0.24759    A:0.15756    G:0.31511

#3: Shig4337       
position  1:    T:0.16077    C:0.29904    A:0.20579    G:0.33441
position  2:    T:0.36013    C:0.23473    A:0.24759    G:0.15756
position  3:    T:0.27974    C:0.24759    A:0.15756    G:0.31511

#4: O157           
position  1:    T:0.16399    C:0.29582    A:0.20257    G:0.33762
position  2:    T:0.36013    C:0.23473    A:0.24759    G:0.15756
position  3:    T:0.27653    C:0.25723    A:0.14469    G:0.32154

#5: EDL933         
position  1:    T:0.16077    C:0.29260    A:0.20579    G:0.34084
position  2:    T:0.36334    C:0.23473    A:0.24437    G:0.15756
position  3:    T:0.27653    C:0.25402    A:0.14148    G:0.32797

Sums of codon usage counts
------------------------------------------------------------------------------
Phe F TTT      20 | Ser S TCT      18 | Tyr Y TAT      15 | Cys C TGT       5
      TTC      25 |       TCC       7 |       TAC       5 |       TGC      10
Leu L TTA      10 |       TCA      24 | *** * TAA       0 | *** * TGA       0
      TTG      57 |       TCG      20 |       TAG       0 | Trp W TGG      35
------------------------------------------------------------------------------
Leu L CTT      30 | Pro P CCT      15 | His H CAT      55 | Arg R CGT      35
      CTC      30 |       CCC      10 |       CAC      69 |       CGC      25
      CTA      10 |       CCA      15 | Gln Q CAA      15 |       CGA      10
      CTG      88 |       CCG      30 |       CAG      20 |       CGG       5
------------------------------------------------------------------------------
Ile I ATT      62 | Thr T ACT      18 | Asn N AAT      35 | Ser S AGT       0
      ATC      29 |       ACC      29 |       AAC      10 |       AGC      15
      ATA      16 |       ACA       5 | Lys K AAA      10 | Arg R AGA       5
Met M ATG      51 |       ACG      20 |       AAG      10 |       AGG       2
------------------------------------------------------------------------------
Val V GTT      13 | Ala A GCT      32 | Asp D GAT      45 | Gly G GGT      35
      GTC      20 |       GCC      58 |       GAC      15 |       GGC      33
      GTA      33 |       GCA      28 | Glu E GAA      53 |       GGA       2
      GTG      67 |       GCG      36 |       GAG      27 |       GGG      28
------------------------------------------------------------------------------


Codon position x base (3x4) table, overall

position  1:    T:0.16141    C:0.29711    A:0.20386    G:0.33762
position  2:    T:0.36077    C:0.23473    A:0.24695    G:0.15756
position  3:    T:0.27846    C:0.25080    A:0.15177    G:0.31897

Codon frequencies under model, for use in evolver:
  0.01659096  0.01494336  0.00904265  0.01900489
  0.01079447  0.00972251  0.00588336  0.01236503
  0.01135638  0.01022861  0.00000000  0.00000000
  0.00724561  0.00652607  0.00000000  0.00829982
  0.03053794  0.02750531  0.01664424  0.03498111
  0.01986871  0.01789561  0.01082914  0.02275954
  0.02090298  0.01882716  0.01139285  0.02394429
  0.01333653  0.01201212  0.00726887  0.01527695
  0.02095352  0.01887269  0.01142040  0.02400219
  0.01363286  0.01227902  0.00743038  0.01561640
  0.01434252  0.01291820  0.00781717  0.01642931
  0.00915083  0.00824208  0.00498752  0.01048224
  0.03470221  0.03125603  0.01891390  0.03975126
  0.02257808  0.02033592  0.01230584  0.02586312
  0.02375338  0.02139450  0.01294642  0.02720942
  0.01515515  0.01365014  0.00826008  0.01736017



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

K12                 
S2457T               0.3366 (0.0029 0.0085)
Shig4337             0.3366 (0.0029 0.0085)-1.0000 (0.0000 0.0000)
O157                 0.0473 (0.0014 0.0304) 0.1419 (0.0043 0.0304) 0.1419 (0.0043 0.0304)
EDL933               0.1955 (0.0070 0.0356) 0.2768 (0.0099 0.0356) 0.2768 (0.0099 0.0356) 1.1037 (0.0055 0.0050)


TREE #  1:  ((1, (2, 3)), (4, 5));   MP score: 13
This is a rooted tree.  Please check!
lnL(ntime:  8  np: 13):  -1358.341127     +0.000000
   6..7     7..1     7..8     8..2     8..3     6..9     9..4     9..5  
  0.01680  0.00348  0.01037  0.00000  0.00000  0.00872  0.00000  0.01811  2.49209  0.99460  0.00000  0.16207 43.18192
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.05749

((1: 0.003478, (2: 0.000004, 3: 0.000004): 0.010373): 0.016797, (4: 0.000004, 5: 0.018112): 0.008716);

((K12: 0.003478, (S2457T: 0.000004, Shig4337: 0.000004): 0.010373): 0.016797, (O157: 0.000004, EDL933: 0.018112): 0.008716);

Detailed output identifying parameters

kappa (ts/tv) =  2.49209


dN/dS for site classes (K=3)

p:   0.99460  0.00000  0.00540
w:   0.16207  1.00000 43.18192

dN & dS for each branch

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

   6..7       0.017    282.7    650.3   0.3943   0.0038   0.0097    2.7    2.5
   7..1       0.003    282.7    650.3   0.3943   0.0008   0.0020    0.6    0.5
   7..8       0.010    282.7    650.3   0.3943   0.0024   0.0060    1.7    1.5
   8..2       0.000    282.7    650.3   0.3943   0.0000   0.0000    0.0    0.0
   8..3       0.000    282.7    650.3   0.3943   0.0000   0.0000    0.0    0.0
   6..9       0.009    282.7    650.3   0.3943   0.0020   0.0050    1.4    1.3
   9..4       0.000    282.7    650.3   0.3943   0.0000   0.0000    0.0    0.0
   9..5       0.018    282.7    650.3   0.3943   0.0041   0.0104    3.0    2.7


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

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

     1 H      1.000**       43.162


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

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

     1 H      0.947         6.492 +- 2.838
     2 S      0.533         3.983 +- 3.717



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.588  0.269  0.099  0.031  0.009  0.003  0.001  0.000  0.000  0.000
w2:   0.057  0.064  0.079  0.095  0.108  0.117  0.122  0.123  0.121  0.116

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.001
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.008
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.004 0.028
 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.013 0.085
 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.001 0.003 0.007 0.039 0.220
 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.001 0.001 0.003 0.007 0.019 0.103 0.448

sum of density on p0-p1 =   1.000000
