Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 385,475 A→C 7.8% intergenic (+15/+20) PP_0320 → / ← ltaE methyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947385,4750AC 7.8% 150.5 / 4.4 51intergenic (+15/+20)PP_0320/ltaEmethyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base A (28/19);  new base C (2/2);  total (30/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 = 9.20e-01

TGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCG  >  NC_002947/385336‑385623
                                                                                                                                           |                                                                                                                                                    
tGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCCGGGGCCGCCGTGCGGccc                                                                                                                                             >  1:305243/1‑149 (MQ=255)
  cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggc                                                                                                                                                  >  1:263575/1‑142 (MQ=255)
  cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCAcggc                                                                                                                                                  <  2:263575/142‑1 (MQ=255)
   aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCgg                                                                                                                                          <  4:251882/149‑1 (MQ=255)
    ggggATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCCGGGGCCGCCGTGCGGCCCCTGc                                                                                                                                         <  2:305243/149‑1 (MQ=255)
         tttCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCCGGGGCCGCCGTGCGGCCCCTGcgatct                                                                                                                                    <  8:34435/149‑6 (MQ=255)
           tCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGctct                                                                                                                                      <  4:348317/144‑1 (MQ=255)
           tCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGctct                                                                                                                                      >  3:348317/1‑144 (MQ=255)
             cGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCa                                                                                                                                <  8:380724/149‑1 (MQ=255)
                   ggCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGccac                                                                                                                           >  7:26131/1‑148 (MQ=255)
                    gcgcGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGccacc                                                                                                                          >  6:26564/1‑148 (MQ=255)
                               ggcggcggAGGAGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCg                                                                                                              >  6:340962/1‑149 (MQ=255)
                                 gggcggTGGTGCACGGGGGCCGGGAGGAGTGCCTGGGGTGGGGGGGGTGGGCTGATGACGTGGCGCGGCTGTTGGGTGGGTTCAGGCTCTAGATCGCCGGAGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGa                                                                                                            <  4:14644/148‑1 (MQ=255)
                                     ggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCtctc                                                                                                                                     <  8:64762/120‑1 (MQ=255)
                                     ggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCtctc                                                                                                                                     >  7:64762/1‑120 (MQ=255)
                                          tGCAAGGGTGCCGCGAGGAGTGCCAGGCGTAGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccgccg                                                                                                   >  1:36560/1‑149 (MQ=255)
                                           gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                  >  1:402685/1‑149 (MQ=255)
                                           gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                  >  7:239873/1‑149 (MQ=255)
                                             aaGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATc                                                                                                                   <  5:105734/130‑1 (MQ=255)
                                             aaGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATc                                                                                                                   >  6:105734/1‑130 (MQ=255)
                                             aaGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                  >  6:240339/1‑147 (MQ=255)
                                              aGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAg                                                                                               >  4:245267/1‑149 (MQ=255)
                                                     cgcgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGc                                                                                        <  2:280964/149‑1 (MQ=255)
                                                      gcgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCa                                                                                       >  1:229396/1‑149 (MQ=255)
                                                           gagTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCcgc                                                                                  >  7:291292/1‑149 (MQ=255)
                                                                      gTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCgg                                                                        >  2:143102/1‑148 (MQ=255)
                                                                      gTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGt                                                                       <  4:369520/149‑1 (MQ=255)
                                                                               gTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCAc                                                                >  4:231807/1‑147 (MQ=255)
                                                                                   gCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGca                                                           <  1:143102/148‑1 (MQ=255)
                                                                                      gatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACAt                                                        >  1:203213/1‑148 (MQ=255)
                                                                                       atgaCCTGGCGGGGGAGAGGGGTGGGTTCAGTCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCTCTCAGCCACCAATGATCGTGCGGATATCGGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACAt                                                        <  6:324485/147‑1 (MQ=255)
                                                                                           ccTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcg                                                  >  5:247406/1‑149 (MQ=255)
                                                                                           ccTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcg                                                  >  5:301345/1‑149 (MQ=255)
                                                                                              ggCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                         >  2:340338/1‑123 (MQ=255)
                                                                                              ggCGCGGCAGAGGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                         <  1:340338/123‑1 (MQ=255)
                                                                                                       gATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCt                                      >  2:224767/1‑149 (MQ=255)
                                                                                                          gggtgggtTCAGGCTCTAGATCCCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGGAGCGCGCGCCACCGCTGc                                    <  3:245267/148‑1 (MQ=255)
                                                                                                                gtTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATa                               >  3:402475/1‑147 (MQ=255)
                                                                                                                    aGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGta                         >  6:85347/1‑149 (MQ=255)
                                                                                                                    aGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGt                          <  6:301345/148‑1 (MQ=255)
                                                                                                                     ggCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGtat                        <  2:36560/149‑1 (MQ=255)
                                                                                                                      gCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCAACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGtat                        >  8:254725/1‑148 (MQ=255)
                                                                                                                          tAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCAc                                          >  7:355063/1‑126 (MQ=255)
                                                                                                                          tAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCAc                                          <  8:355063/126‑1 (MQ=255)
                                                                                                                           aGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCAc                                          >  1:52421/1‑125 (MQ=255)
                                                                                                                           aGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCAc                                          <  2:52421/125‑1 (MQ=255)
                                                                                                                             aTCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAgcgc                >  2:116868/1‑149 (MQ=255)
                                                                                                                                  cggggCCGCCGTGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGcc           >  3:132213/2‑149 (MQ=255)
                                                                                                                                  aGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATg                                                       <  6:356019/105‑1 (MQ=255)
                                                                                                                                  aGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATg                                                       >  5:356019/1‑105 (MQ=255)
                                                                                                                                        cGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGcc     <  5:340962/149‑1 (MQ=255)
                                                                                                                                           aCGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCg  >  5:155263/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
TGCAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCG  >  NC_002947/385336‑385623

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: