Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,429,480 C→T 16.0% G35S (GGC→AGC)  wzzE ← Entobacterial Common Antigen (ECA) polysaccharide chain length modulation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,429,4800CT16.0% 211.0 / 28.8 88G35S (GGC→AGC) wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base C (41/32);  new base T (10/4);  total (52/36)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TCGAGCTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTCAGCGT  >  minE/2429411‑2429548
                                                                     |                                                                    
tCGAGCTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGcc                                                                      >  1:2020785/1‑70 (MQ=255)
tCGAGCTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCa                                                                     <  1:1673006/71‑1 (MQ=255)
tCGAGCTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCa                                                                     <  1:2209017/71‑1 (MQ=255)
   aGCTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCaat                                                                   >  1:65153/1‑70 (MQ=255)
     cTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCaataat                                                                >  1:1136961/1‑71 (MQ=255)
     cTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCaataat                                                                >  1:479359/1‑71 (MQ=255)
       ccACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATc                                                               >  1:3834058/1‑70 (MQ=255)
        cACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATcc                                                              <  1:197803/70‑1 (MQ=255)
          cTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCTAATAATCCAt                                                            >  1:1855627/1‑70 (MQ=255)
          cTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCTAATAATCCAt                                                            >  1:2118356/1‑70 (MQ=255)
           tCCTGACGAGCAAAAAACGTATACGCCAGAGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAg                                                          >  1:1290007/1‑71 (MQ=255)
             cTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCTAATAATCCATAGCt                                                        >  1:2803543/1‑71 (MQ=255)
             cTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCa                                                                     <  1:3475067/58‑1 (MQ=255)
             cTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCt                                                        >  1:395363/1‑71 (MQ=255)
             cTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCt                                                        >  1:2830878/1‑71 (MQ=255)
             cTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCt                                                        >  1:1829799/1‑71 (MQ=255)
               gACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCaat                                                                   >  1:2653848/1‑58 (MQ=255)
                aCGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCaat                                                                   <  1:1073602/57‑1 (MQ=255)
                aCGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTcc                                                     >  1:2819177/1‑71 (MQ=255)
                aCGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTcc                                                     >  1:2861312/1‑71 (MQ=255)
                aCGAGAAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCtt                                                       <  1:1517323/69‑1 (MQ=255)
                     cAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGcc                                                >  1:2033191/1‑71 (MQ=255)
                       aaaaaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCaa                                                                    >  1:3527424/1‑49 (MQ=255)
                       aaaaaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCca                                              >  1:65479/1‑71 (MQ=255)
                        aaaaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCa                                                             >  1:383896/1‑55 (MQ=255)
                        aaaaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCca                                              >  1:3219029/1‑70 (MQ=255)
                        aaaaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGACAATAATCCATAGCTTCCCAGCCcac                                             <  1:1823927/71‑1 (MQ=255)
                         aaaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCTaa                                                                    >  1:2633466/1‑47 (MQ=255)
                          aaCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCcac                                             >  1:955246/1‑69 (MQ=255)
                           aCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAg                                          <  1:2585339/71‑1 (MQ=255)
                           aCGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAg                                          <  1:363480/71‑1 (MQ=255)
                             gTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCa                                                             <  1:3108514/50‑1 (MQ=255)
                             gTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAgg                                         <  1:1137368/70‑1 (MQ=255)
                              tataCGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCcaca                                            <  1:2722100/66‑1 (MQ=255)
                               ataCGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTAc                                      >  1:182714/1‑71 (MQ=255)
                                taCGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACg                                     <  1:981299/71‑1 (MQ=255)
                                taCGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCAAGCCCACAAGGTAc                                      >  1:3341055/1‑70 (MQ=255)
                                taCGCCAGCGCGATTAACGCAAACGCAAGCCCCATGCCAATAATCCATAGCTTCCCa                                                   >  1:2925559/1‑57 (MQ=255)
                                  cGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGa                                    >  1:294967/1‑70 (MQ=255)
                                  cGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGa                                    <  1:1348779/70‑1 (MQ=255)
                                  cGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGa                                    <  1:16093/70‑1 (MQ=255)
                                  cGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGa                                    >  1:1687953/1‑70 (MQ=255)
                                   gCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATa                                                           <  1:2485293/46‑1 (MQ=255)
                                    ccAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGaaa                                  >  1:1295693/1‑70 (MQ=255)
                                      aGCGCGATTAACGCAAACGCCAGCCCCATGCCaat                                                                   <  1:3916643/35‑1 (MQ=255)
                                       gcgcGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAaacaac                              >  1:2830535/1‑71 (MQ=255)
                                       gcgcGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAaacaac                              <  1:107407/71‑1 (MQ=255)
                                          cgATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGaa                                   >  1:1207598/1‑63 (MQ=255)
                                          cgATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCa                           <  1:1289746/71‑1 (MQ=255)
                                           gATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCAc                          >  1:1693980/1‑71 (MQ=255)
                                            aTTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACg                         <  1:1605714/71‑1 (MQ=255)
                                            aTTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACg                         <  1:948920/71‑1 (MQ=255)
                                            aTTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACg                         <  1:1870037/71‑1 (MQ=255)
                                             ttAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACg                         >  1:765288/1‑70 (MQ=255)
                                              tAACGCAAACGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGaa                       >  1:3576378/1‑71 (MQ=255)
                                              tAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCa                                                   >  1:2341114/1‑43 (MQ=255)
                                              tAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGt                                        >  1:2583744/1‑54 (MQ=255)
                                              tAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGaaa                                  >  1:1502394/1‑60 (MQ=255)
                                                aCGCAAACGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAg                                          <  1:1132734/50‑1 (MQ=255)
                                                aCGCAAACGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAg                                          <  1:339835/50‑1 (MQ=255)
                                                aCGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCa                                                   >  1:2760659/1‑41 (MQ=255)
                                                  gCAAACGCCAGCCCCATGCCAATAATCCATAGCTTccc                                                    >  1:3186164/1‑38 (MQ=255)
                                                  gCAAACGCCAGCCCCATGCCAATAATCCATACCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATc                   <  1:414261/71‑1 (MQ=255)
                                                   cAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTAc                                      >  1:1432917/1‑51 (MQ=255)
                                                   cAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATc                   >  1:2413410/1‑70 (MQ=255)
                                                     aaCGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAg                >  1:2345720/1‑71 (MQ=255)
                                                      aCGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGt               >  1:3865639/1‑71 (MQ=255)
                                                      aCGCAAGCCCCATGCAAATAATCCATAGCTTCCCAGcca                                               >  1:2366814/1‑38 (MQ=38)
                                                       cGCCAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGtt              >  1:2304665/1‑71 (MQ=255)
                                                        gCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTc             <  1:1823058/71‑1 (MQ=255)
                                                         ccAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGt               >  1:3057185/1‑68 (MQ=255)
                                                          cAGCCCCATGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCa            <  1:3583621/70‑1 (MQ=255)
                                                          cAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAgg                                         <  1:2306898/41‑1 (MQ=255)
                                                           aGCCCCATGCCAATAATCCATAGCTTCCCAGCCcaca                                            <  1:2076191/37‑1 (MQ=255)
                                                           aGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAtt          <  1:1086085/71‑1 (MQ=255)
                                                            gCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTAc                                      >  1:335699/1‑42 (MQ=255)
                                                            gCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAttt         <  1:3500308/71‑1 (MQ=255)
                                                            gCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAttt         >  1:1347947/1‑71 (MQ=255)
                                                            gCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAttt         >  1:1455557/1‑71 (MQ=255)
                                                            gCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAttt         >  1:2669261/1‑71 (MQ=255)
                                                             ccccATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAtttt        <  1:1568313/71‑1 (MQ=255)
                                                              cccATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAtttt        <  1:3177075/70‑1 (MQ=255)
                                                              cccATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTc       <  1:1312170/71‑1 (MQ=255)
                                                              cccATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTc       <  1:1980441/71‑1 (MQ=255)
                                                               ccATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAatat                    >  1:567615/1‑57 (MQ=255)
                                                                cATGCCAATAATCCATAGCTTCCCAGCCCACAAgg                                         >  1:2669534/1‑35 (MQ=255)
                                                                cATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTCAg     >  1:1242549/1‑71 (MQ=255)
                                                                 cccccAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTCa      >  1:756140/4‑69 (MQ=255)
                                                                 aTGCCAATAATCCATAGCTTCCCAGCCCACAAGGTa                                       >  1:833488/1‑36 (MQ=255)
                                                                  tGCTAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCAtt          <  1:744157/64‑1 (MQ=255)
                                                                    ccAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTCAGCGt  <  1:1975328/70‑1 (MQ=255)
                                                                     cAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAg                <  1:252561/55‑1 (MQ=255)
                                                                     |                                                                    
TCGAGCTCCACTCCTGACGAGCAAAAAACGTATACGCCAGCGCGATTAACGCAAACGCCAGCCCCATGCCAATAATCCATAGCTTCCCAGCCCACAAGGTACGAAACAACCCACGAATATCCAGTTCATTTTCAGCGT  >  minE/2429411‑2429548

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