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.G98.0% 172.9 / ‑3.3 51coding (996/1047 nt)wzzEEntobacterial Common Antigen (ECA) polysaccharide chain length modulation protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (25/25);  total (25/26)
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.58e-01

CCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGAC  >  NC_000913/3969903‑3970146
                                                                                                                            |                                                                                                                          
ccTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCgg                                                                                                            >  2:823350/1‑139 (MQ=255)
  tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCgggg                                                                                                          <  2:195742/139‑1 (MQ=255)
  tGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCgggg                                                                                                          >  1:543576/1‑139 (MQ=255)
   gAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGc                                                                                                         <  1:628892/139‑1 (MQ=255)
          ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTc                                                                                                  >  2:979564/1‑139 (MQ=255)
          ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTc                                                                                                  >  1:779028/1‑139 (MQ=255)
          ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggg‑g                                                                                                                      <  2:15399/118‑1 (MQ=255)
          ggTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggg‑g                                                                                                                      >  1:15399/1‑118 (MQ=255)
              cAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCAt                                                                                              >  2:912593/1‑139 (MQ=255)
                aCCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggg‑gggg                                                                                                                   >  1:689301/1‑114 (MQ=255)
                aCCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggg‑gggg                                                                                                                   <  2:689301/115‑2 (MQ=255)
                     ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccg                                                                                       >  1:204296/1‑139 (MQ=255)
                     ggATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACccg                                                                                       >  1:470184/1‑139 (MQ=255)
                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggg‑gg                                                                                                                     >  2:494797/1‑102 (MQ=255)
                           gcgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTC‑‑ggg‑gg                                                                                                                     <  1:494797/102‑1 (MQ=255)
                            cgTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGc                                                                                >  2:504925/1‑139 (MQ=255)
                             gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                               <  1:504925/139‑1 (MQ=255)
                             gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                               <  2:3003/139‑1 (MQ=255)
                             gTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCt                                                                               <  2:779028/139‑1 (MQ=255)
                              ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                              <  1:912593/139‑1 (MQ=255)
                              ttttCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTc                                                                              <  1:491167/139‑1 (MQ=255)
                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAAcc                                                                                         <  1:273075/118‑1 (MQ=255)
                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAAcc                                                                                         >  2:273075/1‑118 (MQ=255)
                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                            >  2:781763/1‑130 (MQ=255)
                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGa                                                                            <  1:781763/130‑1 (MQ=255)
                                         aTCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATaacaac                                                                   >  2:325684/1‑139 (MQ=255)
                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                  <  1:400601/139‑1 (MQ=255)
                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                  <  2:543576/139‑1 (MQ=255)
                                          tCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                  <  1:823350/139‑1 (MQ=255)
                                                  tGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCGGGGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACtgctgc                                                           >  2:753192/1‑139 (MQ=255)
                                                            ggAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCaa                                                <  1:979564/139‑1 (MQ=255)
                                                             gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGgtg                                                                                                    >  1:435788/1‑86 (MQ=255)
                                                             gaagaaCCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGgtg                                                                                                    <  2:435788/86‑1 (MQ=255)
                                                                   ccGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCt                                         <  2:470184/139‑1 (MQ=255)
                                                                       tAAAACGCGATAGCCCACGTCGTGCCTTCTTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                  <  2:868593/110‑1 (MQ=255)
                                                                       tAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAacaaca                                                                  >  1:868593/1‑110 (MQ=255)
                                                                          aaCGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCt                                  <  2:204296/139‑1 (MQ=255)
                                                                                    cccACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACActgct                                                             >  1:351115/1‑102 (MQ=255)
                                                                                    cccACGTCGTGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACActgct                                                             <  2:351115/102‑1 (MQ=255)
                                                                                        cgtcgtGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACCGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAAt                    <  1:715294/139‑1 (MQ=255)
                                                                                             tGCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtg               <  1:325684/139‑1 (MQ=255)
                                                                                              gCCTTCCTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGAtgt              >  2:436269/1‑139 (MQ=255)
                                                                                                   ccTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGGTAATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACGGCTGTGGTGATCGCAAAGGCGCTCGTCGCTTATTCGAAGAGAATCGATGGGAAAg         >  1:687588/1‑139 (MQ=255)
                                                                                                   ccTGATGATTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAg         >  1:267400/1‑139 (MQ=255)
                                                                                                       atgatTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact     >  2:273289/1‑139 (MQ=255)
                                                                                                       atgatTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTact     >  2:168700/1‑139 (MQ=255)
                                                                                                         gatTATGCGGGGCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactga   <  1:436269/139‑1 (MQ=255)
                                                                                                          tttATGCGGTGCATTGTCG‑GGG‑GGGCTCATCGGGGCTGGTGTCGCCTTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGCCGCTCGCCTCTTATTCGAAGAGAATCGATGTGAAAGCactgac  <  1:574924/138‑1 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                    <  1:868349/66‑1 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACaa                                                                    >  2:868349/4‑69 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                          >  1:660029/4‑79 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                    <  1:20185/98‑1 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTcgccg                                    >  2:20185/4‑101 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                 >  1:54569/4‑120 (MQ=255)
                                                                                                              tgcggggCATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGa                 <  2:54569/117‑1 (MQ=255)
                                                                                                              tgcggggCATTGTC‑‑GGGTGGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCg                                                          <  2:660029/76‑1 (MQ=255)
                                                                                                                     cATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg    <  2:954794/126‑1 (MQ=255)
                                                                                                                     cATTGTCG‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAACAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTactg    >  1:954794/1‑126 (MQ=255)
                                                                                                                            |                                                                                                                          
CCTGAATGTTGGTCCAACCCTGGATCCGCGTTTTCAGACCTATCGCTATTTGCGTACGCCGGAAGAACCGGTAAAACGCGATAGCCCACGTCGTGCCTTCCTGATGATTATGTGGGGCATTGTC‑‑GGG‑GGGCTGATCGGGGCTGGTGTCGCATTAACCCGCCGTTGCTCGAAATAGCAACACTGCTGCGGTGAGCGCAAAGGCGCTCGCCGCTTATTCGAAGAGAATCGATGTGAAAGTACTGAC  >  NC_000913/3969903‑3970146

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