Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,951,349	G→C	13.5%	intergenic (‑128/+33)	oguA ← / ← PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,951,349	0	G	C	13.5%	136.7 / 12.0	51	intergenic (‑128/+33)	oguA/PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base C (3/4);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

AGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCGGGC  >  NC_002947/2951204‑2951497                                                                                                                                                  |                                                                                                                                                      aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgcg                                                                                                                                                       >  4:93926/1‑146 (MQ=255)  ggggTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                                   <  4:1455/149‑1 (MQ=255)  gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTTTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                                   <  8:197118/149‑1 (MQ=39)   ggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                                  <  4:193820/149‑1 (MQ=255)    gTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAagag                                                                                                                                                  <  8:114135/148‑1 (MQ=255)     tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGc                                                                                                                                                >  4:212580/1‑149 (MQ=255)            aTGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACg                                                                                                                                         >  2:220949/1‑149 (MQ=255)            aTGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                                   >  5:100090/1‑139 (MQ=255)            aTGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                                   <  6:100090/139‑1 (MQ=255)             tGAGATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCGCGGTATTGTAAGCAAGCTCTTTACTGATTGTATACAGCTCGAAGTGCCCGGGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACgg                                                                                                                                        <  1:216954/149‑1 (MQ=255)                            cTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                           >  1:70039/1‑147 (MQ=255)                             tGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                           >  6:47805/1‑146 (MQ=255)                               cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGaaa                                                                                                                         <  5:47805/146‑1 (MQ=255)                                                   aCCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGgctgc                                                                                                  <  6:107911/148‑1 (MQ=255)                                                    ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                           >  3:158152/1‑123 (MQ=255)                                                    ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                           <  4:158152/123‑1 (MQ=255)                                                    ccGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGgctgct                                                                                                 >  6:240329/1‑148 (MQ=255)                                                     cGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTca                                                                                               >  3:249012/1‑149 (MQ=255)                                                             tAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTc                                                                                        >  7:266299/1‑148 (MQ=255)                                                               gtagcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGATGTGGCCGTGCCGTCCTCTTGGGGGCTAGAGGCGCGAAGAGGCCGGCACGGGTTGCATCCATGAAAGATCACACCGGCT‑CGCGGCTGCTCACACCCTcct                                                                                      <  4:58169/145‑2 (MQ=255)                                                               ggCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCa                                                                                      <  8:266299/148‑1 (MQ=255)                                                                 cgcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCa                                                                                   <  2:70039/149‑1 (MQ=255)                                                                  gcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCaa                                                                                  >  2:82448/1‑149 (MQ=255)                                                                     cGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCg                                                                               >  3:120689/1‑149 (MQ=255)                                                                          ttGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCGTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTGACCCCCTCCAGCACCCGGGCa                                                                           >  8:102397/1‑148 (MQ=255)                                                                                  aagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                                  >  5:251600/1‑149 (MQ=255)                                                                                  aagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                                  >  6:311735/1‑149 (MQ=255)                                                                                     cTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCggg                                                               <  3:93926/149‑1 (MQ=255)                                                                                        tttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCTCCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCggg                                                               <  1:220949/147‑1 (MQ=255)                                                                                          tACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTa                                                            <  3:155879/147‑1 (MQ=255)                                                                                          tACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTa                                                            >  4:155879/1‑147 (MQ=255)                                                                                          tACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGc                                                          >  6:118598/1‑149 (MQ=255)                                                                                             tGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGtt                                                       >  8:174373/1‑149 (MQ=255)                                                                                                   tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                                 <  5:240329/149‑1 (MQ=255)                                                                                                     taCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                                <  3:63024/148‑1 (MQ=255)                                                                                                     taCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                                <  4:120689/148‑1 (MQ=255)                                                                                                                gTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                                     >  8:184344/1‑148 (MQ=255)                                                                                                                   cccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                                 >  3:214361/1‑149 (MQ=255)                                                                                                                    ccGTGCCGGCCTCATGGCGGCTAGAGCCGCGAAGAGGCCGGCACCGGTTGCCGCCATGAAAGATCACACCGCCT‑CGCCGCTGCTCACACCCACCAGCAACCCGGCAATCCCCACCGGGTTAGCCTTTTTCAACGCAACCGGCAGCAGTg                                >  6:46158/1‑149 (MQ=255)                                                                                                                           ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                         <  5:118598/149‑1 (MQ=255)                                                                                                                           ggCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAgg                         <  7:184344/149‑1 (MQ=255)                                                                                                                               tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCACCCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg                     >  2:97882/1‑149 (MQ=255)                                                                                                                               tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                                     >  5:102878/1‑133 (MQ=255)                                                                                                                               tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                                     <  6:102878/133‑1 (MQ=255)                                                                                                                                 ttcgcggcTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCACCCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt                   >  7:259236/4‑149 (MQ=255)                                                                                                                                    gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCt                 <  6:251600/148‑1 (MQ=255)                                                                                                                                    gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCt                 <  3:212580/148‑1 (MQ=255)                                                                                                                                      ggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTgg               >  6:209858/1‑148 (MQ=255)                                                                                                                                       gcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGaa             <  5:209858/149‑1 (MQ=255)                                                                                                                                       gcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGaa             <  5:311735/149‑1 (MQ=255)                                                                                                                                           gagCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGca          >  4:2274/1‑148 (MQ=255)                                                                                                                                            agCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGcaca        >  6:215127/1‑149 (MQ=255)                                                                                                                                                gcgAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCgg    <  7:174373/149‑1 (MQ=255)                                                                                                                                                  gAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCCGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCGGgc  >  8:227463/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      AGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAAAGATCACACCGCCT‑CGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAGCACACCGGGC  >  NC_002947/2951204‑2951497

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