Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,455,254 G→T 9.2% intergenic (‑136/‑32) 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,2540GT 9.2% 153.0 / 4.8 51intergenic (‑136/‑32)mfd/gapBtranscription‑repair coupling factor/glyceraldehyde‑3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base T (2/3);  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 = 9.46e-01

GCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCG  >  NC_002947/2455106‑2455392
                                                                                                                                                    |                                                                                                                                          
gcAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCg                                                                                                                                            >  2:215293/1‑149 (MQ=255)
  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                                                                                                                >  5:232708/1‑149 (MQ=255)
                            aaGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTccc                                                                                                                <  6:63198/149‑1 (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                                                   >  3:170275/1‑149 (MQ=255)
                                                                                         gCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATc                                                   <  4:170275/149‑1 (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)
                                                                                                                                                    |                                                                                                                                          
GCAGAACTGACACGGGAGATTGCACTCCAAGCGTTGCGGCAAAGAGACCGATTGTACCGGTGACAGTGCCGGCTGTCAGGCTCGAAAGGGCATGAACACTGGCTTGCGTCACTGCGACAGGCACTCATTGCCGGTTCCGCGTTGGGGCGTCATAATGTAGCCCCTTTTTTCTGTCCCTACATGTGGAAGGTTCCCGTGACTCAGAAGCCCGACCAGTGTCTTGGTGAGTGGATTGATCGTGAAGCCCTGGCTGAAGCGATGATCCCGCTTATCGGTCAGCTCTACCG  >  NC_002947/2455106‑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: ‑