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

TTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921319‑2921603
                                                                                                                                             |                                                                                                                                               
ttGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGc                                                                                                                                          <  5:381083/149‑1 (MQ=255)
  gCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCaggg                                                                                                                                              <  5:189803/143‑1 (MQ=255)
  gCTGGGGGAGGCCGGTAACGACAACCTGAACGGGCGCGCTGGCGATGACATCCAGGATGGCGGCAACGGCAATGACACCATGCAGTGCGGAGATGGTAACGACAAATTTTTCGGCGACCTGGGTACCGACACCCTGATCaggg                                                                                                                                              >  6:189803/1‑143 (MQ=255)
            cccgccaaCGACAACCCGACCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACtt                                                                                                                              <  3:160923/143‑1 (MQ=255)
               ggTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTcc                                                                                                                            >  4:186972/1‑148 (MQ=255)
               ggTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTcc                                                                                                                            <  3:186972/148‑1 (MQ=255)
                          ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCg                                                                                                                                        <  6:170747/125‑1 (MQ=255)
                          ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCg                                                                                                                                        >  5:170747/1‑125 (MQ=255)
                                   gggCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCACCGACTTCCTCAACGGTGCGGCAGGTAAc                                                                                                        >  5:16855/1‑148 (MQ=255)
                                       gcTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTc                                                                                                                             <  4:170099/123‑1 (MQ=255)
                                       gcTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTc                                                                                                                             >  3:170099/1‑123 (MQ=255)
                                                          ctGGCGGGACCGGCAATGTCACCCGGCAGGGCGGCGAGGGAACCGACACATTGTTCGGCGACGGGGGCACCGACACCCTCATCAGGGGGGCCGGCAACGACTTCCTCAACGGGGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGt                                                                                <  4:6113/148‑1 (MQ=255)
                                                                      gcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCAt                                                                     >  5:149784/1‑148 (MQ=255)
                                                                                ccTGCAGGGCGGCGATGGTAACGACACATTGTACGGCGACGTGGCTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCAcc                                                           >  1:12930/1‑148 (MQ=255)
                                                                                                 gTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGCGGCAGGTAACGACACCATGATCGCAACAGATGGAAACGATGTCTTc                                         >  6:205864/1‑149 (MQ=255)
                                                                                                 gTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTCGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                         >  5:76085/1‑149 (MQ=255)
                                                                                                  tAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg                                                                                                             <  4:151553/80‑1 (MQ=255)
                                                                                                  tAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAg                                                                                                             >  3:151553/1‑80 (MQ=255)
                                                                                                              gTTCGGCGATTTTGATCCCCACACCCTGATCTGGTGGGCGTGCTACGACTTCGACATCTGTTCGGCTGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCa                             <  1:336150/148‑1 (MQ=255)
                                                                                                                             tACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGgtggtg                                                                                            >  2:145915/1‑70 (MQ=255)
                                                                                                                                ccacacCCTGATCACGGGGCCCGGCACCGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGCTGGGTGCGGCCGGTAACGACACCATGATCGCCACCGCAGGAAACGATGTGATCCAGTTCGATGAAGGATTTGGCAACGACCTg           >  2:262254/3‑148 (MQ=255)
                                                                                                                                 gACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATgg                                                      <  2:41888/104‑1 (MQ=255)
                                                                                                                                 gACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATgg                                                      >  1:41888/1‑104 (MQ=255)
                                                                                                                                 gACACCCTGATCAGCGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATAAGGTGCTGGGTGGGGCCGTTAACGACACCATGATCGCCACAGATGTAAACGATGTCTTCCAGTTCGCTGCAGGCTTTGGCAACGACCtca          >  1:232904/1‑146 (MQ=255)
                                                                                                                                  acacCGTGATCAGTGGGGCCGGCAACGACTTCCTCTACGGTGCGGCAGGTAACGATACGGTGGTGGGTGCGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatc        <  1:162076/149‑1 (MQ=255)
                                                                                                                                  acacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTTGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatc        <  6:76085/149‑1 (MQ=255)
                                                                                                                                      ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc     >  1:26755/1‑148 (MQ=255)
                                                                                                                                      ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc     >  7:216317/1‑148 (MQ=255)
                                                                                                                                      ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGAtcatca    <  2:26755/149‑1 (MQ=255)
                                                                                                                                       cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                         >  7:327729/1‑111 (MQ=255)
                                                                                                                                       cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                         <  8:327729/111‑1 (MQ=255)
                                                                                                                                          aTCAGGGGGTCCGGCAACGACTTCCTCTACGGTGCGTCAGGTAACGATACGGTGGTGGGTGGGGCCTGTAACGACACCATGATCTCCACCGATGGAAAAGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAt  <  5:205864/147‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
TTGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921319‑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: ‑