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	97.7%	139.4 / NA	44	G3494R (GGG→AGG)	PP_2561	hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/16);  minor base T (0/1);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921322‑2921603                                                                                                                                            |                                                                                                                                                ccgggggAGCCCGGTCACGACAACCTGAACGCCGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCgg                                                                                                                                       <  2:399426/147‑1 (MQ=255)  tGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGTCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCATGGGGGCCgg                                                                                                                                       >  7:27356/1‑148 (MQ=255)                 aCGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCaa                                                                                                                        <  3:134007/148‑1 (MQ=255)                  cgtcAACCTTTACGGCGGCGCTGGCTATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGTCGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTTAACg                                                                                                                      <  6:174057/146‑1 (MQ=255)                   gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACgg                                                                                                                     >  4:323703/1‑149 (MQ=255)                     cAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGTCAACGTCAATGACACCCTGCAGGGCTGCGATGGTAACGACACATTGTTCGGCTACGTGTGTACCG‑ACACCCTGATCTGGTGGGCCGGCAACGACTTCCTCAACGGTg                                                                                                                   <  8:27356/149‑1 (MQ=255)                        ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCgg                                                                                                                >  7:6335/1‑149 (MQ=255)                            aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGc                                                                                                                                      <  3:377285/123‑1 (MQ=255)                            aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGc                                                                                                                                      >  4:377285/1‑123 (MQ=255)                            aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg                                                                                                             <  5:206515/148‑1 (MQ=255)                            aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGCGGCCGGCAACGACTTCCTCAACGGTGCGGCAgg                                                                                                            >  8:1003/1‑149 (MQ=255)                               gggggCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑AGACCCTCATCAGGGGGGCCGGCAACGACTAGCGCCACGGTGCGGAAGGTaa                                                                                                         >  8:341221/1‑149 (MQ=255)                                     gcTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATAc                                                                                                   >  2:265702/1‑149 (MQ=255)                                       tGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACgg                                                                                                 >  1:217482/1‑149 (MQ=255)                                           gATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggtggt                                                                                             >  5:236339/1‑149 (MQ=255)                                               aCAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtgggt                                                                                         <  8:6335/149‑1 (MQ=255)                                               aCAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtggg                                                                                          <  6:127762/148‑1 (MQ=255)                                                   ccTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGCTTACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtgggt                                                                                         <  1:347292/146‑1 (MQ=255)                                                   ccTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTgggg                                                                                     >  4:213039/1‑149 (MQ=255)                                                           gcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAAc                                                                             >  4:400633/1‑149 (MQ=255)                                                           gcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAAc                                                                             <  6:236339/149‑1 (MQ=255)                                                             ggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGa                                                                           <  7:1003/149‑1 (MQ=255)                                                                   ggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCAt                                                                     >  4:123628/1‑149 (MQ=255)                                                                          acacCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGGGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGc                                                               <  7:341221/148‑1 (MQ=255)                                                                            acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                            >  4:411125/1‑149 (MQ=255)                                                                            acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                            <  3:123628/149‑1 (MQ=255)                                                                            acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                            >  5:217847/1‑149 (MQ=255)                                                                                               gTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATg                                                                    >  7:329860/1‑122 (MQ=255)                                                                                               gTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                         <  3:323703/149‑1 (MQ=255)                                                                                                tAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                         >  7:371039/1‑148 (MQ=255)                                                                                                             ttCGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTACTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCCGTAAAGACACCATGATCGCCACAGATGAAAACGATCTGTACCAGTTCGCTCCaag                           >  7:380907/1‑147 (MQ=255)                                                                                                               cGGCGACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCt                         >  8:49905/1‑149 (MQ=255)                                                                                                                   gACGTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTgg                     >  4:367822/1‑149 (MQ=255)                                                                                                                      gTGGGTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCa                   >  4:423982/1‑148 (MQ=255)                                                                                                                         ggTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTGCAGTTCGCTGCAGGCTTTGGCAACGa               >  6:216514/1‑149 (MQ=255)                                                                                                                          gTACCG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAc              <  3:411125/149‑1 (MQ=255)                                                                                                                              cG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGa          >  7:370282/1‑149 (MQ=255)                                                                                                                              cG‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGa          >  8:316098/1‑149 (MQ=255)                                                                                                                               g‑ACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGat         <  2:217482/149‑1 (MQ=255)                                                                                                                                 acacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGat         >  4:183687/1‑148 (MQ=255)                                                                                                                                  cacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatca       <  3:213039/149‑1 (MQ=255)                                                                                                                                   acCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatca       >  5:127087/1‑148 (MQ=255)                                                                                                                                      cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCaa   >  1:93226/1‑149 (MQ=255)                                                                                                                                       tGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAt  >  7:257964/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                CTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCG‑ACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921322‑2921603

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