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

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACG  >  NC_000913/3969888‑3970159
                                                                                                                                           |                                                                                                                                      
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                         >  1:376119/1‑139 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCggggggg                                                                                                                                <  2:247528/139‑2 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCggggggg                                                                                                                                <  1:601821/139‑2 (MQ=255)
                 tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCgggg                                                                                                                       <  1:379252/139‑1 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTc                                                                                                               >  2:311350/1‑139 (MQ=255)
                             cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCAt                                                                                                           <  1:672270/139‑1 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGTCATTGTCGGGGGGGCTGATCGGGGc                                                                                                                      <  1:646646/114‑1 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                                                      >  2:646646/1‑114 (MQ=255)
                                                       tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                                 <  1:311350/139‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                                                   <  1:261495/119‑1 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgc                                                                                                   >  2:261495/1‑119 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               >  1:50505/1‑139 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                               <  2:376119/139‑1 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTg                                                                                                                    >  1:550649/1‑92 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTg                                                                                                                    <  2:550649/92‑1 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  1:1052938/1‑111 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                     >  1:564446/1‑139 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCGGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  2:1052938/111‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCAt                                                                                                           <  1:720165/97‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCAt                                                                                                           >  2:720165/1‑97 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         <  2:523168/115‑1 (MQ=255)
                                                                       cGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                         >  1:523168/1‑115 (MQ=255)
                                                                          cGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCa                                                              <  1:55420/139‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  2:413601/70‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  2:835136/70‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  1:835136/1‑70 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  1:413601/1‑70 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGAGCTGATCGGGACTGGTGTCGCATTAAGCCGCCGTTGCTCGAAATAACAACGCTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg                            >  1:722650/1‑139 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                           >  2:270852/1‑74 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                           <  1:270852/74‑1 (MQ=255)
                                                                                                                  ccTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 <  1:692804/64‑1 (MQ=255)
                                                                                                                  ccTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                                 >  2:692804/1‑64 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                                       >  2:1024154/1‑81 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                                       <  1:1024154/81‑1 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTa           >  2:1033210/4‑139 (MQ=255)
                                                                                                                                  ggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTgg      >  1:833459/1‑139 (MQ=255)
                                                                                                                                   gCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGt     >  1:263857/1‑139 (MQ=255)
                                                                                                                                      ttGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTAcg  <  2:263857/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                      
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGGTACG  >  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: