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

GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGG  >  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)
  ccATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gg                                                                                                                                  <  2:48366/139‑1 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                             <  1:623920/139‑2 (MQ=255)
       ttAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                             <  1:382612/139‑2 (MQ=255)
            cacCCTGAATGTTGGTCCAACCCTGGATCCGCGTTCTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGAt                                                                                                                        <  1:782117/139‑1 (MQ=255)
                cTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCggg                                                                                                                    >  1:121477/1‑139 (MQ=255)
                 tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCgggg                                                                                                                   >  1:1210293/1‑139 (MQ=255)
                        tGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGgtgt                                                                                                            <  1:718698/139‑1 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTc                                                                                                           >  1:901104/1‑139 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggggctga                                                                                                                         >  1:40896/1‑120 (MQ=255)
                         ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggggctga                                                                                                                         <  2:40896/125‑6 (MQ=255)
                                cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                             >  2:1091272/1‑113 (MQ=255)
                                cccTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑ggggggg                                                                                                                             <  1:1091272/114‑2 (MQ=255)
                                          gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTg                                                                                          >  2:292376/1‑139 (MQ=255)
                                           cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                         >  2:673928/1‑139 (MQ=255)
                                            gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                                        <  1:251999/139‑1 (MQ=255)
                                                       tATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             <  1:292376/139‑1 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGc                                                                                                                  <  1:311314/101‑1 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGc                                                                                                                  >  2:311314/1‑101 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                     <  1:620748/130‑1 (MQ=255)
                                                        aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                                     >  2:620748/1‑130 (MQ=255)
                                                         tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                           >  2:1006010/1‑139 (MQ=255)
                                                                   cGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGt                                                                 <  2:121477/139‑1 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             >  2:1102690/1‑117 (MQ=255)
                                                                             aagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             <  1:1102690/117‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                              >  1:714831/1‑63 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                              <  2:257002/63‑1 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                              >  1:257002/1‑63 (MQ=255)
                                                                                  ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑gggggg                                                                                                                              <  2:714831/63‑1 (MQ=255)
                                                                                       aaaaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                             <  2:1210293/139‑1 (MQ=255)
                                                                                                          cgtGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                          <  2:901104/139‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                              <  2:591630/69‑1 (MQ=255)
                                                                                                            tGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgcc                                                                                              >  1:591630/1‑69 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGGACTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                       >  1:861092/1‑75 (MQ=255)
                                                                                                             gCCTTCCTGATGATTATGCGGGGCATTGTC‑GGGGGGACTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                                       <  2:861092/75‑1 (MQ=255)
                                                                                                              ccTTCCTGTTGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                             <  2:1002169/68‑1 (MQ=255)
                                                                                                              ccTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACccgccg                                                                                             >  1:1002169/1‑68 (MQ=255)
                                                                                                               cTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTTTTCGAAGAGAATCGATGTGa                     <  1:673928/139‑1 (MQ=255)
                                                                                                               cTTCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGa                     <  1:612089/139‑1 (MQ=255)
                                                                                                                 tCCTGATGATTATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGaaa                   >  2:424782/1‑139 (MQ=255)
                                                                                                                          ttATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgact          >  2:556106/1‑139 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                        >  1:643069/1‑108 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTAt                                        <  2:643069/108‑1 (MQ=255)
                                                                                                                           tATGCGGGGCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactgactg         >  2:295483/1‑139 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTCACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg             >  2:917771/4‑133 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                           >  1:357637/4‑71 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                           <  2:357637/68‑1 (MQ=255)
                                                                                                                             tgcggggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg             <  1:917771/130‑1 (MQ=255)
                                                                                                                                ggggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                             >  1:1218371/1‑114 (MQ=255)
                                                                                                                                ggggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                             <  2:1218371/114‑1 (MQ=255)
                                                                                                                                  ggCATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCAATAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGTCGAGAATCGATGTGAAAGTACTGACTGTATTTgg  >  2:868396/1‑139 (MQ=255)
                                                                                                                                    cATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             >  2:833973/1‑62 (MQ=255)
                                                                                                                                    cATTGTCGGGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                             <  1:833973/62‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                  
GGCCATGTTAAACACCCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑GGGGGG‑CTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGACTGTATTTGG  >  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: ‑