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% 140.9 / 7.9 51intergenic (+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 (28/16);  new base C (3/4);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATG  >  NC_000913/1113401‑1113694
                                                                                                                                                   |                                                                                                                                                  
cAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATg                                                                                                                                                   >  3:496947/1‑149 (MQ=255)
 aaTAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGc                                                                                                                                                  >  1:338496/1‑149 (MQ=255)
              ttGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATATATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttatt                                                                                                                                     <  3:534412/149‑4 (MQ=25)
               tGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttattg                                                                                                                                    <  3:424410/149‑5 (MQ=255)
                  cTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGaa                                                                                                                                 <  4:496947/149‑1 (MQ=255)
                      tGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                             <  2:8312/149‑1 (MQ=255)
                         gcTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTa                                                                                                                          >  2:165819/1‑149 (MQ=255)
                                   tCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGt                                                                                                                >  1:341307/1‑149 (MQ=255)
                                      ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTa                                                                                                             >  4:608509/1‑149 (MQ=255)
                                       gCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTaa                                                                                                            >  3:425601/1‑149 (MQ=255)
                                                     tGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGcatc                                                                                              >  3:101972/1‑149 (MQ=255)
                                                        aaTTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATg                                                                                           <  1:165819/149‑1 (MQ=255)
                                                                       cAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                            <  3:608509/149‑1 (MQ=255)
                                                                        aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                                                                           >  1:162667/1‑149 (MQ=255)
                                                                        aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTg                                                                           >  3:284327/1‑149 (MQ=255)
                                                                         aGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTgg                                                                          >  3:15019/1‑149 (MQ=255)
                                                                          gACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                       >  1:8866/1‑119 (MQ=255)
                                                                          gACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                       <  3:526846/119‑1 (MQ=255)
                                                                          gACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                       >  4:526846/1‑119 (MQ=255)
                                                                          gACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                       <  2:8866/119‑1 (MQ=255)
                                                                            cTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGc                                                                                                                  >  1:435395/1‑106 (MQ=255)
                                                                            cTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGc                                                                                                                  <  2:435395/106‑1 (MQ=255)
                                                                             tttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                                                                      >  1:31799/1‑149 (MQ=255)
                                                                             tttGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                                                                      >  3:314216/1‑149 (MQ=255)
                                                                                  gggTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTACGCCGGTATTCTTTTACGAACGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGgcag                                                                 >  1:547905/1‑149 (MQ=255)
                                                                                         aTAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCattatt                                                          >  1:253074/1‑149 (MQ=255)
                                                                                          tAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                              >  4:168332/1‑80 (MQ=255)
                                                                                          tAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAaga                                                                                                                              <  3:168332/80‑1 (MQ=255)
                                                                                             ggCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                     <  1:374523/70‑1 (MQ=255)
                                                                                             ggCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTa                                                                                                                                     >  2:374523/1‑70 (MQ=255)
                                                                                             ggCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                                                                      <  4:200338/133‑1 (MQ=255)
                                                                                             ggCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                                                                      >  3:200338/1‑133 (MQ=255)
                                                                                              gCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAtt                                                                                                          <  2:603630/96‑1 (MQ=255)
                                                                                              gCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAAtt                                                                                                          >  1:603630/1‑96 (MQ=255)
                                                                                                gTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                             >  3:199820/1‑75 (MQ=255)
                                                                                                gTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                             <  4:199820/75‑1 (MQ=255)
                                                                                                                         acGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCaa                          <  4:314216/149‑1 (MQ=255)
                                                                                                                         acGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCaa                          <  2:31799/149‑1 (MQ=255)
                                                                                                                               cAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAata                    >  1:485013/1‑149 (MQ=255)
                                                                                                                                    aaaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATc               >  3:459325/1‑149 (MQ=255)
                                                                                                                                     aaaTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATcc              >  2:609250/1‑149 (MQ=255)
                                                                                                                                       aTACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCg            >  1:60596/1‑149 (MQ=255)
                                                                                                                                       aTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGa                                                                              >  1:177288/1‑83 (MQ=38)
                                                                                                                                       aTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGa                                                                              <  2:177288/83‑1 (MQ=38)
                                                                                                                                        tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAAc                         >  3:129620/1‑135 (MQ=255)
                                                                                                                                        tACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAAc                         <  4:129620/135‑1 (MQ=255)
                                                                                                                                           cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgtgt        >  1:538409/1‑149 (MQ=255)
                                                                                                                                           cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgtgt        >  3:563937/1‑149 (MQ=255)
                                                                                                                                            ggCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGtgtg       <  2:547905/149‑1 (MQ=255)
                                                                                                                                               gTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCaa    >  2:588634/1‑149 (MQ=255)
                                                                                                                                                 tATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATg  <  2:338496/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAATAAAAACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATG  >  NC_000913/1113401‑1113694

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: ‑