Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	2,921,460	G→A	G3494R (GGG→AGG)	PP_2561 →	hemolysin‑type calcium‑binding bacteriocin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,921,460	0	G	A	93.3%	41.4 / NA	16	G3494R (GGG→AGG)	PP_2561	hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (10/4);  minor base C (0/1);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAA  >  NC_002947/2921316‑2921602                                                                                                                                                 |                                                                                                                                               acGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCaggg                                                                                                                                             >  3:13282/1‑148 (MQ=255)          ggggAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCa                                                                                                                                    >  6:164271/1‑148 (MQ=255)                  cGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTc                                                                                                                            <  4:13282/148‑1 (MQ=255)                         gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACg                                                                                                                     >  5:51372/1‑148 (MQ=255)                                             tGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACg                                                                                                 >  4:49567/1‑148 (MQ=255)                                                                            aaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGCCACCCTGATCAGGGGGGCCGGCAACGACTTCCTCACCGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGa                                                                  >  3:159276/1‑148 (MQ=255)                                                                                    ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                           >  1:28432/1‑147 (MQ=255)                                                                                              ggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACg                                                >  5:164931/1‑148 (MQ=255)                                                                                                    ggTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATgtgt                                          <  5:164271/148‑1 (MQ=255)                                                                                                         cGACACATTGTTCGGCGCCGTGGGCACCGACACCCTGAGCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                      <  8:31558/147‑1 (MQ=255)                                                                                                               aTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCt                               >  1:177825/1‑148 (MQ=255)                                                                                                               aTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCt                               >  4:128763/1‑148 (MQ=255)                                                                                                                             tGGGCACCGACACCCTGATCAGCGGTGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCaa                 <  5:77217/148‑1 (MQ=255)                                                                                                                                 tACCGACACCCTTCTCTGGGGGGCCGGCAACGACTTCCTCAACGGGGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAc             <  8:64535/148‑1 (MQ=255)                                                                                                                                          ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat     >  8:8832/1‑147 (MQ=255)                                                                                                                                            tGATCCGGGGGGCCGGCAAGGTCGTCCTCCCCGGTGCGGCAGGTAACGAGACGGTGGTGGGGGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCaa  <  7:99397/148‑1 (MQ=255)                                                                                                                                              aTCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTg                     >  4:70982/1‑127 (MQ=255)                                                                                                                                                 |                                                                                                                                               ACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAA  >  NC_002947/2921316‑2921602

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