
seed used = 140004485

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


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

ns = 6  	ls = 214
# site patterns = 64
   11    6    1    5    7   13    3    3   14    3    4    5    4   14    1
    4   21    2    1    2    3    6    2    5    5    3    1    1    5    3
    1    3    1    1    2    3    1    5    1    1    1    1    4    1    4
    1    4    1    1    3    1    3    1    1    1    1    1    1    1    1
    1    1    1    1

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

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

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

#1: K12            
position  1:    T:0.21963    C:0.24766    A:0.27103    G:0.26168
position  2:    T:0.49533    C:0.25701    A:0.13551    G:0.11215
position  3:    T:0.32243    C:0.20093    A:0.18224    G:0.29439

#2: S2457T         
position  1:    T:0.22430    C:0.24299    A:0.27570    G:0.25701
position  2:    T:0.49065    C:0.25701    A:0.13551    G:0.11682
position  3:    T:0.31308    C:0.21495    A:0.19159    G:0.28037

#3: Shig4337       
position  1:    T:0.22430    C:0.24299    A:0.27570    G:0.25701
position  2:    T:0.49065    C:0.25701    A:0.13551    G:0.11682
position  3:    T:0.31308    C:0.21495    A:0.19159    G:0.28037

#4: CFT073         
position  1:    T:0.22430    C:0.24299    A:0.27570    G:0.25701
position  2:    T:0.49533    C:0.25234    A:0.13551    G:0.11682
position  3:    T:0.31308    C:0.21495    A:0.19626    G:0.27570

#5: O157           
position  1:    T:0.21963    C:0.24766    A:0.27103    G:0.26168
position  2:    T:0.49533    C:0.25701    A:0.13551    G:0.11215
position  3:    T:0.32243    C:0.20093    A:0.17757    G:0.29907

#6: EDL933         
position  1:    T:0.21963    C:0.24766    A:0.27103    G:0.26168
position  2:    T:0.49533    C:0.25701    A:0.13551    G:0.11215
position  3:    T:0.32243    C:0.20093    A:0.17757    G:0.29907

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


Codon position x base (3x4) table, overall

position  1:    T:0.22196    C:0.24533    A:0.27336    G:0.25935
position  2:    T:0.49377    C:0.25623    A:0.13551    G:0.11449
position  3:    T:0.31776    C:0.20794    A:0.18614    G:0.28816

Codon frequencies under model, for use in evolver:
  0.03550002  0.02323163  0.02079536  0.03219365
  0.01842194  0.01205553  0.01079128  0.01670617
  0.00974291  0.00637587  0.00000000  0.00000000
  0.00823108  0.00538651  0.00000000  0.00746446
  0.03923687  0.02567707  0.02298434  0.03558246
  0.02036109  0.01332454  0.01192721  0.01846471
  0.01076848  0.00704702  0.00630801  0.00976553
  0.00909751  0.00595352  0.00532918  0.00825019
  0.04372108  0.02861159  0.02561112  0.03964902
  0.02268807  0.01484734  0.01329032  0.02057497
  0.01199916  0.00785239  0.00702892  0.01088159
  0.01013722  0.00663392  0.00593823  0.00919307
  0.04147898  0.02714433  0.02429773  0.03761574
  0.02152458  0.01408594  0.01260876  0.01951984
  0.01138382  0.00744971  0.00666846  0.01032356
  0.00961737  0.00629372  0.00563370  0.00872163



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.0898 (0.0063 0.0702)
Shig4337             0.0898 (0.0063 0.0702)-1.0000 (0.0000 0.0000)
CFT073               0.0546 (0.0042 0.0769) 0.1129 (0.0021 0.0185) 0.1129 (0.0021 0.0185)
O157                 0.0000 (0.0000 0.0061) 0.0993 (0.0063 0.0635) 0.0993 (0.0063 0.0635) 0.0598 (0.0042 0.0702)
EDL933               0.0000 (0.0000 0.0061) 0.0993 (0.0063 0.0635) 0.0993 (0.0063 0.0635) 0.0598 (0.0042 0.0702)-1.0000 (0.0000 0.0000)


TREE #  1:  ((1, (2, 3)), (4, (5, 6)));   MP score: 27
This is a rooted tree.  Please check!
lnL(ntime: 10  np: 15):   -983.638417     +0.000000
   7..8     8..1     8..9     9..2     9..3     7..10   10..4    10..11   11..5    11..6  
  0.00000  0.00519  0.07588  0.00000  0.00000  0.00000  0.07665  0.00000  0.00000  0.00000  2.94150  0.98419  0.00000  0.00000 13.63370
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.15774

((1: 0.005192, (2: 0.000004, 3: 0.000004): 0.075883): 0.000000, (4: 0.076649, (5: 0.000004, 6: 0.000004): 0.000000): 0.000000);

((K12: 0.005192, (S2457T: 0.000004, Shig4337: 0.000004): 0.075883): 0.000000, (CFT073: 0.076649, (O157: 0.000004, EDL933: 0.000004): 0.000000): 0.000000);

Detailed output identifying parameters

kappa (ts/tv) =  2.94150


dN/dS for site classes (K=3)

p:   0.98419  0.00000  0.01581
w:   0.00000  1.00000 13.63370

dN & dS for each branch

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

   7..8       0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0
   8..1       0.005    203.9    438.1   0.2156   0.0008   0.0037    0.8    0.4
   8..9       0.076    203.9    438.1   0.2156   0.0117   0.0544   11.1    5.1
   9..2       0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0
   9..3       0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0
   7..10      0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0
  10..4       0.077    203.9    438.1   0.2156   0.0119   0.0550   11.2    5.2
  10..11      0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0
  11..5       0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0
  11..6       0.000    203.9    438.1   0.2156   0.0000   0.0000    0.0    0.0


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

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

    70 A      1.000**       13.634
   114 G      1.000**       13.634
   125 L      1.000**       13.633


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

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

    70 A      0.931         6.456 +- 2.901
   114 G      0.926         6.429 +- 2.921



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.982  0.018  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
w2:   0.044  0.064  0.083  0.097  0.108  0.115  0.119  0.122  0.124  0.124

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.001
 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.033
 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.022 0.942

sum of density on p0-p1 =   1.000000
