Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	129,229	C→A	18.1%	P569P (CCC→CCA)	polA →	DNA polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	129,229	0	C	A	18.1%	61.7 / 2.4	22	P569P (CCC→CCA)	polA	DNA polymerase I
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base A (2/2);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

GCCAGGGTGCACTGGTCGATGCCGAGCTGTTGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGCCAGAGG  >  NC_002947/129105‑129359                                                                                                                             |                                                                                                                                   gcCAGGGTGCACTGGTCGATGCCGAGCTGTTGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTggg                                                                                                                         >  1:103329/1‑136 (MQ=255)             tGGTCGATGCCGAGCTGTTGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTa                                                                                                             <  1:40939/136‑1 (MQ=255)                   aTGCCGAGCTGTTGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACaa                                                                                                       >  2:444584/1‑136 (MQ=255)                   aTGCCGAGCTGTTGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACaa                                                                                                       >  1:423584/1‑136 (MQ=255)                              ttGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGc                                                                                            <  2:423584/136‑1 (MQ=255)                                        agagCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAg                                                                                  >  1:127851/1‑136 (MQ=255)                                                         gTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGgcc                                                                 >  2:115467/1‑136 (MQ=255)                                                         gTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGgcc                                                                 >  1:471060/1‑136 (MQ=255)                                                                   cTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCAc                                                       <  1:115467/136‑1 (MQ=255)                                                                   cTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCAc                                                       <  1:444584/136‑1 (MQ=255)                                                                          ggagctgCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGTCCAGCCATCCACCGCCGaa                                                <  2:510073/136‑1 (MQ=255)                                                                               tgCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGt                                           >  1:373119/1‑136 (MQ=255)                                                                                gCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTg                                          >  2:1576/1‑136 (MQ=255)                                                                                    gCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTcg                                      >  1:397539/1‑136 (MQ=255)                                                                                                 ccGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTCGACGAACTGGCCtt                         <  2:138748/136‑1 (MQ=255)                                                                                                  cGGTGAGACGTTCAACCTCGGCTCGCCAAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACAGCCGAAGCAGTGCTCGACGAACTGGCCTTg                        >  1:462150/1‑136 (MQ=255)                                                                                                  cGGTGAGACGTTCAACCTCGGCTCGCCAAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACAGCCGAAGCAGTGCTCGACGAACTGGCCTTg                        >  1:421151/1‑136 (MQ=255)                                                                                                                  cTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGc        <  1:1576/136‑1 (MQ=255)                                                                                                                  cTCGGCTCGCCAAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACAGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGc        <  2:421151/136‑1 (MQ=255)                                                                                                                  cTCGGCTCGCCAAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACAGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGc        <  2:462150/136‑1 (MQ=255)                                                                                                                   tCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGcc       <  2:103329/136‑1 (MQ=255)                                                                                                                        tCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGCCAGAgg  <  2:373119/136‑1 (MQ=255)                                                                                                                             |                                                                                                                                   GCCAGGGTGCACTGGTCGATGCCGAGCTGTTGAAAGTACAGAGCGGTGAGCTGGGCGTGAAGATGGCTGAGCTGGAGCTGCGTGCCTACGAGCTCGCCGGTGAGACGTTCAACCTCGGCTCGCCCAAGCAGCTGGGCACGATCCTTTACGACAAGCTGGGCATGCCGGTGTTGAGCAAGACCGCCAAGGGCCAGCCATCCACCGCCGAAGCAGTGCTCGACGAACTGGCCTTGCAAGGCTACCCGCTGCCAGAGG  >  NC_002947/129105‑129359

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