
seed used = 140016377

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


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

ns = 5  	ls = 298
# site patterns = 77
   10    6   12   16    2    7    6    6   15    1    6    4   13   19    7
    8    3    2    4    1    7    3   10    2    3    2    4    1    8    2
    4    5    4    1    2    2    5    7    4    1    1    1    1    1    1
    3    1    4    8    2    2   10    1    1    5    1    7    2    1    3
    1    1    1    1    1    1    1    1    1    1    1    1    1    1    1
    1    1

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

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

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

#1: K12            
position  1:    T:0.12416    C:0.24832    A:0.24497    G:0.38255
position  2:    T:0.29866    C:0.26510    A:0.26510    G:0.17114
position  3:    T:0.21812    C:0.20470    A:0.18456    G:0.39262

#2: S2457T         
position  1:    T:0.12416    C:0.24832    A:0.24161    G:0.38591
position  2:    T:0.29866    C:0.26846    A:0.26174    G:0.17114
position  3:    T:0.20470    C:0.20805    A:0.19128    G:0.39597

#3: Shig4337       
position  1:    T:0.12416    C:0.24832    A:0.24161    G:0.38591
position  2:    T:0.29866    C:0.26846    A:0.26174    G:0.17114
position  3:    T:0.20470    C:0.20805    A:0.19128    G:0.39597

#4: CFT073         
position  1:    T:0.12416    C:0.24832    A:0.24497    G:0.38255
position  2:    T:0.29866    C:0.26510    A:0.26174    G:0.17450
position  3:    T:0.19799    C:0.22483    A:0.18121    G:0.39597

#5: O157           
position  1:    T:0.12416    C:0.24832    A:0.24497    G:0.38255
position  2:    T:0.29866    C:0.26510    A:0.26510    G:0.17114
position  3:    T:0.21477    C:0.20470    A:0.18792    G:0.39262

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


Codon position x base (3x4) table, overall

position  1:    T:0.12416    C:0.24832    A:0.24362    G:0.38389
position  2:    T:0.29866    C:0.26644    A:0.26309    G:0.17181
position  3:    T:0.20805    C:0.21007    A:0.18725    G:0.39463

Codon frequencies under model, for use in evolver:
  0.00789661  0.00797303  0.00710695  0.01497809
  0.00704484  0.00711302  0.00634036  0.01336248
  0.00695612  0.00702343  0.00000000  0.00000000
  0.00454277  0.00458673  0.00000000  0.00861661
  0.01579322  0.01594606  0.01421390  0.02995618
  0.01408969  0.01422604  0.01268072  0.02672495
  0.01391223  0.01404687  0.01252101  0.02638836
  0.00908554  0.00917346  0.00817699  0.01723322
  0.01549443  0.01564438  0.01394499  0.02938944
  0.01382312  0.01395690  0.01244081  0.02621934
  0.01364903  0.01378112  0.01228413  0.02588913
  0.00891365  0.00899991  0.00802229  0.01690718
  0.02441547  0.02465175  0.02197392  0.04631063
  0.02178189  0.02199268  0.01960370  0.04131533
  0.02150756  0.02171570  0.01935680  0.04079499
  0.01404575  0.01418168  0.01264118  0.02664162



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.0910 (0.0045 0.0497)
Shig4337             0.0910 (0.0045 0.0497)-1.0000 (0.0000 0.0000)
CFT073               0.0204 (0.0015 0.0736) 0.0514 (0.0045 0.0881) 0.0514 (0.0045 0.0881)
O157                 0.0000 (0.0000 0.0177) 0.1119 (0.0045 0.0404) 0.1119 (0.0045 0.0404) 0.0204 (0.0015 0.0736)


TREE #  1:  ((1, (2, 3)), (4, 5));   MP score: 28
This is a rooted tree.  Please check!
lnL(ntime:  8  np: 13):  -1344.570460     +0.000000
   6..7     7..1     7..8     8..2     8..3     6..9     9..4     9..5  
  0.00000  0.01025  0.04279  0.00000  0.00000  0.00000  0.06111  0.00401  4.31202  0.99645  0.00000  0.00000 27.28836
Note: Branch length is defined as number of nucleotide substitutions per codon (not per neucleotide site).

tree length =   0.11817

((1: 0.010253, (2: 0.000004, 3: 0.000004): 0.042786): 0.000000, (4: 0.061108, 5: 0.004014): 0.000000);

((K12: 0.010253, (S2457T: 0.000004, Shig4337: 0.000004): 0.042786): 0.000000, (CFT073: 0.061108, O157: 0.004014): 0.000000);

Detailed output identifying parameters

kappa (ts/tv) =  4.31202


dN/dS for site classes (K=3)

p:   0.99645  0.00000  0.00355
w:   0.00000  1.00000 27.28836

dN & dS for each branch

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

   6..7       0.000    261.6    632.4   0.0968   0.0000   0.0000    0.0    0.0
   7..1       0.010    261.6    632.4   0.0968   0.0009   0.0095    2.5    0.6
   7..8       0.043    261.6    632.4   0.0968   0.0038   0.0395   10.3    2.4
   8..2       0.000    261.6    632.4   0.0968   0.0000   0.0000    0.0    0.0
   8..3       0.000    261.6    632.4   0.0968   0.0000   0.0000    0.0    0.0
   6..9       0.000    261.6    632.4   0.0968   0.0000   0.0000    0.0    0.0
   9..4       0.061    261.6    632.4   0.0968   0.0055   0.0564   14.8    3.5
   9..5       0.004    261.6    632.4   0.0968   0.0004   0.0037    1.0    0.2


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

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

   296 N      1.000**       27.288


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

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

   296 N      0.928         4.567 +- 2.823



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.998  0.002  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000
w2:   0.201  0.169  0.135  0.107  0.088  0.075  0.065  0.058  0.053  0.049

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.004
 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.005 0.991

sum of density on p0-p1 =   1.000000
