Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,922,483 G→A intergenic (‑43/+574) atpI ← / ← rsmG ATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,922,4830GA98.1% 189.2 / ‑3.6 53intergenic (‑43/+574)atpI/rsmGATP synthase, membrane‑bound accessory factor/16S rRNA m(7)G527 methyltransferase, SAM‑dependent; glucose‑inhibited cell‑division protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/37);  total (15/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCT  >  NC_000913/3922336‑3922630
                                                                                                                                                   |                                                                                                                                                   
gAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCatat                                                                                                                                                  >  1:908627/1‑151 (MQ=255)
gAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCAta                                                                                                                                                   <  2:368484/150‑1 (MQ=255)
  aGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATaa                                                                                                                                                <  1:300116/151‑1 (MQ=255)
    gggTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATaaaa                                                                                                                                              <  2:54006/151‑1 (MQ=255)
          ttGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActtt                                                                                                                                          <  1:755059/149‑1 (MQ=255)
           tGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttct                                                                                                                                        <  1:305683/150‑1 (MQ=255)
            gAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttctt                                                                                                                                       <  1:299072/150‑1 (MQ=255)
            gAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAActttctt                                                                                                                                       <  2:588339/150‑1 (MQ=255)
             aGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTa                                                                                                                                     <  2:797304/151‑1 (MQ=255)
                       gCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTagag                                                                                                                           <  1:63102/151‑1 (MQ=255)
                             cACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGc                                                                                                                     >  2:747402/1‑151 (MQ=255)
                               cTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCAt                                                                                                                   >  1:873577/1‑151 (MQ=255)
                               cTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCAt                                                                                                                   >  1:641284/1‑151 (MQ=255)
                                   cTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCaa                                                                                                                >  2:511051/1‑150 (MQ=255)
                                     aTCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGGGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGCCTTAGAGTCAGCCATCaaaa                                                                                                              >  2:986053/1‑150 (MQ=255)
                                       caccacCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCTAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGc                                                                                                           <  2:174601/151‑1 (MQ=255)
                                         ccaccaGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCg                                                                                                          <  1:82460/150‑1 (MQ=255)
                                                  aCTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtat                                                                                                >  1:93965/1‑151 (MQ=255)
                                                     gAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtata                                                                                               <  2:376180/149‑1 (MQ=255)
                                                     gAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtata                                                                                               <  1:545203/149‑1 (MQ=255)
                                                             aGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGccc                                                                                     <  1:58470/151‑1 (MQ=255)
                                                                  cTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTg                                                                                >  1:369709/1‑151 (MQ=255)
                                                                      cGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGAttt                                                                            <  2:928238/151‑1 (MQ=255)
                                                                       gAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTc                                                                           <  1:834835/151‑1 (MQ=255)
                                                                        aGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCa                                                                          <  1:450008/151‑1 (MQ=255)
                                                                          cAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCaaa                                                                        <  1:399355/151‑1 (MQ=255)
                                                                             cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtat                                                                                                >  2:819014/1‑124 (MQ=255)
                                                                             cGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATtat                                                                                                <  1:819232/124‑1 (MQ=255)
                                                                                   gACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTa                                                               <  1:643797/151‑1 (MQ=255)
                                                                                         cGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaa                                                          <  2:438586/150‑1 (MQ=255)
                                                                                          gAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCaaa                                                         >  1:812115/1‑150 (MQ=255)
                                                                                              gACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTg                                                    >  1:209766/1‑151 (MQ=255)
                                                                                                    gACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATaaa                                              <  2:683584/151‑1 (MQ=255)
                                                                                                    gACATCACGGTTTTAACTCCTGATCCATTGGAGGTATAGCGCGTGTCAAATAGAACTTTCTTTAAGGCATAGAGGCAATCATCAAAAAGCGGTCAAATTATCCGGGGCGCCCCCGAGATTGGCGAAAGTAAGTAGCCAAAAGGTCAATata                                              <  2:498164/151‑3 (MQ=255)
                                                                                                     aCATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAAt                                             <  1:911452/151‑1 (MQ=255)
                                                                                                      cATTACGTTTTTCACTCCTGCCCCCTTCGAGGTATGCCGCGCGTCATATAACACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGCTCAAATTATTCTGTGCGCCCCCGTGATTTAAAACAATATGTAGCTAAAAGGtttcta                                                >  2:407099/1‑142 (MQ=255)
                                                                                                      cATCACGTTGTTCACTCCTGCTCACTTCGAGGTATGCCGCGTGTCATATAAAACGTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGATTAACAAAt                                             <  2:914567/150‑1 (MQ=255)
                                                                                                         cACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGtt                                          <  1:584035/150‑1 (MQ=255)
                                                                                                          aCGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTa                                        <  1:486036/151‑1 (MQ=255)
                                                                                                           cGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTaa                                       <  2:34160/151‑1 (MQ=255)
                                                                                                               tttCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAata                                    <  1:948467/150‑1 (MQ=255)
                                                                                                                  cACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAtttt                                <  2:820522/151‑1 (MQ=255)
                                                                                                                   aCTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttttt                               <  1:632824/151‑1 (MQ=255)
                                                                                                                    cTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATAttt                                 <  2:939378/148‑1 (MQ=255)
                                                                                                                    cTCCTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTc                              >  1:273818/1‑151 (MQ=255)
                                                                                                                      ccTGCTCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCa                            <  2:804491/151‑1 (MQ=255)
                                                                                                                           tCCCTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGc                        >  2:426640/1‑150 (MQ=255)
                                                                                                                            cccTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCAt                      <  1:753554/151‑1 (MQ=255)
                                                                                                                              cTTCGAGGTATGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATAc                    >  2:431430/1‑151 (MQ=255)
                                                                                                                                 cGAGGTATGTCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAAt                 <  2:437001/151‑1 (MQ=255)
                                                                                                                                       aTGCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGa            <  1:8324/150‑1 (MQ=255)
                                                                                                                                         gCCGCGTGTCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACt          <  1:602593/150‑1 (MQ=255)
                                                                                                                                                gtCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCt  >  1:563801/1‑151 (MQ=255)
                                                                                                                                                gtCATATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCt  <  1:209180/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAGGGGTCTTTGAGGCTGAACAGCAATCCACTTGCTATCACCACCAGTAACTGAACGAGCAGAAGCTTCCGAGCAACGTTTCGACTCACGAGCGACACAGACATCACGTTTTTCACTCCTGCTCCCTTCGAGGTATGCCGCGTGTCGTATAAAACTTTCTTTAAGGCTTAGAGTCAAGCATCAAAAAGCGGTCAAATTATACGGTGCGCCCCCGTGATTTCAAACAATAAGTAGCCAAAAGGTGAATAAATGTTTAAATATTTTTCCAGTGCATACAATTGCGACTTTTCTGCT  >  NC_000913/3922336‑3922630

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