Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,032 (G)6→7 coding (1002/1047 nt) 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
*NC_0009133,970,0261.G97.1% 120.3 / ‑3.1 35coding (996/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (17/17);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACG  >  NC_000913/3969888‑3970159
                                                                                                                                           |                                                                                                                                      
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                         >  1:376119/1‑139 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggg‑gg                                                                                                                                <  2:247528/139‑2 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggg‑gg                                                                                                                                <  1:601821/139‑2 (MQ=255)
                 tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCgggg                                                                                                                       <  1:379252/139‑1 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTc                                                                                                               >  2:311350/1‑139 (MQ=255)
                             cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCAt                                                                                                           <  1:672270/139‑1 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGTCATTGTCGGGGGG‑GCTGATCGGGGc                                                                                                                      <  1:646646/114‑1 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGc                                                                                                                      >  2:646646/1‑114 (MQ=255)
                                                       tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                 <  1:311350/139‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                                                   <  1:261495/119‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                                                   >  2:261495/1‑119 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               >  1:50505/1‑139 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               <  2:376119/139‑1 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTg                                                                                                                    >  1:550649/1‑92 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTg                                                                                                                    <  2:550649/92‑1 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  1:1052938/1‑111 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                     >  1:564446/1‑139 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCGGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  2:1052938/111‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCAt                                                                                                           <  1:720165/97‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCAt                                                                                                           >  2:720165/1‑97 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         <  2:523168/115‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         >  1:523168/1‑115 (MQ=255)
                                                                          cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCa                                                              <  1:55420/139‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  2:413601/70‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  2:835136/70‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  1:835136/1‑70 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  1:413601/1‑70 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGAGCTGATCGGGACTGGTGTCGCATTAAGCCGCCGTTGCTCGAAATAACAACGCTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                            >  1:722650/1‑139 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                           >  2:270852/1‑74 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                           <  1:270852/74‑1 (MQ=255)
                                                                                                                  ccTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  1:692804/64‑1 (MQ=255)
                                                                                                                  ccTGATGATTATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  2:692804/1‑64 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                                       >  2:1024154/1‑81 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                                       <  1:1024154/81‑1 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa           >  2:1033210/4‑139 (MQ=255)
                                                                                                                                  ggCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTgg      >  1:833459/1‑139 (MQ=255)
                                                                                                                                   gCATTGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGt     >  1:263857/1‑139 (MQ=255)
                                                                                                                                      ttGTCGGGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg  <  2:263857/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGG‑GCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACG  >  NC_000913/3969888‑3970159

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑