Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,400,882 T→G 8.0% intergenic (+39/‑130) ybgF → / → queE TolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,400,8820TG 8.0% 124.0 / 3.6 50intergenic (+39/‑130)ybgF/queETolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base T (22/24);  new base G (2/2);  total (24/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 = 4.57e-01

ATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACGC  >  NC_002947/1400746‑1401029
                                                                                                                                        |                                                                                                                                                   
aTACCGACAAGGTAAAAGGCATCCTGCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCCCTCTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTc                                                                                                                                        <  2:72770/150‑1 (MQ=255)
   ccGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGtt                                                                                                                                     <  1:39450/150‑1 (MQ=255)
           gTAAAAGGCATCCTTCCGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCAc                                                                                                                             <  2:315973/150‑1 (MQ=255)
             aaaaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTGAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGGTTTTTCGTTAGAATCACTg                                                                                                                           >  1:458126/1‑150 (MQ=255)
                   cATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCttt                                                                                                                     <  2:631709/150‑1 (MQ=255)
                     tCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGTGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCttttt                                                                                                                   <  1:1154591/150‑1 (MQ=255)
                        ttCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGAGCCCTAAAGAAGCTCACACGGTTAGCGCTGTAGCGCGCGAAAGAAACCCGCGCCAGTCGCGGGGTTTTTCGTTAGAATCACTGCCCTTTTTTCg                                                                                                                >  1:448861/1‑150 (MQ=255)
                             cagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAac                                                                                                           >  2:987034/1‑150 (MQ=255)
                                   gTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTcgcc                                                                                                                                                 <  2:442344/106‑2 (MQ=255)
                                   gTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTcgcc                                                                                                                                                 >  1:442344/1‑105 (MQ=255)
                                       cccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCtgctgc                                                                                                 <  1:1052223/150‑1 (MQ=255)
                                       cccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCtgctgc                                                                                                 >  2:1052223/1‑150 (MQ=255)
                                        ccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACCCGCTGCTGCg                                                                                                >  1:654079/1‑150 (MQ=255)
                                            tACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAACCCCGCGCCAGTCGCGGGCCTTCTCGTCAAAATCACTGCCGTTGTTTGTTAAACACGCTGGTGCGGATa                                                                                            >  2:132367/1‑150 (MQ=255)
                                                ccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATaa                                                                                           <  1:432204/147‑1 (MQ=255)
                                                ccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATaa                                                                                           >  2:432204/1‑147 (MQ=255)
                                                ccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGc                                                                                        >  2:640503/1‑150 (MQ=255)
                                                 cGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACg                                                                                                        >  1:420182/1‑133 (MQ=255)
                                                 cGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACg                                                                                                        <  2:420182/133‑1 (MQ=255)
                                                       ctctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTcgc                                                                                                                                                  <  2:4507/85‑1 (MQ=255)
                                                       ctctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTcgc                                                                                                                                                  >  1:4507/1‑85 (MQ=255)
                                                        tctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCTCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                                                                <  2:154277/150‑1 (MQ=255)
                                                                  aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACa                                                                      <  2:1013106/150‑1 (MQ=255)
                                                                        cccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                                <  1:547422/150‑1 (MQ=255)
                                                                         ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTg                                                               >  1:1110541/1‑150 (MQ=255)
                                                                         ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTg                                                               >  1:355035/1‑150 (MQ=255)
                                                                          cAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCAAGCCCGCGTCAGTTCCTGa                                                              <  2:379996/150‑1 (MQ=255)
                                                                            gCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                                                                <  1:1086262/130‑1 (MQ=255)
                                                                            gCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                                                                >  2:1086262/1‑130 (MQ=255)
                                                                                             tctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCCCTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                                  >  1:466574/1‑127 (MQ=255)
                                                                                             tctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTTGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                           >  2:1051773/1‑150 (MQ=255)
                                                                                             tctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAACCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                                  <  2:466574/127‑1 (MQ=255)
                                                                                               tAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTa                                         >  1:717766/1‑150 (MQ=255)
                                                                                                     ttGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAGGCACGCTGCTGCGGATAACGCTATGTNGAGTCCNCGACAGTGCCTGNCGGAGGCGGACAGCCTGTTTAGCTGTc                                   <  1:1051773/150‑1 (MQ=255)
                                                                                                           ctctACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTc                                                                                  >  1:1203696/4‑97 (MQ=255)
                                                                                                           cTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTc                                                                                  <  2:1203696/97‑1 (MQ=255)
                                                                                                           cTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGccc                             <  2:654079/150‑1 (MQ=255)
                                                                                                               aCCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTgg                         >  1:851947/1‑150 (MQ=255)
                                                                                                                ccGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGc                        <  1:132367/150‑1 (MQ=255)
                                                                                                                       gaaaaaaCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCCACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCt                                      <  1:475465/125‑1 (MQ=255)
                                                                                                                       gaaaaaaCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCCACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCt                                      >  2:475465/5‑129 (MQ=255)
                                                                                                                       aaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATg                 <  1:640503/150‑1 (MQ=255)
                                                                                                                       aaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTGTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATg                 >  1:26574/1‑150 (MQ=255)
                                                                                                                           aaaCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAg             >  2:837610/1‑150 (MQ=255)
                                                                                                                            aaCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGa            <  1:837610/150‑1 (MQ=255)
                                                                                                                             aCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCa                                  >  1:1071794/1‑127 (MQ=255)
                                                                                                                             aCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCa                                  <  2:1071794/127‑1 (MQ=255)
                                                                                                                                cgcgCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGacacc        >  1:170826/1‑149 (MQ=255)
                                                                                                                                     cAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACg   <  2:748971/150‑1 (MQ=255)
                                                                                                                                      aGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACGc  <  1:987034/150‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
ATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACGC  >  NC_002947/1400746‑1401029

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