Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,137,907 Δ1 bp coding (537/942 nt) flgJ → putative peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,137,9070A.98.2% 228.5 / ‑2.1 55coding (537/942 nt)flgJputative peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (36/18);  minor base C (0/1);  total (36/19)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACG  >  NC_000913/1137790‑1138040
                                                                                                                     |                                                                                                                                     
cTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACg                                                                                                                >  1:266571/1‑140 (MQ=255)
     atgatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCaaa                                                                                                           >  1:458672/1‑140 (MQ=255)
     atgatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCaaa                                                                                                           >  2:29897/1‑140 (MQ=255)
       gatTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATc                                                                                                         >  2:106278/1‑140 (MQ=255)
           cGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATccgcc                                                                                                     >  2:477165/1‑140 (MQ=255)
                           tAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGAAAACCCGCCGCGAAAACGGcgagcc                                                                                     <  2:887630/140‑1 (MQ=255)
                                  ttCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCtat                                                                              >  1:792884/1‑140 (MQ=255)
                                   tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCtata                                                                             >  1:1134147/1‑140 (MQ=255)
                                   tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCtata                                                                             >  1:763522/1‑140 (MQ=255)
                                   tCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCtata                                                                             >  2:778569/1‑140 (MQ=255)
                                    ccTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccgag                                                                                  >  2:1060899/1‑134 (MQ=255)
                                    ccTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGcgagccgag                                                                                  <  1:1060899/134‑1 (MQ=255)
                                       cgcgCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCt                                                                         >  2:1197544/1‑140 (MQ=255)
                                         cgcAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCTCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGt                                                                       >  1:316046/1‑140 (MQ=255)
                                         cgcAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGt                                                                       >  2:1386524/1‑140 (MQ=255)
                                         cgcAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGCCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGt                                                                       <  2:792884/140‑1 (MQ=255)
                                            aaCTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTg                                                                    <  2:266571/140‑1 (MQ=255)
                                            aaCTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGACAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTg                                                                    >  1:1341565/1‑140 (MQ=255)
                                               tctcGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGgtg                                                                 <  2:763522/140‑1 (MQ=255)
                                                  cGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCa                                                              <  2:762750/140‑1 (MQ=255)
                                                  cGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCa                                                              >  2:168539/1‑140 (MQ=255)
                                                      gcccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGc                                                          >  2:904342/1‑140 (MQ=255)
                                                        ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCaa                                                             <  1:1309797/135‑1 (MQ=255)
                                                        ccgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCaa                                                             >  2:1309797/1‑135 (MQ=255)
                                                         cgcccAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGCCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCctc                                                       >  1:293781/1‑140 (MQ=255)
                                                             cAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGc                                                   >  2:465084/1‑140 (MQ=255)
                                                               aCTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCaa                                                 >  2:179365/1‑140 (MQ=255)
                                                                cTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAAc                                                <  1:477165/140‑1 (MQ=255)
                                                                    cAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGa                                            >  2:56315/1‑140 (MQ=255)
                                                                    cAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGa                                            <  1:1197544/140‑1 (MQ=255)
                                                                     aaGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGAAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGaa                                           <  1:828828/140‑1 (MQ=255)
                                                                          aGCAAGGCGGTGTGCCACATCATTTGATCCCCGCTCAGACGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGc                                      <  2:1172614/140‑1 (MQ=255)
                                                                          aGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGc                                      >  1:1172614/1‑140 (MQ=255)
                                                                           gCAAAGCGGTGTGCCACATCATTTGATCCTCTCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGcc                                     <  2:316046/140‑1 (MQ=255)
                                                                           gCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGcc                                     <  2:1134147/140‑1 (MQ=255)
                                                                                 cGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTAc                               >  2:1278193/1‑140 (MQ=255)
                                                                                 cGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTAc                               >  2:545833/1‑140 (MQ=255)
                                                                                     gtgCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGaa                           >  2:879299/1‑140 (MQ=255)
                                                                                      tgCCACATCATTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGaaa                          >  1:804520/1‑140 (MQ=255)
                                                                                           catcatTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATcacc                     >  2:638784/1‑140 (MQ=255)
                                                                                               atTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGa                 <  2:228754/140‑1 (MQ=255)
                                                                                               atTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGa                 <  1:168539/140‑1 (MQ=255)
                                                                                               atTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGa                 <  1:179365/140‑1 (MQ=255)
                                                                                               atTTGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGa                 <  1:291405/140‑1 (MQ=255)
                                                                                                tttGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGAc                >  2:847521/1‑140 (MQ=255)
                                                                                                tttGATCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGAc                >  2:276694/1‑140 (MQ=255)
                                                                                                tttGACCCTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGAc                <  2:86070/140‑1 (MQ=255)
                                                                                                  tGATCCTCGCCCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTg              <  2:561347/140‑1 (MQ=255)
                                                                                                      ccTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGa                 <  2:584066/133‑1 (MQ=255)
                                                                                                      ccTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGa                 >  1:584066/1‑133 (MQ=255)
                                                                                                      ccTCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAata          >  1:842843/1‑140 (MQ=255)
                                                                                                        tCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATg        >  1:349398/1‑140 (MQ=255)
                                                                                                        tCGCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATg        >  2:1249937/1‑140 (MQ=255)
                                                                                                          gCTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGaa      >  2:350256/1‑140 (MQ=255)
                                                                                                           cTCAGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGaaa     >  1:969959/1‑140 (MQ=255)
                                                                                                              aGGCGGCCTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACg  >  1:744468/1‑140 (MQ=255)
                                                                                                                     |                                                                                                                                     
CTACGATGATTCGCTGCCGGGTGACAGTAAAGCATTCCTCGCGCAACTCTCGCTGCCCGCCCAACTGGCAAGCCAGCAAAGCGGTGTGCCACATCATTTGATCCTCGCTCAGGCGGCACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACG  >  NC_000913/1137790‑1138040

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