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

ACACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATTT  >  NC_002947/2921314‑2921605
                                                                                                                                                  |                                                                                                                                                 
acacGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACTGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                                 >  7:320493/1‑149 (MQ=255)
    gTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagggggg                                                                                                                                             >  4:130424/1‑149 (MQ=255)
           ggggAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCaa                                                                                                                                      <  3:227671/149‑1 (MQ=255)
                ggCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACt                                                                                                                                 <  3:130424/149‑1 (MQ=255)
                          gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACgg                                                                                                                       >  6:168522/1‑149 (MQ=255)
                          gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATCACACCCTGCAGGTCGGCGATGGTAACGAAACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACgg                                                                                                                       >  6:229803/1‑149 (MQ=255)
                               ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCgg                                                                                                                  >  3:267102/1‑149 (MQ=255)
                                  gAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg                                                                                                               >  3:120070/1‑149 (MQ=255)
                                                      aCAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtgggt                                                                                           <  5:168522/149‑1 (MQ=255)
                                                       cAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGgtggtg                                                                                              >  1:106205/1‑145 (MQ=255)
                                                       cAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGgtggtg                                                                                              <  2:106205/145‑1 (MQ=255)
                                                        aGCCTGGATGGCGGCAACGGCAATGACACACTGCAGGGCGGAGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCAGGCAACGACTTCCTAAACGGAGCGGCAGGTAACGAGACGGAGgagggtgg                                                                                         >  6:287275/1‑149 (MQ=255)
                                                           atgtaTGGCGTCAACGGCAATGACACCATGCATGGAGGCTATGGTATCTAAACATTGTTCGGCGACGTGGGTACCGACACCGTGATCTGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGc                                                                                      <  8:320493/145‑1 (MQ=255)
                                                              gATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTgggg                                                                                       <  5:175917/145‑1 (MQ=255)
                                                              gATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTgggg                                                                                       >  6:175917/1‑145 (MQ=255)
                                                                        acggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACAcc                                                                         <  4:120070/149‑1 (MQ=255)
                                                                          ggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCAt                                                                       >  8:329238/1‑149 (MQ=255)
                                                                           gcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCAt                                                                       >  3:337532/1‑148 (MQ=255)
                                                                              aTGACAAACTTCAGGGCGGCGATGGTAACGACACATTTTTCTTCGACGTGGCTACCGACACCCTTATCAGGGGGGCCGTCAACGACTTCCTCAACGGTGCTGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATc                                                                   <  2:19237/149‑1 (MQ=255)
                                                                                  cacCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCa                                                               >  7:291082/1‑149 (MQ=255)
                                                                                   acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                              >  6:179339/1‑149 (MQ=255)
                                                                                   acCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                              >  8:162800/1‑149 (MQ=255)
                                                                                         cAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGAt                                                          <  4:267102/147‑1 (MQ=255)
                                                                                          aGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGAAGGTAACGATACGGTGGTGGGTGGGGCAGGTAACGACCACATGATCGACACCGATgg                                                        >  3:73789/1‑148 (MQ=255)
                                                                                                gcgATGGTAACGACACATTGTTCGGCGACGTGGGGACCGACACCGTTAGCAGGGGGGCCGGCAACTACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGACCGGTAACGACACCATGATCGCCACCGATGGAAACGAt                                                 <  5:229803/149‑1 (MQ=255)
                                                                                                         aCGACACATTGTTCGGCGACGTGGGTACAGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                           >  5:144792/1‑146 (MQ=255)
                                                                                                         aCGACACATTGTTCGGCGACGTGGGTACAGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                         <  6:144792/148‑1 (MQ=255)
                                                                                                              acaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCATCAACGGTGCGGAAGGTAACGATACGTTGGTGGGTGGGGACGGTAACGACACCATGATCGCCAACGATGGAAACTATGTGTTCCAGTTCGc                                   >  4:73465/1‑149 (MQ=255)
                                                                                                                    ttCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAgg                             >  4:44370/1‑149 (MQ=255)
                                                                                                                        gCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCttt                         <  7:162800/149‑1 (MQ=255)
                                                                                                                               gggTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAAc                   >  3:351737/1‑148 (MQ=255)
                                                                                                                                 gTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAc                <  4:351737/149‑1 (MQ=255)
                                                                                                                                  tACCGTCACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAcc               >  2:139288/1‑149 (MQ=255)
                                                                                                                                     cGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGc                                   <  4:12842/126‑1 (MQ=255)
                                                                                                                                     cGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGc                                   >  3:12842/1‑126 (MQ=255)
                                                                                                                                       acacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGAAGGCTTTGGCAACGACCtcag           >  8:282042/1‑145 (MQ=255)
                                                                                                                                         cccatGATCAGGGTGGCCTTCAACGACTTCCTAAACTGTTCGGCAGGTAACGTTACGGTGGTGGGTGGGGCCGGTAACTAAACCATGATCGCCACCGATGGAAACTATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat        <  3:73465/145‑1 (MQ=255)
                                                                                                                                         acCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGaaa                                                     >  2:213956/1‑104 (MQ=255)
                                                                                                                                         acCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGaaa                                                     <  1:213956/104‑1 (MQ=255)
                                                                                                                                            cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCaa     >  5:208611/1‑149 (MQ=255)
                                                                                                                                               aTCAGGGTGTCCGTCAACGACTTCCTCAACGGTTCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAttt  <  5:287275/149‑1 (MQ=255)
                                                                                                                                               aTCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAttt  >  8:206171/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACACGTTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATTT  >  NC_002947/2921314‑2921605

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