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	96.0%	76.5 / ‑6.4	25	G3494R (GGG→AGG)	PP_2561	hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/11);  minor base C (0/1);  total (13/12)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATC  >  NC_002947/2921316‑2921600                                                                                                                                                 |                                                                                                                                              acGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCTATGACAGCCTGGATGGCGGCAACGGAAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGAGACGTGGGTACCGACACCCTGATCagggg                                                                                                                                           >  1:230980/1‑149 (MQ=255)  cGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCaggggg                                                                                                                                          <  1:104993/149‑1 (MQ=255)       cTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCg                                                                                                                                      <  4:280370/148‑1 (MQ=255)           gggAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCaggggg                                                                                                                                          >  7:37840/1‑140 (MQ=255)           gggAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCaggggg                                                                                                                                          <  8:37840/140‑1 (MQ=255)           gccccgCCGTTATCGACAACCGCAACGGGGCCGCTGGCGATGACAGCCTGGATGGCGTCAACTGCAATGACACCCTGCAGTGCGTCGTTGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑Ac                                                                                                                                <  2:41694/144‑1 (MQ=255)                         gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGa                                                                                                                              >  4:47981/1‑137 (MQ=255)                         gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGa                                                                                                                              <  3:47981/137‑1 (MQ=255)                         gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACgg                                                                                                                  >  4:320327/1‑149 (MQ=255)                         gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACgg                                                                                                                  >  4:181456/1‑149 (MQ=255)                           cAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTg                                                                                                                <  3:320327/149‑1 (MQ=255)                                  aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAgg                                                                                                         >  2:21529/1‑149 (MQ=255)                                  aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAg                                                                                                          <  7:29633/148‑1 (MQ=255)                                   aCGGCGGCGCTGGCGATGACAGCCTTGATTTCTGCAACGGCAATGACACCCTGCATGGCTGCGATTGTAACGACACATTTTTCGGCGACGTGTGTAACGACACCCTGATCAGGGTGGCCGGAA‑ACGACTTCCTCAACGGTGCGGCAGGt                                                                                                        <  2:230980/149‑1 (MQ=255)                                       gggCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                             <  6:174704/109‑1 (MQ=255)                                       gggCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                             >  5:174704/1‑109 (MQ=255)                                         gcgcTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCATACGACTTCCTCAACGGTGCGGCAGGTAACGa                                                                                                   <  1:226757/149‑1 (MQ=255)                                                             gagggcgtcAACGCCAACTAAAGCCGGCAGCGCGGCTATGGTAACGACACATTGGTTGGCGACGTGTGCACCGACCCCCTGATCCGGCTGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCgg                                                                              <  1:323369/146‑1 (MQ=255)                                                                                  acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                         >  1:180878/1‑149 (MQ=255)                                                                                    ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAcc                                                        >  3:162227/1‑148 (MQ=255)                                                                                                       aaCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAtgt                                          >  3:350124/1‑143 (MQ=255)                                                                                                       aaCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAtgt                                          <  4:350124/143‑1 (MQ=255)                                                                                                             acaTTGTTCGGCGACGTGGGTACCGACAACCTGATCAGGGGGGCCGGCA‑ACGACTTc                                                                                                                          >  8:77335/1‑57 (MQ=255)                                                                                                                   ttCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAg                         <  3:181456/148‑1 (MQ=255)                                                                                                                        cGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTg                   <  2:180878/149‑1 (MQ=255)                                                                                                                                    cGACACCCTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGa       >  5:166799/1‑149 (MQ=255)                                                                                                                                         cccTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat   >  2:232005/1‑148 (MQ=255)                                                                                                                                          ccTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc  >  3:136127/1‑148 (MQ=255)                                                                                                                                          ccTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTGTGGCTACGCCGTGatcatc  >  4:248407/1‑148 (MQ=255)                                                                                                                                           cTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGtt                                       <  4:3737/110‑1 (MQ=255)                                                                                                                                           cTGATCAGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGtt                                       >  3:3737/1‑110 (MQ=255)                                                                                                                                                 |                                                                                                                                              ACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCA‑ACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATC  >  NC_002947/2921316‑2921600

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