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.G95.1% 135.4 / ‑0.0 41coding (996/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (20/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGG  >  NC_000913/3969888‑3970155
                                                                                                                                           |                                                                                                                                  
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                     >  2:251999/1‑139 (MQ=255)
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                     >  2:382612/1‑139 (MQ=255)
ggCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                     >  1:48366/1‑139 (MQ=255)
ggCCATGTTAAACACCCTAAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTc                                                                                                                                     >  2:718698/1‑139 (MQ=255)
  ccATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgg                                                                                                                                  <  2:48366/139‑1 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCggggggg                                                                                                                             <  1:623920/139‑2 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCggggggg                                                                                                                             <  1:382612/139‑2 (MQ=255)
            cacCCTGAATGTTGGTCCAACCCTGGATCCGCGTTCTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGAt                                                                                                                        <  1:782117/139‑1 (MQ=255)
                cTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCggg                                                                                                                    >  1:121477/1‑139 (MQ=255)
                 tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCgggg                                                                                                                   >  1:1210293/1‑139 (MQ=255)
                        tGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGgtgt                                                                                                            <  1:718698/139‑1 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTc                                                                                                           >  1:901104/1‑139 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgggggggctga                                                                                                                         >  1:40896/1‑120 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgggggggctga                                                                                                                         <  2:40896/125‑6 (MQ=255)
                                cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCggggggg                                                                                                                             >  2:1091272/1‑113 (MQ=255)
                                cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCggggggg                                                                                                                             <  1:1091272/114‑2 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                          >  2:292376/1‑139 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                         >  2:673928/1‑139 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                        <  1:251999/139‑1 (MQ=255)
                                                       tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             <  1:292376/139‑1 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                                                  <  1:311314/101‑1 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGc                                                                                                                  >  2:311314/1‑101 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                     <  1:620748/130‑1 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                     >  2:620748/1‑130 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                           >  2:1006010/1‑139 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                 <  2:121477/139‑1 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             >  2:1102690/1‑117 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             <  1:1102690/117‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgggggg                                                                                                                              >  1:714831/1‑63 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgggggg                                                                                                                              <  2:257002/63‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgggggg                                                                                                                              >  1:257002/1‑63 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCgggggg                                                                                                                              <  2:714831/63‑1 (MQ=255)
                                                                                       aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                             <  2:1210293/139‑1 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                          <  2:901104/139‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                              <  2:591630/69‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                              >  1:591630/1‑69 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGACTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                       >  1:861092/1‑75 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGACTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                       <  2:861092/75‑1 (MQ=255)
                                                                                                              ccTTCCTGTTGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                             <  2:1002169/68‑1 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                             >  1:1002169/1‑68 (MQ=255)
                                                                                                               cTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTTTTCGAAGAGAATCGATGTGa                     <  1:673928/139‑1 (MQ=255)
                                                                                                               cTTCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                     <  1:612089/139‑1 (MQ=255)
                                                                                                                 tCCTGATGATTATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaaa                   >  2:424782/1‑139 (MQ=255)
                                                                                                                          ttATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact          >  2:556106/1‑139 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                        >  1:643069/1‑108 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                        <  2:643069/108‑1 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg         >  2:295483/1‑139 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTCACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg             >  2:917771/4‑133 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                           >  1:357637/4‑71 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                           <  2:357637/68‑1 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg             <  1:917771/130‑1 (MQ=255)
                                                                                                                                ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                             >  1:1218371/1‑114 (MQ=255)
                                                                                                                                ggggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                             <  2:1218371/114‑1 (MQ=255)
                                                                                                                                  ggCATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCAATAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGTCGAGAATCGATGTGAAAGTACTGACTGTATTTgg  >  2:868396/1‑139 (MQ=255)
                                                                                                                                    cATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             >  2:833973/1‑62 (MQ=255)
                                                                                                                                    cATTGTCGGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             <  1:833973/62‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                  
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTCGGGGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGG  >  NC_000913/3969888‑3970155

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