Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,951,342	G→T	12.5%	intergenic (‑121/+40)	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,342	0	G	T	12.5%	150.1 / 11.9	56	intergenic (‑121/+40)	oguA/PP_2585	8‑oxoguanine deaminase/alpha‑ketoglutaric semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (28/21);  new base T (4/3);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

CAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAG  >  NC_002947/2951203‑2951487                                                                                                                                            |                                                                                                                                                   cAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgcg                                                                                                                                             >  4:102952/1‑147 (MQ=255)   gggTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTAATGATTGTATACAGCTCGAAGTACCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAaga                                                                                                                                         <  8:181558/149‑1 (MQ=39)         gTAGATGACCTGCTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATCGGCGCTCGGTACTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTACCCGTGCCGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCg                                                                                                                                    <  7:233786/148‑1 (MQ=37)            gATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCAc                                                                                                                                >  8:36653/1‑149 (MQ=255)               gACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgc                                                                                                                                              <  6:301288/132‑1 (MQ=255)               gACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCcgc                                                                                                                                              >  5:301288/1‑132 (MQ=255)                  aTGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGtt                                                                                                                           >  7:29044/1‑148 (MQ=255)                      tGCAAGCCTGCGCAATTCGCATCGGAGATTACCTATGTCATAGGCGCTCGGTATTGTAAGAAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGc                                                                                                                      <  8:276224/149‑1 (MQ=255)                       gCAAGCCTGCGCAATACGCATCGGCGAATACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGc                                                                                                                      >  6:159037/1‑148 (MQ=255)                        cAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGcc                                                                                                                     <  5:314560/148‑1 (MQ=255)                         aaGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCAt                                                                                                                   <  7:36653/149‑1 (MQ=255)                           gCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGGATTGTAAGCAAGCTTTTTCCTGATTGTATACAGCTCGACGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                 <  7:230747/149‑1 (MQ=255)                            ccTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGa                                                                                                                 >  5:279659/1‑148 (MQ=255)                                cgcAATACGCATCGGCTTTTACCGTTGTCATAGGCGCTCGTTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAg                                                                                                             <  8:29044/147‑1 (MQ=255)                                cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTg                                                                                                                          >  8:243846/1‑135 (MQ=255)                                cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTg                                                                                                                          <  7:243846/135‑1 (MQ=255)                                cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGa                                                                                                            <  2:281727/148‑1 (MQ=255)                                cgcAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGAt                                                                                                           <  6:279659/149‑1 (MQ=255)                                        gCATCGGCGATTACCGATGTGATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGc                                                                                                   >  5:84453/1‑149 (MQ=255)                                        gCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGCGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGc                                                                                                   >  8:134898/1‑149 (MQ=255)                                                        aTGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTcaca                                                                                   >  4:120038/1‑149 (MQ=255)                                                            cATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCCGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCcccccca                                                                               >  8:112118/1‑148 (MQ=255)                                                             aTAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTc                                                                              >  6:3976/1‑149 (MQ=255)                                                                                   aagcTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGc                                                        >  2:207253/1‑149 (MQ=255)                                                                                     gcTTTATACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTTGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGCGCAATCCCCAGCgg                                                      >  6:219336/1‑149 (MQ=255)                                                                                         tttACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGt                                                    <  3:120038/147‑1 (MQ=255)                                                                                            aCTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGc                                                <  5:159037/148‑1 (MQ=255)                                                                                            aCTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGc                                                <  7:112118/148‑1 (MQ=255)                                                                                             cTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGTGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTcg                                                                                               <  7:220051/100‑1 (MQ=255)                                                                                             cTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGTGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTcg                                                                                               >  8:220051/1‑100 (MQ=255)                                                                                                aTTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCCGCGGGATAGCGtttt                                           >  8:125560/1‑149 (MQ=255)                                                                                                  tGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGttttt                                          >  8:46074/1‑148 (MQ=255)                                                                                                   gTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTc                                         <  7:46074/148‑1 (MQ=255)                                                                                                     ataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                       <  7:64395/148‑1 (MQ=255)                                                                                                     ataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                       >  8:64395/1‑148 (MQ=255)                                                                                                              gAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCgg                              >  2:64194/1‑148 (MQ=255)                                                                                                              gAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgc                             <  5:3976/149‑1 (MQ=255)                                                                                                                  tGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGgcag                           >  3:226143/1‑147 (MQ=255)                                                                                                                     ccGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTg                      >  2:219550/1‑149 (MQ=255)                                                                                                                          ccGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCggg                                                     <  4:250622/113‑1 (MQ=255)                                                                                                                          ccGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCggg                                                     >  3:250622/1‑113 (MQ=255)                                                                                                                          ccGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGt                   >  1:168785/1‑147 (MQ=255)                                                                                                                           cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                      <  4:289701/127‑1 (MQ=255)                                                                                                                           cGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAAc                                      >  3:289701/1‑127 (MQ=255)                                                                                                                           cGGCCTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAg                >  2:158686/1‑149 (MQ=255)                                                                                                                              ccTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCTGGCAATCCCCTGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                       <  5:189997/139‑1 (MQ=38)                                                                                                                              ccTCTTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGCGTTAGCGTTTTTCAACGCATCCGGCAGCAGt                       >  6:189997/1‑139 (MQ=39)                                                                                                                               ctctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg              >  2:259938/1‑147 (MQ=255)                                                                                                                               ctctTCGCGGCTCTAGCCGCCAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg              >  2:137453/1‑147 (MQ=255)                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGAGAAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAggg              >  1:240246/1‑147 (MQ=255)                                                                                                                                tctTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAg           >  1:99592/1‑149 (MQ=255)                                                                                                                                  ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt         >  1:140026/1‑149 (MQ=255)                                                                                                                                  ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt         <  3:102952/149‑1 (MQ=255)                                                                                                                                     gcggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCt       <  4:226143/148‑1 (MQ=255)                                                                                                                                       ggcTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGa    <  7:134898/149‑1 (MQ=255)                                                                                                                                          tAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAg  >  4:163435/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                                   CAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTTCTGGAAG  >  NC_002947/2951203‑2951487

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