Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,087,768 G→C 7.5% R259R (CGC→CGG)  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,7680GC 7.5% 162.1 / 2.1 51R259R (CGC→CGG) metR‑IIhomocysteine‑binding DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base C (2/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

AACCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG  >  NC_002947/3087621‑3087914
                                                                                                                                                   |                                                                                                                                                  
aaCCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGgc                                                                                                                                                   <  4:251616/149‑1 (MQ=255)
         tGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCga                                                                                                                                          <  2:339459/149‑2 (MQ=255)
                aCGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATa                                                                                                                                   <  7:248726/149‑1 (MQ=255)
                        acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcgc                                                                                                                            >  2:291251/1‑148 (MQ=255)
                        acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGAATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcgc                                                                                                                            >  3:364707/1‑148 (MQ=255)
                           tGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCtt                                                                                                                                               <  4:308584/126‑1 (MQ=255)
                           tGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCtt                                                                                                                                               >  3:308584/1‑126 (MQ=255)
                             ccAGCTTGGTGAACGAGATGTCCTTGGCCGTGAGAATAAAATCACGCATGTTCGGCACCTCCAGCATGTCACTGCGCAGTTCAGCATAGAGGTTTGCCTTCAGGCCGTGCTAGCCAAGGCGCTTGGCTGtt                                                                                                                                        <  4:35756/131‑1 (MQ=255)
                             ccAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACAGCGCAGTGCAGCATAGAGGTTTGCCTCCAGGCCGTGCACGCCAAGCCGCTTGGCGGtt                                                                                                                                        >  3:35756/1‑131 (MQ=255)
                                      tgAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTg                                                                                                              >  8:200920/1‑148 (MQ=255)
                                         aCGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACg                                                                                                          >  3:17078/1‑149 (MQ=255)
                                         aCGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACg                                                                                                          >  7:158313/1‑149 (MQ=255)
                                                   cTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCAc                                                                                                           <  6:120550/138‑1 (MQ=255)
                                                   cTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCAc                                                                                                           >  5:120550/1‑138 (MQ=255)
                                                       gCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACt                                                                                                                       <  2:47774/122‑1 (MQ=255)
                                                       gCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACt                                                                                                                       >  1:47774/1‑122 (MQ=255)
                                                            gAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGgcaggcag                                                                                        >  7:239300/1‑148 (MQ=255)
                                                             aGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGgcagc                                                                                       <  4:17078/148‑1 (MQ=255)
                                                             aGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGgcagc                                                                                       <  8:158313/148‑1 (MQ=255)
                                                                        aCGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGgcaggcag                                                                                        >  6:95062/1‑136 (MQ=255)
                                                                        aCGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGgcaggcag                                                                                        <  5:95062/136‑1 (MQ=255)
                                                                        aCGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGcc                                                                           >  6:349698/1‑149 (MQ=255)
                                                                         cGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGccc                                                                          >  6:77876/1‑149 (MQ=255)
                                                                          gCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCg                                                                         <  5:349698/149‑1 (MQ=255)
                                                                              gTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGcagcag                                                                                     <  4:118673/133‑1 (MQ=255)
                                                                              gTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGcagcag                                                                                     >  3:118673/1‑133 (MQ=255)
                                                                                aCGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGc                                                                    >  5:258075/1‑148 (MQ=255)
                                                                                       cTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGtt                                                                                                                                        >  8:138608/1‑73 (MQ=255)
                                                                                       cTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGtt                                                                                                                                        <  7:138608/73‑1 (MQ=255)
                                                                                          cAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGc                                                          <  5:77876/148‑1 (MQ=255)
                                                                                          acgCATGTCACTGCGCAGTGCAGCATATAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGc                                                          <  1:136716/146‑1 (MQ=255)
                                                                                            gCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc                                                       <  7:200920/149‑1 (MQ=255)
                                                                                                   aCTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCatcatcatc                                                 <  6:258075/148‑1 (MQ=255)
                                                                                                    cTGCGCAGTGCAGCATAGAGGTTTGCCTGTAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCAc                                               >  3:43658/1‑149 (MQ=255)
                                                                                                               agcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTAAATAGCCGCGCGAACTGTATTCTTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTg                                    <  5:87238/149‑1 (MQ=255)
                                                                                                                  aTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAg                                 >  7:249824/1‑149 (MQ=255)
                                                                                                                    agagGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGc                                                          <  6:310561/122‑1 (MQ=255)
                                                                                                                    agagGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGc                                                          >  5:310561/1‑122 (MQ=255)
                                                                                                                    agagGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgta                               >  8:316721/1‑149 (MQ=255)
                                                                                                                        gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc                                                       <  7:233281/121‑1 (MQ=255)
                                                                                                                        gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatc                                                       >  8:233281/1‑121 (MQ=255)
                                                                                                                              cttcAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCtgtg                     <  4:364707/146‑1 (MQ=255)
                                                                                                                               tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGggcagg                                                                                           >  6:361325/1‑78 (MQ=255)
                                                                                                                               tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGggcagg                                                                                           <  5:361325/78‑1 (MQ=255)
                                                                                                                               tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGaa                                  <  7:304264/135‑1 (MQ=255)
                                                                                                                               tGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGaa                                  >  8:304264/1‑135 (MQ=255)
                                                                                                                                   ggCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGc                >  2:271989/1‑149 (MQ=255)
                                                                                                                                    gCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCt               >  6:128783/1‑149 (MQ=255)
                                                                                                                                        tGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGAc           <  7:92534/149‑1 (MQ=255)
                                                                                                                                           ccGCCAAGCGCCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGca                                                         <  2:362109/99‑1 (MQ=255)
                                                                                                                                           ccGCCAAGCGCCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGca                                                         >  1:362109/2‑100 (MQ=255)
                                                                                                                                                 aGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCggg  <  8:239300/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG  >  NC_002947/3087621‑3087914

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: ‑