Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,400,880 A→T 8.2% intergenic (+37/‑132) 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,8800AT 8.2% 133.3 / 4.6 49intergenic (+37/‑132)ybgF/queETolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base T (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

AGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACGC  >  NC_002947/1400731‑1401029
                                                                                                                                                     |                                                                                                                                                     
aGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCa                                                                                                                                                       <  1:550219/150‑1 (MQ=255)
               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                                  <  2:1071794/127‑1 (MQ=255)
                                                                                                                                            aCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCa                                  >  1:1071794/1‑127 (MQ=255)
                                                                                                                                               cgcgCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGacacc        >  1:170826/1‑149 (MQ=255)
                                                                                                                                                    cAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACg   <  2:748971/150‑1 (MQ=255)
                                                                                                                                                     aGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACGc  <  1:987034/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTACGC  >  NC_002947/1400731‑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: ‑