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

TGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921323‑2921603
                                                                                                                                         |                                                                                                                                                  
tGGGGCCGGCCGGTTCCGACAACATGTTCTGGGGCGCTGGCGATGACAGCCTGGATTGCGGCAACGGCAATTACACCCTGCAGGGCTGCGATTGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGcc                                                                                                                                            <  2:381181/146‑1 (MQ=255)
            gTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCCACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCATCGACTTcc                                                                                                                               >  8:272145/1‑147 (MQ=255)
                    aaCCTGAACGGTGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGt                                                                                                                       <  8:294133/147‑1 (MQ=255)
                     aCCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTg                                                                                                                      <  4:363955/147‑1 (MQ=255)
                       cTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCg                                                                                                                    >  6:172155/1‑147 (MQ=255)
                          aaCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCa                                                                                                                 >  3:157137/1‑147 (MQ=255)
                                    cTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGCCACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑ACCGATa                                                                                                    >  6:231891/1‑147 (MQ=255)
                                        cGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCCGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACggt                                                                                                <  1:170583/147‑1 (MQ=255)
                                          aTGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACggtgg                                                                                              >  2:349237/1‑147 (MQ=255)
                                                 ccTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGgtgggtgg                                                                                       >  7:75390/1‑147 (MQ=255)
                                                  cTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGtgggtggg                                                                                      >  4:356839/1‑147 (MQ=255)
                                                  cTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGtgggtggg                                                                                      >  5:255272/1‑147 (MQ=255)
                                                                 ggcaaTGACACCCTGCAGGGCGGCGATGGTACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCa                                                                      >  6:328185/1‑147 (MQ=255)
                                                                    aaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGGGGCAGGGTGAAAATATCCGGTGGGGGGTGGGGCCGGTAACGACAcc                                                                       >  1:323344/1‑147 (MQ=255)
                                                                    aaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATg                                                                    >  5:125313/1‑147 (MQ=255)
                                                                            ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCTGTAACGACACCATGATCGCCAc                                                            >  5:366949/1‑147 (MQ=255)
                                                                            ccTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAc                                                            <  4:157137/147‑1 (MQ=255)
                                                                             cTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGt                                                                                                                       <  1:205666/90‑1 (MQ=255)
                                                                             cTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGt                                                                                                                       >  2:205666/1‑90 (MQ=255)
                                                                             cTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGcccc                                                            >  8:53685/1‑144 (MQ=255)
                                                                                 aGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGa                                                         <  3:356839/145‑1 (MQ=255)
                                                                                   ggcggcGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATgg                                                      >  3:80937/1‑146 (MQ=255)
                                                                                         gATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAt                                               >  8:238431/1‑147 (MQ=255)
                                                                                             gTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCTACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATgtgt                                           >  6:93492/1‑147 (MQ=255)
                                                                                             gTAACGACAAATTGTTCGGCGACGTGGGTACAGACCCCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCCGTAACGACACCATCATCGCCACCGATGGAAACGATgtgt                                           >  7:358998/1‑147 (MQ=255)
                                                                                                aCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGtt                                          >  4:373951/1‑145 (MQ=255)
                                                                                                  gagactttGGTCGGGGACGTGGGGACCGTCAGGGTGATGTGGGGGGCCGGCAACGACTTCCTCTACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGGGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATTTGTTCCAg                                      <  3:343366/141‑1 (MQ=255)
                                                                                                    cacaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGtt                                    >  3:258455/1‑147 (MQ=255)
                                                                                                           ttCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGc                              <  6:138733/146‑1 (MQ=255)
                                                                                                              ggCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGc                          >  2:379877/1‑147 (MQ=255)
                                                                                                                 gACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCttt                       <  5:172155/147‑1 (MQ=255)
                                                                                                                  aCGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTg                      >  4:237977/1‑147 (MQ=255)
                                                                                                                       ggTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAAc                 >  3:320248/1‑147 (MQ=255)
                                                                                                                       ggTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGCTCCAGTTCGCTGCAGGCTTTGGCaa                  >  8:81686/1‑146 (MQ=255)
                                                                                                                       ggTACCGACACCCTGATCAGGGGCGCCGGCAACGACTTCCTCAACCGTGCGGCAGGT‑‑AACGAAACGCTCGTGGCTGGCCCCGCTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCACTTCGCAGCAGGCATTGGaaa                  >  5:364348/1‑146 (MQ=255)
                                                                                                                               cacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGGGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGat         <  6:366949/147‑1 (MQ=255)
                                                                                                                                 cccTGATCAGGGGGGCCGGCAACGACTTCCTCATCGGTGCGGCAGGT‑‑AAGGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatca       <  6:125313/147‑1 (MQ=255)
                                                                                                                                  ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat      <  4:258455/147‑1 (MQ=255)
                                                                                                                                     gATCAGGGGGGTCGGCAACGACTTCCTCTACGCTGCGGCAGTT‑‑AACTATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCaa   <  8:358998/147‑1 (MQ=255)
                                                                                                                                      aTCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAt  >  6:322842/1‑147 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
TGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGT‑‑AACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921323‑2921603

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: