Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,419,779 C→T 30.0% A78A (GCG→GCA wzyE ← predicted Wzy protein involved in ECA polysaccharide chain elongation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,419,7790CT30.0% 60.4 / 26.5 63A78A (GCG→GCAwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base C (18/26);  new base T (9/10);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAC  >  minE/2419717‑2419841
                                                              |                                                              
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCTCGTAGAAGCAGCCCGCAGa                                                           <  1:2736753/68‑1 (MQ=255)
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaa                                                          >  1:2551417/1‑69 (MQ=255)
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGa                                                           <  1:1441557/68‑1 (MQ=255)
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGa                                                           <  1:421042/68‑1 (MQ=255)
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGa                                                           <  1:3108685/68‑1 (MQ=255)
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGa                                                           <  1:1602996/68‑1 (MQ=255)
tCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGa                                                           <  1:1691744/68‑1 (MQ=255)
 cAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaa                                                          >  1:1531824/1‑68 (MQ=255)
 cAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaa                                                          >  1:2271070/1‑68 (MQ=255)
 cAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaa                                                          >  1:1964051/1‑68 (MQ=255)
  aGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaaa                                                         >  1:934159/1‑68 (MQ=255)
  aGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaaa                                                         >  1:3379892/1‑68 (MQ=255)
  aGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaaa                                                         >  1:3219796/1‑68 (MQ=255)
  aGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGaaa                                                         >  1:3098149/1‑68 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:1877820/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:879019/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:673072/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:669740/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:490287/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:3262695/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:2256135/68‑1 (MQ=255)
    cAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGc                                                       <  1:1693793/68‑1 (MQ=255)
                  aGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAg                                        <  1:2500907/69‑1 (MQ=255)
                   ggCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAg                                        <  1:479380/68‑1 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCTTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:1493226/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:1271858/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:369554/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:1508142/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:3043480/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:1634182/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:2093304/1‑69 (MQ=255)
                      gggTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:2089291/1‑69 (MQ=255)
                      ggcttttGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGAttt                                    >  1:1465754/4‑69 (MQ=255)
                                ggTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGgcgc                          <  1:2596823/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:2445886/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:1641417/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:2655802/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:894829/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:52918/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:692683/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:737864/69‑1 (MQ=255)
                                      aTAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    >  1:324606/1‑69 (MQ=255)
                                       tAGTAAACCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACAcc                    <  1:1460867/68‑1 (MQ=255)
                                         gTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCaa                  >  1:2215931/1‑68 (MQ=255)
                                         gTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAAc                 >  1:347372/1‑69 (MQ=255)
                                         gTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAAc                 >  1:1737687/1‑69 (MQ=255)
                                         gTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAAc                 >  1:535507/1‑69 (MQ=255)
                                         gTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAAc                 >  1:83102/1‑69 (MQ=255)
                                         gTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAAc                 >  1:1574787/1‑69 (MQ=255)
                                           aaaCCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACa                <  1:1909261/68‑1 (MQ=255)
                                           aaaCCGCGTAGAAGCAGCCTGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACa                <  1:1959756/68‑1 (MQ=255)
                                               cgcgTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCaaa           <  1:1727605/69‑1 (MQ=255)
                                                   tAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGa       <  1:2165794/69‑1 (MQ=255)
                                                   tAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGa       <  1:2223603/69‑1 (MQ=255)
                                                    aGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGa       <  1:1950922/68‑1 (MQ=255)
                                                       agcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTg                                 >  1:3424591/1‑39 (MQ=255)
                                                       agcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAAt    >  1:52690/1‑68 (MQ=255)
                                                       agcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAAt    >  1:3410501/1‑68 (MQ=255)
                                                        gcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAc  <  1:2028398/69‑1 (MQ=255)
                                                        gcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAc  <  1:3343413/69‑1 (MQ=255)
                                                        gcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAc  <  1:3339226/69‑1 (MQ=255)
                                                        gcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAc  <  1:2256910/69‑1 (MQ=255)
                                                        gcagcCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAc  <  1:2519922/69‑1 (MQ=255)
                                                              |                                                              
TCAGCAACGCGTTTGCGTAGGCGGGTTTTGTAGGTGACATAGTAAACCGCGTAGAAGCAGCCCGCAGAAAGCAACGCCTGCAACAAGATTTCTGGCGGCGCGACACCAACATCAAAGCGAAATAC  >  minE/2419717‑2419841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: