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	95.2%	133.4 / 0.6	43	coding (10480/10827 nt)	PP_2561	hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/16);  minor base . (1/1);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

ACACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921314‑2921603                                                                                                                                                    |                                                                                                                                                acacGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCagg                                                                                                                                               <  6:258508/149‑1 (MQ=255)  cccggtGCTGGGCGCCGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCaggg                                                                                                                                              <  3:14342/144‑1 (MQ=255)     gTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCagggggg                                                                                                                                           >  2:53682/1‑149 (MQ=255)          tGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCgg                                                                                                                                       >  2:198644/1‑148 (MQ=255)          tGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCgg                                                                                                                                       >  2:10947/1‑148 (MQ=255)          cgggggACGCCGGCAACCCCAACCTGCACGGGGGCCCTGGCGATGCCAGCCTGGATGGCGGCAACCGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGCGGTACCGACACCCTGATCAGGGGGGCCgg                                                                                                                                       <  1:198644/148‑1 (MQ=255)                         aCGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCaa                                                                                                                        <  1:53682/148‑1 (MQ=255)                           gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGCTGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCATCAACgg                                                                                                                     >  1:252703/1‑149 (MQ=255)                              aaCCTGAACGGGGGAGCTGCCGATGACAGCCTGGATGGCGGCAACCGCAATGACACCCTGCAGGCCCGCGATGGAAACGACACATTGTACGCCCACGTC‑GGTACCGACACCCTCATCAGGGGGGCCGGCAACCACTTCCTCAACGCTGc                                                                                                                  >  2:304729/1‑149 (MQ=255)                                ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCgg                                                                                                                >  5:294803/1‑149 (MQ=255)                                                 gCGAGGACCGCCTGGATGGCGGCACCGGCAATGACACCCTGCAGGGCGGCGATGGTACCGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCCGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggtg                                                                                               <  1:170111/149‑1 (MQ=255)                                                       aCAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtgggt                                                                                         >  5:92191/1‑149 (MQ=255)                                                         agtctcgATTGCGGCAAAGGCAATGACATCGGGCTGTGGGGCCATGCTAAGGACACATTGTGCGGCGACGAG‑GGTATCTTCTGCCTGTTCAGGGGGGCCGGCATCGACTTCCTCAACGGTGCGGCAGGTTACGATACGGTGgtgggtgg                                                                                       <  8:308966/143‑1 (MQ=255)                                                         aGCCTGGATTGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGGACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGgtgggtgg                                                                                       <  2:237648/149‑1 (MQ=255)                                                         aGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGgtgggtgg                                                                                       <  6:92191/149‑1 (MQ=255)                                                         aGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGgtgggtgg                                                                                       >  5:290396/1‑149 (MQ=255)                                                                  ggcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACAGATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGgtggtg                                                                                            >  1:74585/1‑135 (MQ=255)                                                                  ggcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGgtggtg                                                                                            <  2:74585/135‑1 (MQ=255)                                                                    cggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAAc                                                                             >  2:222469/1‑148 (MQ=255)                                                                             caaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATC‑AGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGAt                                                                  >  5:35802/1‑149 (MQ=255)                                                                                      ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATa                                                                                                    >  4:272586/1‑107 (MQ=255)                                                                                      ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATa                                                                                                    <  3:272586/107‑1 (MQ=255)                                                                                            gggCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGa                                                     >  1:333709/1‑148 (MQ=255)                                                                                            gggCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGa                                                     >  1:338862/1‑148 (MQ=255)                                                                                                 gcgATGGTAACGGCACATTGTTCGGCGAGGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGGGCGGCAGGGAACGATGCGGTGGTGGGTGGGGCCGGGAGCGAGACCAGGATCGCCAGCGAGGGGAAGGag                                               >  4:174940/1‑148 (MQ=255)                                                                                                 gcgATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAt                                               <  6:294803/149‑1 (MQ=255)                                                                                                 gcgATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGAACGCCACCGATGGAAACGAt                                               <  1:10947/149‑1 (MQ=255)                                                                                                   gATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAtgt                                             <  1:222469/149‑1 (MQ=255)                                                                                                         aaCGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCAGGCTACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTTACGACACCATGATCGCAACCGATGGAAACGAGGTGtta                                         >  3:249492/1‑146 (MQ=255)                                                                                                                caTTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAt                                               <  2:234678/134‑1 (MQ=255)                                                                                                                caTTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAt                                               >  1:234678/1‑134 (MQ=255)                                                                                                                caTTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCt                                >  4:119496/1‑149 (MQ=255)                                                                                                                 aTTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGAt                                                        <  5:284876/124‑1 (MQ=255)                                                                                                                 aTTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGAt                                                        >  6:284876/1‑124 (MQ=255)                                                                                                                 aTTGTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCt                                >  7:145473/1‑148 (MQ=255)                                                                                                                    gTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCa                             <  2:333709/148‑1 (MQ=255)                                                                                                                    gTTCGGCGACGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCa                             <  2:338862/148‑1 (MQ=255)                                                                                                                        ggCGACGTG‑GGTACCGACACCCTGATC‑AGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCttt                       <  6:35802/149‑1 (MQ=255)                                                                                                                         gCGTCGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACg                                                                                                  >  4:123638/1‑74 (MQ=255)                                                                                                                         gCGTCGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACg                                                                                                  <  3:123638/74‑1 (MQ=255)                                                                                                                             cGTG‑GGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGc                    >  2:6131/1‑148 (MQ=255)                                                                                                                                g‑ggTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGTCCGGTAACGACACCATGATCGCAACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCACCg                >  4:187569/1‑149 (MQ=255)                                                                                                                                             ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCTGTTCGCTGCAGGCGTTGGCAACGACCTGatcatc     >  4:320255/1‑148 (MQ=255)                                                                                                                                             ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc     >  6:326468/1‑148 (MQ=255)                                                                                                                                             ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc     >  2:24221/1‑148 (MQ=255)                                                                                                                                                 aTCAGGGGGGCCGGCAACGACTTCCTCATCGGTGCGGCAGGTAACGATACGGTTGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACTATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAt  <  3:187569/147‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                ACACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTG‑GGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921314‑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: ‑