Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,949,794	A→T	7.4%	intergenic (+59/+54)	PP_2583 → / ← oguA	acetyltransferase/8‑oxoguanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,949,794	0	A	T	7.4%	157.2 / 2.6	54	intergenic (+59/+54)	PP_2583/oguA	acetyltransferase/8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base T (2/2);  total (26/28)
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.13e-01

CAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCACGA  >  NC_002947/2949652‑2949942                                                                                                                                                |                                                                                                                                                     cAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGc                                                                                                                                                >  1:239479/1‑149 (MQ=255)   gAAACTGGGCTACGGATACGTGAAACCGTCGCCGTCGGACCATGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAg                                                                                                                                                   <  8:129026/144‑1 (MQ=255)   gAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACAACGGCGGCATGGCCTGGCATGAGATGGTAATGGACATTTGAGTGGACATGACATCTTAGAGTGCGCTACACGCAGCTTGAC‑TAGGTGGCAGCGAAGCCCAAAgtaga                                                                                                                                               >  7:129026/1‑144 (MQ=255)       cTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCa                                                                                                                                          >  5:117581/1‑149 (MQ=255)        tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGGCATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCAt                                                                                                                                         >  8:320700/1‑149 (MQ=255)        tGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCAt                                                                                                                                         >  7:208561/1‑149 (MQ=255)          ggCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATg                                                                                                                                        <  5:192694/148‑1 (MQ=255)               ccGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGTGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCt                                                                                                                                    <  8:122709/147‑1 (MQ=255)               ccGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCtt                                                                                                                                   <  8:208561/148‑1 (MQ=255)                cGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCttca                                                                                                                                 <  3:35299/149‑1 (MQ=255)                        gAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGc                                                                                                                                       >  8:97116/1‑135 (MQ=255)                        gAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGc                                                                                                                                       <  7:97116/135‑1 (MQ=255)                        gAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCACCCCAAAGGGCTGCGCATGCCCTTCATTCACCCg                                                                                                                         >  7:83343/1‑149 (MQ=255)                         aaaCCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCgg                                                                                                                        >  4:280533/1‑149 (MQ=255)                             cgtcgccgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGGCATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGAt                                                                                                                    <  7:320700/149‑1 (MQ=255)                             cgtcgccgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGAt                                                                                                                    <  6:262337/149‑1 (MQ=255)                                   cgtcgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGGGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAg                                                                                                              <  3:246789/149‑1 (MQ=255)                                      cgGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCt                                                                                                            >  4:276779/1‑148 (MQ=255)                                        gACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGAt                                                                                                         <  2:239479/149‑1 (MQ=255)                                                        ccTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGt                                                                                           <  1:328696/147‑1 (MQ=255)                                                        ccTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTAATTAAACCGGCGAAGCCGAGCAGATGACAGGCCTCAAGTg                                                                                          >  2:328696/1‑148 (MQ=255)                                                                      gTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTg                                                                            <  3:276779/148‑1 (MQ=255)                                                                       tCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGAt                                                                          >  4:134639/1‑149 (MQ=255)                                                                         aTGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                        >  5:29204/1‑149 (MQ=255)                                                                         aTGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATgg                                                                        >  3:132925/1‑149 (MQ=255)                                                                          tGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCGTCATTCACCCGGCGATGGCGAGCTGATGCCAGGCCTCAAGTGCCAATGACCAGCTGATggg                                                                       >  5:53579/1‑149 (MQ=255)                                                                            gACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGc                                                                      <  3:280533/148‑1 (MQ=255)                                                                             aCCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCCTTAGTTGGAACCGAACCCAAAAGGGATGCGCATGCCCTTAATAAACCCGCCGATGCCGACATGATGCCAGGCAACAATTGCCAATGACAACATCATGGGcc                                                                     >  1:81803/1‑149 (MQ=255)                                                                                  tGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCt                                                                 <  8:231308/147‑1 (MQ=255)                                                                                    aGTGGACATGAGATCGTAGTGTGCGCTACCGGCTGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGAATGGCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGtt                                                             <  8:83343/149‑1 (MQ=255)                                                                                      tGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGc                                                                  <  4:51900/142‑1 (MQ=255)                                                                                      tGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGc                                                                  >  3:51900/1‑142 (MQ=255)                                                                                                atCGTAGTGTGCGCTACCCGCGGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAg                                                 >  7:226991/1‑149 (MQ=255)                                                                                                atCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgt                                                           <  4:86094/139‑1 (MQ=255)                                                                                                atCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtgt                                                           >  3:86094/1‑139 (MQ=255)                                                                                                      gtgtGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtg                                                            <  2:318800/132‑1 (MQ=255)                                                                                                      gtgtGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTtg                                                            >  1:318800/1‑132 (MQ=255)                                                                                                             cTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAgg                                    <  3:134639/149‑1 (MQ=255)                                                                                                                 ccGCAGCTTGCC‑TAGGTGGCAGCGCCGCCCAAAGGGCTGCGCAAGCCCATCATTAACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGa                                >  2:302771/1‑149 (MQ=255)                                                                                                                    cAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACgg                             >  7:211380/1‑149 (MQ=255)                                                                                                                     aGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGt                            <  4:132925/149‑1 (MQ=255)                                                                                                                      ggttGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTg                           <  8:226991/147‑1 (MQ=255)                                                                                                                        ttGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGc                         >  5:6884/1‑149 (MQ=255)                                                                                                                                agggcatgcgcagcCCTTTGGGCTGCGCTTGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGc                         <  2:281760/135‑1 (MQ=255)                                                                                                                                agggcatgcgcagcCCTTTGGGCTGCGCTTGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGCCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGc                         >  1:281760/8‑142 (MQ=255)                                                                                                                                 gTGGCAGCGCAGCCCAAAGGGCTGGGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGACCTCAAGTTCCACTGACCAGCTGTAGGGCCAGCTGGGTGTGGCGTTGCAGTACGCGCGGGCGGTCGACGGTGGGCAGCTggg                 >  8:160166/1‑148 (MQ=255)                                                                                                                                   gcatgcgcagcCCTTTGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTAcgcgc                                          <  1:217915/118‑1 (MQ=255)                                                                                                                                   gcatgcgcagcCCTTTGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTAcgcgc                                          >  2:217915/5‑122 (MQ=255)                                                                                                                                      agcgcagcCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGACCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCt            <  7:160166/149‑1 (MQ=255)                                                                                                                                        atccgcCCAAAGGGCTGCGCATGTCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCt            <  3:105982/143‑1 (MQ=255)                                                                                                                                         gcagcCCAAAGGGCTGCGCATGCCCTGCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTc         <  6:117581/149‑1 (MQ=255)                                                                                                                                          cagcCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTca        <  8:211380/149‑1 (MQ=255)                                                                                                                                              ccAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTAc                                              >  6:3154/1‑107 (MQ=255)                                                                                                                                              ccAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTAc                                              <  5:3154/107‑1 (MQ=255)                                                                                                                                                aaaGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCacga  <  8:134466/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CAGAAACTGGGCTACCGATACGTGAAACCGTCGCCGTCGGACCACGGCGGCATGGCCTGGCATGAGATGGTCATGGACCTTTGAGTGGACATGACATCGTAGTGTGCGCTACCCGCAGCTTGCC‑TAGGTGGCAGCGCAGCCCAAAGGGCTGCGCATGCCCTTCATTCACCCGGCGATGCCGAGCTGATGCCAGGCCTCAAGTGCCACTGACCAGCTGATGGGCCAGCTGGTTGTGGCGTTCCAGTACGCGCGGCAGGTCGACGGTGGCCAGCTGGCCGTCTTTCACCACGA  >  NC_002947/2949652‑2949942

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