Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,787,623	G→A	12.2%	intergenic (+70/‑55)	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_002947	2,787,623	0	G	A	12.2%	116.4 / 8.3	51	intergenic (+70/‑55)	ahpC/ahpF	peroxiredoxin/alkylhydroperoxide reductase small subunit/alkyl hydroperoxide reductase subunit F
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base A (3/3);  total (25/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

TGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGC  >  NC_002947/2787483‑2787767                                                                                                                                              |                                                                                                                                                  tgtgCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCAGGCGaggt                                                                                                                                           >  5:219122/1‑146 (MQ=255)   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                                            <  1:262693/149‑1 (MQ=255)                                                                                                tAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCgg                                            <  6:354078/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                                              <  1:170678/122‑1 (MQ=255)                                                                                                                         ctaccGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGc                                              >  2:170678/1‑122 (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                                   >  4:263136/3‑122 (MQ=39)                                                                                                                                    aaaaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAg                                   <  3:263136/120‑1 (MQ=39)                                                                                                                                      aaaCGCCCAGGCGAGGTCGCCCGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTcgacga      >  7:262269/1‑149 (MQ=255)                                                                                                                                          gCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGgc  <  2:88120/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  TGTGCCCGGCCAAGTGGAAAGAAGGCGAAGCCACTCTGGCGCCGTCGCTGGACCTCGTAGGC‑AAGATCTAAGGTCATGCCCATAACGCGGGCGGTAAGCAGCCGCCAAAGCTGAAGTACCTACCGCCCCGTCAAACGCCCGGGCGACCTCGCCTGGGCGTTTTTGTATCCGAGTTTTGAAAAAGGAATCGCCCGTATGTTGGACGCCACGCTTAAATCGCAACT‑GAAAACCTACTTGGAGCGGGTCACCCAGCCGATCGAGATCGTTGCTTCCCTCGACGACGGC  >  NC_002947/2787483‑2787767

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