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_0029472,921,4600GA97.6% 131.9 / NA 42G3494R (GGG→AGG) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/16);  minor base T (0/1);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATT  >  NC_002947/2921322‑2921604
                                                                                                                                          |                                                                                                                                                 
cTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCgg                                                                                                                                         <  5:358895/149‑1 (MQ=255)
cTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGCGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTTATCAGGGGGGCAgg                                                                                                                                         >  5:223242/1‑149 (MQ=255)
 tGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGCGTACCGACACCCTGATCAGGGGGGCCgg                                                                                                                                         >  5:95265/1‑148 (MQ=255)
      catcccGCTAACCACAACCTCAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGAGGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGa                                                                                                                                   <  6:184330/145‑1 (MQ=255)
       aGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATCACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACAATGTTCTGCGACGAGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGAc                                                                                                                                  >  8:260281/1‑149 (MQ=255)
        ggCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACt                                                                                                                                 >  1:192940/1‑149 (MQ=255)
              tAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                                 >  4:339353/1‑127 (MQ=255)
              tAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                                 <  3:339353/127‑1 (MQ=255)
              tAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCa                                                                                                                          <  4:102731/149‑1 (MQ=255)
                       ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCTATGACACCGTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTTGGTACCGACACCGTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCgg                                                                                                                 <  7:260281/149‑1 (MQ=255)
                              gggggCGCTGGCGTTGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTaa                                                                                                          >  5:115249/1‑149 (MQ=255)
                                  gcgcTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGa                                                                                                       >  5:91071/1‑148 (MQ=255)
                                     cTGGCGATGACAGCCTGGATTGTGTCAACTGAAATGACACCCTTCAGGGCGTCTATGGTAACGACACATTGTTCGGCTACGTGGGTACAGACACCCTGATGTGGGTGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACg                                                                                                   <  6:95265/149‑1 (MQ=255)
                                        gCGATGACAGCCTGGATGGCGGCAACGGCAATGACACACTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACggtg                                                                                                <  7:303701/149‑1 (MQ=255)
                                              aCAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTggtggg                                                                                           <  6:115249/148‑1 (MQ=255)
                                                aGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGgtgggtgg                                                                                        >  1:326765/1‑149 (MQ=255)
                                                aGCCTGGATGGCGGCAACGGCAATGACACACTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTggtgggtg                                                                                         <  4:156972/148‑1 (MQ=255)
                                                aGCATGGCTGCCTGCAACGGCAATTCCACCCTGCAGCGCGCGGCTGGTAAGGACACATTGTTCGTCGACGTTGGTACCGACACCATGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGgtgggtgg                                                                                        <  6:73995/149‑1 (MQ=255)
                                                           cggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGt                                                                                 <  3:189494/145‑1 (MQ=255)
                                                           cggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGt                                                                                 >  4:189494/1‑145 (MQ=255)
                                                               aacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGCCTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGacac                                                                         >  4:110512/1‑149 (MQ=255)
                                                                acggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGacac                                                                         >  4:279436/1‑148 (MQ=255)
                                                                    caaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCGGGTTACGCCACCATg                                                                     >  4:310363/1‑148 (MQ=255)
                                                                     aaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACAATGAt                                                                   >  7:164368/1‑149 (MQ=255)
                                                                     aaTGACACCCTGCAGGGCGGCGATGCTAACGACACTTTGTTCGGCGACTTTTTTAACCACACCCTGCTCAGGGGGGCCGGCAACTACTTC‑CTAAACTGTGCTTCAGTTAACGATACGGTGGTGTTTTTGGCCGGTAACGACACCATGAt                                                                   <  2:192940/149‑1 (MQ=255)
                                                                              cTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGCCTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGCCACCATGATCGCCGCCGa                                                          >  4:269002/1‑149 (MQ=255)
                                                                                 cAGGGCGTCGATGGTCACGACACATTGTTCGGCGACGTGGGTACCGACACCCTTATCAGGGGGGCCGGCAACGTCTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGAt                                                         <  3:310363/147‑1 (MQ=255)
                                                                                  aGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGa                                                      <  2:326765/149‑1 (MQ=255)
                                                                                      cggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAAc                                                   >  7:108031/1‑148 (MQ=255)
                                                                                           aTGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAtgtg                                             >  4:83406/1‑149 (MQ=255)
                                                                                               tAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                          <  3:83406/148‑1 (MQ=255)
                                                                                                      acaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCAGTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCg                                   >  1:309977/1‑149 (MQ=255)
                                                                                                        aTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTg                                >  8:325026/1‑149 (MQ=255)
                                                                                                                               acacGCTGATCAGGGGGGCCGGCAACGACTTC‑GTCAACGGTGCGGCAGGTAACGATAGGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGTTGTGTTCCAGTTCGCTGCa                              <  4:138182/128‑1 (MQ=255)
                                                                                                                               acacCCTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGCTGCCGCAGGTAACGATACGGTGGAGGGTGGGGCCCGTAACGACACCATGATCGCCACCGATCGAAACGATCTGTTCCAGTACGCTGCa                              >  3:138182/1‑128 (MQ=255)
                                                                                                                               acacCCTCATCAGGGGCGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatc         <  7:241199/149‑1 (MQ=255)
                                                                                                                                   ccTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACCCAATGATCGCAACCGATGGAAACGATGTGTTCCATTTCGCTGCAGGATTTGGCAACGACAGGatcatc      >  5:110849/1‑148 (MQ=255)
                                                                                                                                   ccTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc      >  4:1604/1‑148 (MQ=255)
                                                                                                                                   ccTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc      >  4:124926/1‑148 (MQ=255)
                                                                                                                                   ccTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc      >  2:187819/1‑148 (MQ=255)
                                                                                                                                    cTGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCaa    >  7:73696/1‑149 (MQ=255)
                                                                                                                                     tGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAgg                            >  6:206235/1‑124 (MQ=255)
                                                                                                                                     tGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAgg                            <  5:206235/124‑1 (MQ=255)
                                                                                                                                     tGATCAGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGT‑ACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAtt  >  4:355578/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTC‑CTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATT  >  NC_002947/2921322‑2921604

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