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

ACACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATTT  >  NC_002947/2921314‑2921605                                                                                                                                                    |                                                                                                                                                  acacGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCagg                                                                                                                                                 >  6:203502/1‑149 (MQ=255)  cacGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCaggg                                                                                                                                                <  5:203502/149‑1 (MQ=255)               caGGCCCGAAACGACAACCTGAACGGGGGCGCTGGCCAAGACAGCCTGGATGGCGGCAACCGCAATGACACCCTGCAGGGCGGACATGGTAACGACACATTGTTCGGCTACCTGGGTACCGCCACCCAGA‑TCAGCGCGGCCGGCAACGa                                                                                                                                   >  6:333040/2‑149 (MQ=255)                   ccGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTc                                                                                                                               <  2:363965/149‑1 (MQ=255)                    cGGTAACGACAACCTGAACGGGGGCGCAGGCGATGACAGCCTGCATGGAGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACTACTTcc                                                                                                                              >  1:18848/1‑149 (MQ=255)                         ccGACATCCTGATCTGGGGCGCTGGCGATGACCTCCTGGATGTCGGCAACGTCAATGACACCCTGCAGGGCGGCGATGGTTACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCaa                                                                                                                          <  8:171819/147‑1 (MQ=255)                            aCAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGt                                                                                                                      >  7:242859/1‑149 (MQ=255)                              aaCCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGc                                                                                                                    >  1:333109/1‑149 (MQ=255)                                  tGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCgg                                                                                                                  >  7:64512/1‑147 (MQ=255)                                               tGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACgg                                                                                                   <  1:85434/149‑1 (MQ=255)                                                   gATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTcc                                                                                                                              >  2:47425/1‑118 (MQ=255)                                                   gATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTcc                                                                                                                              <  1:47425/118‑1 (MQ=255)                                                    aTGACAGCCTGGATGGCGGCAACGGCAATGACCCCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCCACGGTGCGGCAGGTACCGATCCGggggtg                                                                                              >  2:366888/1‑149 (MQ=255)                                                     tGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGtggtgg                                                                                             <  8:64512/149‑1 (MQ=255)                                                           ccTTGATGGCGGCAACGGGAATGACACCCTTCATGGCGGCGTTTGTAACGACAAATTGTTCGGCGACGTTGGTAACGTAACCCTGATTCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGa                                                                                                        <  1:299282/133‑1 (MQ=255)                                                           ccTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGa                                                                                                        >  2:299282/1‑132 (MQ=255)                                                            cTGGATGGCGGCAAGGGCATTGACACCGTGCTGGGCGGCGATGGTAACGTCACATTGTTCGGCGAGGTGGGTTGCGTGACCGTGA‑TGAGGGGGGCGGGGAACGACTTCGTCTTCGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGc                                                                                      <  5:333040/149‑1 (MQ=255)                                                                   gcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTaa                                                                                <  5:318937/148‑1 (MQ=255)                                                                   gcggcAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAAc                                                                               <  2:333109/149‑1 (MQ=255)                                                                     ggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGa                                                                             >  3:241052/1‑149 (MQ=255)                                                                           ggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCAt                                                                       >  3:152020/1‑149 (MQ=255)                                                                             caaTGAAACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGTGGCAGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGa                                                                     >  3:165733/1‑149 (MQ=255)                                                                                tGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCg                                                                  >  7:221410/1‑149 (MQ=255)                                                                                    acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                              >  2:18353/1‑149 (MQ=255)                                                                                          cAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATgg                                                        <  8:221410/149‑1 (MQ=255)                                                                                                    aTGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAtgtg                                              <  4:152020/149‑1 (MQ=255)                                                                                                     tGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATgtgt                                             <  3:194531/149‑1 (MQ=255)                                                                                                         aaCGACACTATCCTCGGCGACGTAGTTACCTACACCCTCA‑GCCGGTCTGCCGGCAACTACATCCTCAACGGTTCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                         <  2:399703/149‑1 (MQ=255)                                                                                                             acacATTGTTCGGCTACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGTTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACCATGTGTTCCAGtt                                      >  2:154370/1‑148 (MQ=255)                                                                                                                  ttGTTCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTg                                 >  8:247202/1‑148 (MQ=255)                                                                                                                     ttCGGCGACGTGGGTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAgg                             >  1:267312/1‑149 (MQ=255)                                                                                                                                gggTACCGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCCGGCTTTGGCAACg                  >  6:277652/1‑149 (MQ=255)                                                                                                                                      cGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGa            <  8:242859/149‑1 (MQ=255)                                                                                                                                      cGACACCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTCGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGAAGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAcc               >  7:103619/1‑146 (MQ=255)                                                                                                                                      agACACCCTGA‑TCAGGGGTGCCGGCATCGACTTCCTCATCGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACGGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAcc               <  8:103619/145‑1 (MQ=255)                                                                                                                                        acacCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGat           >  3:195678/1‑148 (MQ=255)                                                                                                                                          acCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat        <  2:267312/149‑1 (MQ=255)                                                                                                                                          acCCTGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatca         >  5:141886/1‑148 (MQ=255)                                                                                                                                              tGA‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAt    >  8:57656/1‑149 (MQ=255)                                                                                                                                                a‑TCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAttt  >  7:24540/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ACACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGA‑TCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATTT  >  NC_002947/2921314‑2921605

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