Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,787,629 C→G 9.4% intergenic (+76/‑49) ahpC → / → ahpF peroxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,787,6290CG 9.4% 156.5 / 7.1 53intergenic (+76/‑49)ahpC/ahpFperoxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base C (23/25);  new base G (2/3);  total (25/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 = 7.58e-01

TGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAA  >  NC_002947/2787485‑2787772
                                                                                                                                                 |                                                                                                                                                
tgCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCATCAAACGCCCGGGCGACCTCg                                                                                                                                              >  7:75374/1‑149 (MQ=255)
 gCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGc                                                                                                                                             >  5:6265/1‑149 (MQ=255)
   ccGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCt                                                                                                                                           <  6:6265/149‑1 (MQ=255)
    cGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCAGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTg                                                                                                                                          >  1:207026/1‑149 (MQ=255)
      gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTggg                                                                                                                                        >  2:236898/1‑149 (MQ=255)
      gCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCggg                                                                                                                                        <  6:219122/149‑1 (MQ=255)
          aGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGAGGTCGCCCGGGCGtt                                                                                                                                    <  4:17776/149‑1 (MQ=255)
            tGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGtttt                                                                                                                                  >  2:262693/1‑149 (MQ=255)
              gAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGCTAAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGttttt                                                                                                                                 >  1:78818/1‑149 (MQ=255)
                     ggCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGa                                                                                                                         <  3:218406/149‑1 (MQ=255)
                      gCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAg                                                                                                                        >  1:340659/1‑149 (MQ=255)
                                  tGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGaaaaa                                                                                                              <  3:355564/147‑1 (MQ=255)
                                                gACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAg                                                                                                             <  4:88578/134‑1 (MQ=255)
                                                gACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAg                                                                                                             >  3:88578/1‑134 (MQ=255)
                                                gACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGt                                                                                              >  7:33352/1‑149 (MQ=255)
                                                 aCCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCATCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGtt                                                                                             <  8:75374/149‑1 (MQ=255)
                                                   cTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTg                                                                                            >  4:46511/1‑148 (MQ=255)
                                                           c‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACg                                                                                   >  2:232170/1‑149 (MQ=255)
                                                             aaGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGTCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGc                                                                                  <  2:207026/149‑1 (MQ=255)
                                                                  cTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTaaa                                                                             >  8:236394/1‑149 (MQ=255)
                                                                     aGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                          <  7:236394/149‑1 (MQ=255)
                                                                      ggTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCg                                                                          >  5:195785/1‑148 (MQ=255)
                                                                        tCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCaa                                                                       >  5:354078/1‑149 (MQ=255)
                                                                            gCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑g                                                                   >  4:3971/1‑148 (MQ=255)
                                                                              ccATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑Gaa                                                                 >  1:88120/1‑148 (MQ=255)
                                                                                    cgcgGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTAc                                                          <  8:33352/149‑1 (MQ=255)
                                                                                     gcgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACt                                                         >  7:178545/1‑149 (MQ=255)
                                                                                      cgggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACtt                                                        <  2:78818/149‑1 (MQ=255)
                                                                                       gggcggTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTg                                                       <  3:46511/149‑1 (MQ=255)
                                                                                             tAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCgg                                                 <  6:354078/149‑1 (MQ=255)
                                                                                             tAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCgg                                                 <  1:262693/149‑1 (MQ=255)
                                                                                               agcagcCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGt                                               >  1:185797/1‑149 (MQ=255)
                                                                                                       ccAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGc                                       <  1:23384/149‑1 (MQ=255)
                                                                                                            gCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                  <  3:3971/149‑1 (MQ=255)
                                                                                                              tGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCga                                <  6:195785/149‑1 (MQ=255)
                                                                                                               gAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATc                                  >  3:23393/1‑146 (MQ=255)
                                                                                                                aaGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACTCGAAAACCTACTTGGAGCGGGTCACCCAGCCGATCg                                 <  1:236898/147‑1 (MQ=255)
                                                                                                                     cctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGtt                         <  1:232170/149‑1 (MQ=255)
                                                                                                                      ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGc                                                   >  2:170678/1‑122 (MQ=255)
                                                                                                                      ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGc                                                   <  1:170678/122‑1 (MQ=255)
                                                                                                                        accGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCGCGTATGTTGGACGCCACGCTTAATTCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCt                      >  7:157947/1‑149 (MQ=255)
                                                                                                                          cGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATc                            <  7:319678/141‑1 (MQ=255)
                                                                                                                          cGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTAGGAGCGGGTCACCCAGCCGATCGAGATc                            >  8:319678/1‑141 (MQ=255)
                                                                                                                            ccccGACAAACCCCCCCCCGACCTCCCCAGCGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTccc                  <  2:140897/149‑1 (MQ=255)
                                                                                                                             cccGTCAACCGGCCGGTCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCt                 <  8:178545/149‑1 (MQ=255)
                                                                                                                              gcGTCAAACGGCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTAGTTGGAGCGGGTCAGCCAGCCGATCGAGATTGTTGCTTCCCtg                <  2:340659/148‑2 (MQ=255)
                                                                                                                                 aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAg                                        <  3:263136/120‑1 (MQ=39)
                                                                                                                                 aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAg                                        >  4:263136/3‑122 (MQ=39)
                                                                                                                                   aaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgacga           >  7:262269/1‑149 (MQ=255)
                                                                                                                                       gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgc       <  2:88120/149‑1 (MQ=255)
                                                                                                                                            ggCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa  <  4:23393/149‑1 (MQ=255)
                                                                                                                                                aCCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa  >  7:82965/1‑145 (MQ=255)
                                                                                                                                                aCCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGaa  <  8:82965/145‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGCGCGAA  >  NC_002947/2787485‑2787772

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