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

CGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTT  >  NC_002947/2951188‑2951480                                                                                                                                                   |                                                                                                                                                    cGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTggc                                                                                                                                                   <  5:247874/149‑1 (MQ=255) cGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTggc                                                                                                                                                   <  7:296958/149‑1 (MQ=255)   aggagCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTggcgg                                                                                                                                                 <  1:132466/149‑1 (MQ=255)                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                                                                                        >  8:220051/1‑100 (MQ=255)                                                                                                            cTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGTGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTcg                                                                                        <  7:220051/100‑1 (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                                >  8:64395/1‑148 (MQ=255)                                                                                                                    ataCAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCaa                                <  7:64395/148‑1 (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  <  3:102952/149‑1 (MQ=255)                                                                                                                                                 ttGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGtt  >  1:140026/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGAGGAGCGTTTTTGCAGGGTTTGTAGATGACATGGTGCAAGCCTGCGCAATACGCATCGGCGATTACCGATGTCATAGGCGCTCGGTATTGTAAGCAAGCTTTTTACTGATTGTATACAGCTCGAAGTGCCCGTGCCGGCCTCTTGGCGGCTAGAGCCGCGAAGAGGCCGGCACGGGTTGCAGCCATGA‑AAGATCACACCGCCTCGCGGCTGCTCACACCCTCCAGCAACCGGGCAATCCCCAGCGGGTTAGCGTTTTTCAACGCATCCGGCAGCAGTGCGTCAGGGTAGTT  >  NC_002947/2951188‑2951480

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