Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,921,460	G→A	100%	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	100.0%	147.0 / NA	44	G3494R (GGG→AGG)	PP_2561	hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/20);  total (24/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATTT  >  NC_002947/2921333‑2921605                                                                                                                                |                                                                                                                                                  cGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTcc                                                                                                                              >  2:341984/1‑149 (MQ=255) cGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTc                                                                                                                               >  5:140546/1‑148 (MQ=255) cGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTc                                                                                                                               <  6:140546/148‑1 (MQ=255)  ggTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCt                                                                                                                             <  7:351825/149‑1 (MQ=255)             ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                                 <  8:251581/118‑1 (MQ=255)              cTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCgg                                                                                                                  >  3:104559/1‑148 (MQ=255)              cTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGc                                                                                                                 <  8:81245/149‑1 (MQ=255)                  aCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACTTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGt                                                                                                             <  3:285468/149‑1 (MQ=255)                      gggCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACg                                                                                                         <  1:341984/149‑1 (MQ=255)                            tGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACgg                                                                                                   >  7:105096/1‑149 (MQ=255)                                gATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggtggt                                                                                               >  2:362486/1‑149 (MQ=255)                                     cAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtgggtg                                                                                          <  4:17584/149‑1 (MQ=255)                                                gcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAAc                                                                               <  1:362486/149‑1 (MQ=255)                                                  ggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGa                                                                             >  5:182698/1‑149 (MQ=255)                                                             tGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCg                                                                  >  7:31282/1‑149 (MQ=255)                                                                  cccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAcc                                                             >  8:222444/1‑149 (MQ=255)                                                                  cccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGGGGGGCCGGTACCGACACCATGATCGCCAcc                                                             >  6:301371/1‑149 (MQ=255)                                                                   ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCTGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAcc                                                             >  8:833/1‑148 (MQ=255)                                                                           gcggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAAc                                                    >  1:12081/1‑149 (MQ=255)                                                                                gATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAtgt                                               >  4:281089/1‑149 (MQ=255)                                                                                      aaCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCTGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                           >  2:254605/1‑147 (MQ=255)                                                                                       aCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                         <  8:105096/148‑1 (MQ=255)                                                                                           cacaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCg                                    <  8:31282/149‑1 (MQ=255)                                                                                              aTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg                                                                                                               >  2:14180/1‑71 (MQ=255)                                                                                              aTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg                                                                                                               <  1:14180/71‑1 (MQ=255)                                                                                              aTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTg                                 >  3:125817/1‑149 (MQ=255)                                                                                                  ttCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCa                                                                        >  8:76481/1‑106 (MQ=255)                                                                                                  ttCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCa                                                                        <  7:76481/106‑1 (MQ=255)                                                                                                      gCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCttt                         >  2:161409/1‑149 (MQ=255)                                                                                                        gACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGc                                                                                      <  2:35016/86‑1 (MQ=255)                                                                                                        gACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGc                                                                                      >  1:35016/1‑86 (MQ=255)                                                                                                            tGGGTACCTACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGGGGGGCCGGTAACGACACCATGATCGCCACCTATGGAAACGATGTGTTACAGTTCTCTGCAGGCTTTGGCAAc                   <  1:254605/149‑1 (MQ=255)                                                                                                                tACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAcc               <  4:125817/149‑1 (MQ=255)                                                                                                                    gACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGaaa                                                     <  2:215090/107‑1 (MQ=255)                                                                                                                    gACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGaaa                                                     >  1:215090/1‑107 (MQ=255)                                                                                                                         ccTGATCAGGGGGGCTGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTAGCAACGACCTGAtcatca      <  7:833/149‑1 (MQ=255)                                                                                                                         ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGCTGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAAc                   >  8:283544/1‑136 (MQ=255)                                                                                                                         ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAAc                                                    >  7:52034/1‑103 (MQ=255)                                                                                                                         ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAAc                                                    <  8:52034/103‑1 (MQ=255)                                                                                                                         ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCATTTCGCTGCATGCTTTGGCAAc                   <  7:283544/136‑1 (MQ=255)                                                                                                                         ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc       >  3:172085/1‑148 (MQ=255)                                                                                                                          cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATAc                                                                                                     <  5:180200/53‑1 (MQ=255)                                                                                                                          cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATAc                                                                                                     >  6:180200/1‑53 (MQ=255)                                                                                                                             aTCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAttt  >  2:267097/1‑149 (MQ=255)                                                                                                                             aTCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAttt  <  7:222444/149‑1 (MQ=255)                                                                                                                                |                                                                                                                                                  CGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATTT  >  NC_002947/2921333‑2921605

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