Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,951,334	G→C	23.0%	intergenic (‑113/+48)	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,334	0	G	C	23.0%	74.5 / 17.9	39	intergenic (‑113/+48)	oguA/PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (14/16);  new base C (6/3);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTT  >  NC_002947/2951186‑2951480                                                                                                                                                     |                                                                                                                                                    gACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTg                                                                                                                                                     >  8:120/1‑149 (MQ=255)   cGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTgg                                                                                                                                                    <  6:2534/148‑1 (MQ=255)       gagCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCt                                                                                                                                               >  6:4398/1‑149 (MQ=255)            tttttGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGc                                                                                                                                          >  8:22911/1‑149 (MQ=255)                tGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgc                                                                                                                                       <  3:98189/148‑1 (MQ=255)                  cAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTATTCGCGGCTCTAGCcgcca                                                                                                                                     >  2:306967/1‑146 (MQ=255)                   aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgcg                                                                                                                                      >  7:21609/1‑146 (MQ=255)                   aGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAag                                                                                                                                   <  4:139232/149‑1 (MQ=255)                    ggggTTGTAGATGACATGGTGCAAGCATGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                  <  6:268741/149‑1 (MQ=39)                    gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAaga                                                                                                                                  <  4:265924/149‑1 (MQ=255)                       tttGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGc                                                                                                                               <  5:4398/149‑1 (MQ=255)                         tGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGcc                                                                                                                              <  1:306967/148‑1 (MQ=255)                                  cATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGtt                                                                                                                    >  5:79853/1‑149 (MQ=255)                                   aTGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTg                                                                                                                   <  6:79853/149‑1 (MQ=255)                                   aTGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTg                                                                                                                   >  2:45775/1‑149 (MQ=255)                                      gTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACCGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑Ca                                                                                                                <  7:120/148‑1 (MQ=255)                                              cTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGaaa                                                                                                       <  7:22911/149‑1 (MQ=255)                                                 cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGa                                                                                                     <  1:45775/148‑1 (MQ=255)                                                               gCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTcgcg                                                                                      <  1:200976/149‑1 (MQ=255)                                                               gCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTcgcg                                                                                      >  1:80610/1‑149 (MQ=255)                                                                 gATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGgc                                                                                    <  6:190382/149‑1 (MQ=255)                                                                          tGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCacac                                                                           >  2:48728/1‑149 (MQ=255)                                                                            tCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACAcc                                                                          >  5:40515/1‑148 (MQ=255)                                                                              aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑c                                                                                                                 >  5:72179/1‑107 (MQ=255)                                                                              aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑c                                                                                                                 <  6:72179/107‑1 (MQ=255)                                                                                    gcTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCaa                                                                 >  2:93952/1‑149 (MQ=255)                                                                                               tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATc                                                       <  1:93952/148‑1 (MQ=255)                                                                                               tAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATc                                                       <  6:40515/148‑1 (MQ=255)                                                                                                aagcaagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGACAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATc                                                       <  1:253718/148‑1 (MQ=255)                                                                                                        tttttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGACGCGAAGAGGCAGGCACGGGTTG‑CAGCCATGAAAGATAACAACGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                             >  7:22257/1‑149 (MQ=255)                                                                                                                 aTTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGtttt                                    >  2:308255/1‑149 (MQ=255)                                                                                                                     tataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                >  8:175336/1‑149 (MQ=255)                                                                                                                       taCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                               <  8:21609/148‑1 (MQ=255)                                                                                                                            cTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATcc                         <  7:175336/149‑1 (MQ=255)                                                                                                                               gAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                       >  5:176977/1‑148 (MQ=255)                                                                                                                                aaGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                       >  3:107571/1‑147 (MQ=255)                                                                                                                                           ccGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGTCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCTGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTc           <  8:22257/148‑1 (MQ=255)                                                                                                                                              gCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg       >  5:122186/1‑149 (MQ=255)                                                                                                                                              gCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTGTTCAACGCATCCGGCAGCAGTGCGTCAggg       >  3:197903/1‑149 (MQ=39)                                                                                                                                                ctctTCGCGGCTCTAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg       >  5:150992/1‑147 (MQ=255)                                                                                                                                                ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCACCGCATCCGGCAGCAGTGCGTCAggg       >  8:180241/1‑147 (MQ=39)                                                                                                                                                     ggcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt  <  6:176977/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GACGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTG‑CAGCCATGAAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTT  >  NC_002947/2951186‑2951480

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑