Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,913,051 G→A 9.1% A691T (GCC→ACC)  PP_2561 → hemolysin‑type calcium‑binding bacteriocin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,913,0510GA 9.1% 121.1 / 3.2 44A691T (GCC→ACC) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (2/2);  total (21/23)
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.48e-01

CCGCGGTCGTGATACCGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGCCTATGCCAACGTGA  >  NC_002947/2912903‑2913198
                                                                                                                                                    |                                                                                                                                                   
ccGCGGTCGTGATACCGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTg                                                                                                                                                     <  1:11204/149‑1 (MQ=255)
   cgGTCGTGATACCGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTgccg                                                                                                                                                  <  2:306671/149‑1 (MQ=255)
      tCGTGATACCGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTgccgcc                                                                                                                                                >  2:67598/1‑148 (MQ=255)
         tGATACCGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACg                                                                                                                                            >  3:17121/1‑149 (MQ=255)
              ccGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCAc                                                                                                                                        <  7:349534/148‑1 (MQ=255)
                      cccAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCAcagcag                                                                                                                               <  7:22280/149‑1 (MQ=255)
                      cccAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCAcagca                                                                                                                                <  2:204048/148‑1 (MQ=255)
                       ccAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTa                                                                                                                                              <  7:54561/133‑1 (MQ=255)
                        cAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAgcagca                                                                                                                             >  4:368167/1‑149 (MQ=255)
                                  gAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGc                                                                                                                   <  2:164192/149‑1 (MQ=255)
                                    ggCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGcca                                                                                                                 <  6:294662/149‑1 (MQ=255)
                                      cccGCCGCTAGTTCTATTCCTCCACCGGAGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTTCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGccac                                                                                                                <  8:176877/148‑1 (MQ=37)
                                                  tCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTACCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGgccgc                                                                                                    >  6:169265/1‑148 (MQ=255)
                                                  tCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTACCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGgccgc                                                                                                    >  8:211157/1‑148 (MQ=255)
                                                    tATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCaa                                                                                                 >  8:281319/1‑149 (MQ=255)
                                                     aTGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAg                                                                                                >  8:206522/1‑149 (MQ=255)
                                                          gccACCGGCGACAGCCAGCTCAAGCCGTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAgcgcg                                                                                            >  8:130233/1‑148 (MQ=255)
                                                          gccACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTAc                                                                                                                                             <  4:100463/99‑1 (MQ=255)
                                                          gccACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAgcgcg                                                                                            <  7:130233/148‑1 (MQ=255)
                                                            cACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTCATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGCCCACCACCGAGGCCGCCAACCGcgccg                                                                                         >  8:118525/1‑149 (MQ=255)
                                                                 gCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCgccgccgcc                                                                                    >  4:171274/1‑149 (MQ=255)
                                                                   gACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGt                                                                                  <  3:178737/149‑1 (MQ=255)
                                                                     cagccagcTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCACCATTACCGGGGCCACCACCGAGGCCGCCAACCGCGCCGCCGCCGTgg                                                                                >  5:238011/1‑149 (MQ=255)
                                                                           ggtCAAGCGCTACTTGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCTTTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGGCAAGCGCGCGGCCGCCGTGGCACTgg                                                                          <  7:118525/147‑1 (MQ=255)
                                                                               aaGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTAc                                                                       <  2:207496/148‑1 (MQ=255)
                                                                               aaGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTAc                                                                       >  1:207496/1‑148 (MQ=255)
                                                                               aaGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACt                                                                      >  8:19989/1‑149 (MQ=255)
                                                                                   ccTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGc                                                                  <  7:19989/149‑1 (MQ=255)
                                                                                    cTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCg                                                                 <  1:67503/149‑1 (MQ=255)
                                                                                       cATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTACCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGt                                                               <  7:211157/148‑1 (MQ=255)
                                                                                       cATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTACCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGt                                                               <  5:169265/148‑1 (MQ=255)
                                                                                           aGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGt                                                                                  >  8:21563/1‑125 (MQ=255)
                                                                                           aGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGt                                                                                  <  7:21563/125‑1 (MQ=255)
                                                                                           aGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACgg                                                          >  3:72596/1‑149 (MQ=255)
                                                                                               gggTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCa                                                      >  3:361031/1‑149 (MQ=255)
                                                                                                ggTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCAc                                                     <  4:17121/149‑1 (MQ=255)
                                                                                                     aCTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTc                                                <  3:171274/149‑1 (MQ=255)
                                                                                                      cTTTGCGGACCACCTCTACCCCCCGGGGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTc                                                <  6:238011/148‑1 (MQ=38)
                                                                                                            cGACCACCTCAAGCACCCGGCGTCGCTGATAAACTTCATTGCCGCCTACGGCACCCACAGCAGAATTAACGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCCTGGCAAGGGTACATAGAGGTGAGGGCGCCACCCTCgccagcg                                         >  8:126672/1‑143 (MQ=255)
                                                                                                              accaccTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTc                                       <  1:47158/149‑1 (MQ=255)
                                                                                                                      aaGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTc                                 <  7:206522/147‑1 (MQ=255)
                                                                                                                      aaGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTc                                 <  3:368167/147‑1 (MQ=255)
                                                                                                                                    gatGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTg                 >  5:102276/3‑149 (MQ=255)
                                                                                                                                       gATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGcc               >  2:372586/1‑148 (MQ=255)
                                                                                                                                       gATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGcc               >  4:178754/1‑148 (MQ=255)
                                                                                                                                        aTCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGCCTa             >  2:323824/1‑149 (MQ=255)
                                                                                                                                                   tGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGCCTATGCCAACGTGa  >  7:326263/1‑149 (MQ=255)
                                                                                                                                                    gccgccTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGCCTAt            >  8:199491/1‑138 (MQ=255)
                                                                                                                                                    gccgccTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGCCTAt            <  7:199491/138‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGCGGTCGTGATACCGGCATCCCCAGCCTGAACGAGGCCCGCCGCGAGTTCTATGCCGCCACCGGCGACAGCCAGCTCAAGCCCTACATGAGCTGGGTCGACTTTGCCGACCACCTCAAGCACCCGGCGTCGCTGATCAACTTCATTGCCGCCTACGGCACCCACAGCAGCATTACCGGGGCCACCACCGAGGCCGCCAAGCGCGCCGCCGCCGTGGCACTGGTACTTGGCGGTGACGGCGCACCCCTCGACCGCCTCGACTTCCTCAACGGCACCGGTGCCTATGCCAACGTGA  >  NC_002947/2912903‑2913198

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