Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,113,548 T→C 13.3% intergenic (+142/‑31) opgH → / → yceK OPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,113,5480TC13.3% 120.8 / 9.8 45intergenic (+142/‑31)opgH/yceKOPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (21/18);  new base C (3/3);  total (24/21)
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.96e-01

GGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAA  >  NC_000913/1113439‑1113692
                                                                                                             |                                                                                                                                                
ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTa                                                                                                           >  1:58116/1‑149 (MQ=255)
           ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                            >  4:1547062/1‑121 (MQ=255)
           ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                            <  3:1547062/121‑1 (MQ=255)
           ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAt                                                                                                                                          >  4:654789/1‑107 (MQ=255)
           ttCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAt                                                                                                                                          <  3:654789/107‑1 (MQ=255)
               tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                   >  4:855662/1‑110 (MQ=255)
               tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                   <  3:855662/110‑1 (MQ=255)
                   aTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATgg                                                                                        <  1:1450483/149‑1 (MQ=255)
                    tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCaaa                                                                                                                   >  1:538416/1‑121 (MQ=255)
                    tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCaaa                                                                                                                   <  2:538416/121‑1 (MQ=255)
                    tttATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGt                                                                                       >  3:710014/1‑149 (MQ=255)
                      tATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGa                                                                                     <  2:58116/149‑1 (MQ=255)
                       aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGaaa                                                                                                                              >  3:691771/1‑107 (MQ=255)
                       aTTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGaaa                                                                                                                              <  4:691771/107‑1 (MQ=255)
                        ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                   >  4:391082/1‑101 (MQ=255)
                        ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                   <  3:391082/101‑1 (MQ=255)
                                 cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                          <  4:710014/149‑1 (MQ=255)
                                  aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGt                                                                                                                         >  2:1697361/1‑101 (MQ=255)
                                  aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGt                                                                                                                         <  1:1697361/101‑1 (MQ=255)
                                  aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcat                                                                                             <  3:1246953/129‑1 (MQ=255)
                                  aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcat                                                                                             >  4:1246953/1‑129 (MQ=255)
                                  aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGCCGTTATCCCGGTTTTTTTTTACGAAAAAGGAATCAAATGCGTTTAATTTTGGTGAGAATAAAGGTCACCTTGCGGAggg                                                                         >  1:1826052/1‑147 (MQ=37)
                                                      aGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattatt                                                        <  2:1446458/146‑1 (MQ=255)
                                                      aGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattatt                                                        >  1:1446458/1‑146 (MQ=255)
                                                                                acaacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGa                                                                                                 <  2:1749199/79‑1 (MQ=255)
                                                                                acaacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGa                                                                                                 >  1:1749199/1‑79 (MQ=255)
                                                                                          aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggg                                       >  4:627044/1‑127 (MQ=255)
                                                                                          aaTAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCggg                                       <  3:627044/127‑1 (MQ=255)
                                                                                              aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                                                                                   >  1:1135960/1‑63 (MQ=255)
                                                                                              aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTgtggt                                                                                                   <  2:1135960/63‑1 (MQ=255)
                                                                                              aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATc             >  1:1574596/1‑149 (MQ=255)
                                                                                                aaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATccc           >  3:1430498/1‑149 (MQ=255)
                                                                                                 aTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCAc                                                <  1:140475/111‑1 (MQ=255)
                                                                                                 aTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTGTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCAc                                                >  2:140475/1‑111 (MQ=255)
                                                                                                  tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATg                            <  3:392401/130‑1 (MQ=255)
                                                                                                  tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCgg         >  4:772924/1‑149 (MQ=255)
                                                                                                     cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc                                                              <  4:1540250/93‑1 (MQ=255)
                                                                                                     cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc                                                              >  3:1540250/1‑93 (MQ=255)
                                                                                                     cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc                                                              <  2:1711941/93‑1 (MQ=255)
                                                                                                     cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcagc                                                              >  1:1711941/1‑93 (MQ=255)
                                                                                                      ggCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATAc                                           >  3:232356/1‑111 (MQ=255)
                                                                                                      ggCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATAc                                           <  4:232356/111‑1 (MQ=255)
                                                                                                       gCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAat                   <  3:1687968/134‑1 (MQ=255)
                                                                                                       gCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAat                   >  4:1687968/1‑134 (MQ=255)
                                                                                                         gTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCaa  >  3:965773/1‑149 (MQ=255)
                                                                                                             |                                                                                                                                                
GGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAA  >  NC_000913/1113439‑1113692

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: