Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,455,255 T→A 7.9% intergenic (‑137/‑31) mfd ← / → gapB transcription‑repair coupling factor/glyceraldehyde‑3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,455,2550TA 7.9% 135.9 / 2.8 51intergenic (‑137/‑31)mfd/gapBtranscription‑repair coupling factor/glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (21/26);  new base A (2/2);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

AGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCG  >  NC_002947/2455108‑2455392
                                                                                                                                                   |                                                                                                                                         
aGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTTGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTc                                                                                                                                          <  2:208744/149‑1 (MQ=255)
    cTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATaa                                                                                                                                      >  5:63198/1‑149 (MQ=255)
                  tGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCttttt                                                                                                                        >  2:332763/1‑149 (MQ=255)
                   gCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCTACATTATGACGccc                                                                                                                              >  1:219044/1‑142 (MQ=255)
                   gCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCTACATTATGACGccc                                                                                                                              <  2:219044/142‑1 (MQ=255)
                    cACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCtttttt                                                                                                                       <  1:332763/148‑1 (MQ=255)
                    cACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTc                                                                                                                      <  1:180500/149‑1 (MQ=255)
                       tCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCTACATTATGACGccccaac                                                                                                                          <  4:197569/142‑4 (MQ=255)
                       tCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCTACATTATGACGccccaac                                                                                                                          >  3:197569/1‑139 (MQ=255)
                       tCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGt                                                                                                                   >  4:241377/1‑149 (MQ=255)
                        gcAAGCGTTGCGGCATAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATTAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCt                                                                                                                     <  6:117404/145‑1 (MQ=255)
                        ccAAGAGATGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGAGTCATAATGTAGCCCCTTTATTCt                                                                                                                     >  5:117404/1‑146 (MQ=255)
                         cAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGt                                                                                                                   <  6:314397/147‑1 (MQ=255)
                          aaGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTccc                                                                                                                <  6:63198/149‑1 (MQ=255)
                          aaGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTccc                                                                                                                >  5:232708/1‑149 (MQ=255)
                              gTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACa                                                                                                            <  3:371500/149‑1 (MQ=255)
                                 gCGGCAACGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGACAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGGTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACAtgt                                                                                                         <  2:222040/149‑1 (MQ=255)
                                      aaaGTGACCGATTGTACCGGTGACAGTGCCGGCTATCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAg                                                                                                    >  1:99096/1‑149 (MQ=255)
                                        agagACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGt                                                                                                  >  3:229404/1‑149 (MQ=255)
                                              cGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTAACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTTTCCCTACATGTGGAAGGTTCACGt                                                                                            >  7:266685/1‑149 (MQ=255)
                                                    tACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCa                                                                                      <  1:301349/149‑1 (MQ=255)
                                                      ccGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAg                                                                                     >  3:203399/1‑148 (MQ=255)
                                                                tGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGAc                                                                            >  6:336899/1‑147 (MQ=255)
                                                                tGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCa                                                                          <  3:241377/149‑1 (MQ=255)
                                                                         tCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCtt                                                                   <  5:189131/147‑1 (MQ=255)
                                                                           aGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGcc                                                                                <  6:356876/132‑1 (MQ=255)
                                                                           aGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGcc                                                                                >  5:356876/1‑132 (MQ=255)
                                                                             gCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAg                                                             >  6:308565/1‑149 (MQ=255)
                                                                              cTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGt                                                            >  3:196384/1‑149 (MQ=255)
                                                                               tCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTg                                                           >  7:252644/1‑149 (MQ=255)
                                                                                  aaaGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGAt                                                        <  4:196384/149‑1 (MQ=255)
                                                                                      ggCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGAt                                                    <  6:232708/149‑1 (MQ=255)
                                                                                       gCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATc                                                   <  4:170275/149‑1 (MQ=255)
                                                                                       gCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATc                                                   >  3:170275/1‑149 (MQ=255)
                                                                                          tGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTg                                                <  2:99096/149‑1 (MQ=255)
                                                                                            aaCACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGaa                                              <  5:336899/149‑1 (MQ=255)
                                                                                              cacTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTc                                                                                                <  5:276728/97‑1 (MQ=255)
                                                                                              cacTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTc                                                                                                >  6:276728/1‑97 (MQ=255)
                                                                                               acTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGcc                                           >  4:345830/1‑149 (MQ=255)
                                                                                                        cGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGaa                                  <  1:215293/149‑1 (MQ=255)
                                                                                                        cGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGaa                                  <  5:308565/149‑1 (MQ=255)
                                                                                                                 gACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATc                         <  4:229404/149‑1 (MQ=255)
                                                                                                                      gCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGccc                                          <  8:162540/127‑1 (MQ=255)
                                                                                                                         cTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTAt                 <  3:345830/149‑1 (MQ=255)
                                                                                                                          tCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCg                                                  <  5:319969/115‑1 (MQ=255)
                                                                                                                          tCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCg                                                  >  6:319969/1‑115 (MQ=255)
                                                                                                                            aTTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCgg              >  6:61625/1‑149 (MQ=255)
                                                                                                                               ggcGGTTCCGCGTTGGTGCTTCATAATGTAGCGCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCAGGACCAGTGTGTTGGTGAGTGGATTTATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTc            <  8:266685/146‑1 (MQ=255)
                                                                                                                                        gcgTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCg  >  5:169224/1‑149 (MQ=255)
                                                                                                                                          gTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGccc                                                                               >  4:113044/1‑70 (MQ=255)
                                                                                                                                          gTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGccc                                                                               <  3:113044/70‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
AGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCG  >  NC_002947/2455108‑2455392

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