Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 262,378 T→G 15.8% A77A (GCA→GCC)  U712_01225 ← Homocysteine S‑methyltransferase ybgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881262,3780TG15.8% 123.4 / 4.0 58A77A (GCA→GCC) U712_01225Homocysteine S‑methyltransferase ybgG
Reads supporting (aligned to +/- strand):  ref base T (13/35);  new base G (8/2);  total (21/37)
Fisher's exact test for biased strand distribution p-value = 2.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTATAGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGCTCTTCCATCAGGATTTTTGCCGACCACA  >  CP006881/262241‑262505
                                                                                                                                         |                                                                                                                               
ccTATAGAAGCCGCAATGATCGGTTTTGGCCGAGTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGCCAGCTCGAGGCGCCGGCGTGCTTCCGCGTCGCGCCGCCCGCGTg                                                                                                                                >  8:307244/1‑139 (MQ=255)
  tataGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGCCCAAAATTCAGCACGTGCTTCAGCTGCGATGCTGACGGACACCTCGATGAGCCGGCGTGGTTCCGCATCGCTCCGCCCGCGGGCg                                                                                                                              >  8:345587/1‑139 (MQ=255)
  tataGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACTAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCCCTTAGCTAAGCCCGCGTGCg                                                                                                                              >  8:101017/1‑139 (MQ=255)
  tataGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGGTCCGCTTCGCGCCCCCGGCGTGCg                                                                                                                              >  2:257558/1‑139 (MQ=255)
  tataGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCGTCCCCGTCGGGCCGCCCCCGGGCg                                                                                                                              >  5:89434/1‑139 (MQ=255)
     aGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGGTCCGCGTGGCGCCGCCCGCGGGCg                                                                                                                              >  3:263686/1‑136 (MQ=255)
     aGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCa                                                                                                                           <  2:411600/139‑1 (MQ=255)
      gacgCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCg                                                                                                                              <  4:326714/132‑1 (MQ=255)
      gAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCg                                                                                                                              <  4:410348/135‑1 (MQ=255)
        aGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACTGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCa                                                                                                                           <  3:186828/136‑1 (MQ=255)
         gCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGTTCCGCGTCGCTCCGCCCGCGTGCGGCAAAGc                                                                                                                       >  8:406003/1‑139 (MQ=255)
         gCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGTTCCGCGTAGCTCAGCCCGCCTGCGGCCAAGc                                                                                                                       >  6:419848/1‑139 (MQ=255)
         gCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGGTCCCCTTGGCGGCGCCCGCGGGGGGCacacg                                                                                                                       >  7:361572/1‑135 (MQ=255)
         gCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCGCAGCCCGCGTGCGGCAAAGc                                                                                                                       >  1:283499/1‑139 (MQ=255)
         gCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCGTCCGCTTGGCTCCGCCCCCGGGCGGCCAAGc                                                                                                                       >  7:99379/1‑139 (MQ=255)
         gCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCGTCCGCGTCGCGCACCCCGCCTGCGGCCAAGc                                                                                                                       >  8:214444/1‑139 (MQ=255)
         gCCGCAATCAACTTTTATGGCCGATTCAAACGGTTCTCTTCAACCGACCAAAATACATCACGTGCTTCAGCTCAGATCCTGACGGACAACTCCATGAGCCCGCGTACTTACCCATAGCCCAGCCCCCCTGGGGCaatcc                                                                                                                       >  6:310161/1‑136 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGGTTCGCTGCGCGCCGCCCCCCGGCGGCCAAGc                                                                                                                       >  3:428909/1‑136 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGc                                                                                                                       <  3:343140/136‑1 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGc                                                                                                                       <  3:41750/136‑1 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCtt                                                                                                                    <  6:391949/139‑1 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCtt                                                                                                                    <  1:288517/139‑1 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCGGCGATGCTGAGGGCCAACTCGATGAGCCGGCGTGCTTCCGCGTCCCGCCGCCCCCGTGCGGCcaacc                                                                                                                       >  3:368088/1‑134 (MQ=255)
            gCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCGTCCGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCtt                                                                                                                    <  1:321858/139‑1 (MQ=255)
             cAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGGTGCCCCTGGCGCCGCCCCCGGGCGGCGAAGcc                                                                                                                      >  3:156324/1‑136 (MQ=255)
                tGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCCCAGCCCGCGTGCGGCAAAGCCtt                                                                                                                    >  4:67729/1‑135 (MQ=255)
                              cGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGc                                                                                                  <  1:222569/139‑1 (MQ=255)
                                   cAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGc                                                                                              <  2:283499/138‑1 (MQ=255)
                                     aaCGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGt                                                                                           <  2:43502/139‑1 (MQ=255)
                                       cGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCgg                                                                                            <  3:47853/136‑1 (MQ=255)
                                          ttCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTg                                                                                      <  1:86555/139‑1 (MQ=255)
                                            ctctTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGca                                                                                    <  5:254455/139‑1 (MQ=255)
                                                    gCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGGCTCCGCGTCGCACAGCCCCCCGGCGGCAAAGCCTTGATATCTATACTGATAGCTGGCGGTGCTTTCCCaaactgc                                                                            >  7:347184/1‑134 (MQ=255)
                                                           aaaTTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCt                                                                        <  3:363957/136‑1 (MQ=255)
                                                              ttCATCACGTGCTTCAGCTGCGCTGCTGACGGACAACTAGATGAGCCGGCGTGCGTCCGCGTCGCTCAGCCCGCGGGAGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGcc                                                                     <  3:66806/136‑1 (MQ=255)
                                                              ttCATCACGTGCTTCAGCCGCGATGCTGACGGACAACTCGATGAGCCGGCGAGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGAAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGcc                                                                     <  3:72735/136‑1 (MQ=255)
                                                                   caCGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTACCCTTGGCGCAACCCGCGTTCGGCAAAGCCATAAAAAGTAGACTGATAGCAGGCCGGGATTGCACAAACGGCGCCTGCCGCAAaat                                                              >  5:161978/1‑138 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaataat                                                           <  1:79167/139‑1 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaataat                                                           <  5:409831/139‑1 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaataat                                                           <  5:419848/139‑1 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaataat                                                           <  1:239606/139‑1 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaataat                                                           <  8:417766/139‑1 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaat                                                              <  3:217418/136‑1 (MQ=255)
                                                                     cGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTGCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAaataat                                                           <  1:101167/139‑1 (MQ=255)
                                                                     acttcTTCAGCAGCGAAGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCAGAATCGGCGCCTGCCGCAAaat                                                              <  3:275195/132‑1 (MQ=255)
                                                                         cTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGt                                                       <  5:102971/139‑1 (MQ=255)
                                                                             aGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGt                                                       <  4:156324/135‑1 (MQ=255)
                                                                              gCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTAt                                                     <  3:67729/136‑1 (MQ=255)
                                                                              gCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGaa                                                  <  7:137126/139‑1 (MQ=255)
                                                                              gCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGaa                                                  <  5:364857/139‑1 (MQ=255)
                                                                              gCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGaa                                                  <  1:339241/139‑1 (MQ=255)
                                                                                   gATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACt                                                <  4:58570/136‑1 (MQ=255)
                                                                                          aCGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAAt                                        <  7:300156/137‑1 (MQ=255)
                                                                                                     gATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGctc                            <  7:214444/138‑1 (MQ=255)
                                                                                                                         cTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGCTCTTCCATCAGGATTTTTGCCGa       <  1:273613/139‑1 (MQ=255)
                                                                                                                         cTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAAAGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGCTCTTCCATCAGGATTTTTTCCGa       >  6:346618/1‑139 (MQ=255)
                                                                                                                              cTCAGCCCGCGGGCGGCAAAGCCTTCAACAGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGCTCTTCCATCAGGATTTTTGCCGACcaca  <  6:52383/139‑1 (MQ=255)
                                                                                                                               tCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGCTCTTCCATCAGGATTTTTGCCGACca    >  4:408595/1‑136 (MQ=255)
                                                                                                                                         |                                                                                                                               
CCTATAGAAGCCGCAATGATCGGTTTTGGCCGATTCAAACGGTTCTCTTCAAGCGACCAAAATTCATCACGTGCTTCAGCTGCGATGCTGACGGACAACTCGATGAGCCGGCGTGCTTCCGCTTCGCTCAGCCCGCGTGCGGCAAAGCCTTCAAATGTAGACTGATAGCTGGCGGTGATTGCACAATCGGCGCCTGCCGCAAAATAATCAGTATGAACTTGTTTAATCAGCTCCGGCTCTTCCATCAGGATTTTTGCCGACCACA  >  CP006881/262241‑262505

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