Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,087,770 G→C 7.3% R259G (CGC→GGC) ‡ metR‑II ← homocysteine‑binding DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,087,7700GC 7.3% 172.8 / 2.6 54R259G (CGC→GGC) ‡metR‑IIhomocysteine‑binding DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base C (2/2);  total (26/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.14e-01

CGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG  >  NC_002947/3087624‑3087914
                                                                                                                                                  |                                                                                                                                                
cGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCt                                                                                                                                                <  4:43618/149‑1 (MQ=255)
          gCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCtt                                                                                                                                               <  6:29623/140‑1 (MQ=255)
          gCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCtt                                                                                                                                               >  5:29623/1‑140 (MQ=255)
            cACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACa                                                                                                                                     >  8:33857/1‑148 (MQ=255)
            cACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACa                                                                                                                                     >  7:254922/1‑148 (MQ=255)
             aCGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATa                                                                                                                                   <  5:157381/149‑1 (MQ=255)
                 tAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcg                                                                                                                               >  3:136522/1‑149 (MQ=255)
                     acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcgc                                                                                                                            >  2:351634/1‑148 (MQ=255)
                     acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCgagc                                                                                                                                        >  1:22620/1‑133 (MQ=255)
                     acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCgagc                                                                                                                                        <  2:22620/136‑4 (MQ=255)
                           cAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCACCAGCATGTCACTGCGCAGTGCAGCAAAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAAcgg                                                                                                                      >  2:320620/1‑146 (MQ=255)
                              cGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTAtt                                                                                                                  <  5:358193/149‑1 (MQ=255)
                                  gtgAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGt                                                                                                               <  6:300052/148‑1 (MQ=255)
                                       cGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACgg                                                                                                         >  3:148795/1‑149 (MQ=255)
                                             gTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGt                                                                                                   >  2:270981/1‑149 (MQ=255)
                                                cTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTgg                                                                                                <  8:254922/149‑1 (MQ=255)
                                                           gCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGgcagca                                                                                      <  4:136522/148‑1 (MQ=255)
                                                           gCAGAAAATCACGCATGTACGGCACCTCCAGCATGCCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGgcaggcag                                                                                        >  2:106626/1‑146 (MQ=255)
                                                            cAGAAAATCACGCATGTGCGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGa                                                                                    >  5:281279/1‑149 (MQ=255)
                                                            cAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGa                                                                                    >  6:199943/1‑149 (MQ=255)
                                                                aaaTCACGCATGTACGGCACCTCCAGCATGCCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGAcgcc                                                                                <  1:106626/149‑1 (MQ=255)
                                                                 aaTCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCCGCAGAcgccg                                                                               >  7:330932/1‑149 (MQ=255)
                                                                         aTGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACt                                                                                                                       >  8:101794/1‑101 (MQ=255)
                                                                         aTGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCgg                                                                        >  6:364300/1‑148 (MQ=255)
                                                                         aTGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTGGGCGGTTACATAGCCGCGCGAACt                                                                                                                       <  7:101794/101‑1 (MQ=255)
                                                                          tGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa                                                                      <  5:199943/149‑1 (MQ=255)
                                                                          tGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa                                                                      <  5:364300/149‑1 (MQ=255)
                                                                          tGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa                                                                      <  1:351634/149‑1 (MQ=255)
                                                                           gTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGaa                                                                      <  8:13055/148‑1 (MQ=255)
                                                                                   ccTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGt                                                             <  8:283173/149‑1 (MQ=255)
                                                                                              tCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatcatcat                                                  <  4:148795/149‑1 (MQ=255)
                                                                                                     gcAGTGCAGCATAGAGGGTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTc                                           <  3:299087/149‑1 (MQ=255)
                                                                                                        gTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGt                                        >  7:33861/1‑149 (MQ=255)
                                                                                                         tGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCTGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGtt                                       <  6:281279/149‑1 (MQ=255)
                                                                                                                tAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgt                                <  1:270981/149‑1 (MQ=255)
                                                                                                                 ngagGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgt                                >  4:315161/2‑148 (MQ=255)
                                                                                                                     gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAATCAACCa                                                               >  3:51/1‑113 (MQ=255)
                                                                                                                     gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCa                                                               <  4:51/113‑1 (MQ=255)
                                                                                                                      tttGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgtacgtac                          <  7:33857/149‑1 (MQ=255)
                                                                                                                         gCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCACGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCt                        <  6:180165/148‑1 (MQ=255)
                                                                                                                            tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTg                                    >  6:33961/1‑133 (MQ=255)
                                                                                                                            tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTg                                    <  5:33961/133‑1 (MQ=255)
                                                                                                                            tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGc                    >  8:284760/1‑149 (MQ=255)
                                                                                                                             gCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCg                   >  4:7713/1‑149 (MQ=255)
                                                                                                                              cAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGggcagg                                                                                           >  6:41496/1‑76 (MQ=255)
                                                                                                                                 gCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc                                                       <  7:276841/109‑1 (MQ=255)
                                                                                                                                 gCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc                                                       >  8:276841/1‑109 (MQ=255)
                                                                                                                                  ccGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTc              <  3:7713/149‑1 (MQ=255)
                                                                                                                                        ccGCCAAGCGCCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGAc           <  7:50103/145‑1 (MQ=255)
                                                                                                                                        ccGCCAAGCGCCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGAc           >  8:50103/2‑146 (MQ=255)
                                                                                                                                          gCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcgg      >  2:166346/1‑149 (MQ=255)
                                                                                                                                            cAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCt                                     >  4:220763/1‑116 (MQ=255)
                                                                                                                                            cAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCt                                     <  3:220763/116‑1 (MQ=255)
                                                                                                                                              aGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCggg  <  3:315161/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG  >  NC_002947/3087624‑3087914

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: