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,4600GA90.5% 53.0 / ‑5.8 21G3494R (GGG→AGG) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (13/6);  minor base C (0/2);  total (13/8)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATT  >  NC_002947/2921320‑2921604
                                                                                                                                            |                                                                                                                                                
ggcggggggCGGCCGGTTACGACAACCTGTACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGc                                                                                                                                           <  6:76258/144‑1 (MQ=255)
              ggTCACGACCACCTGCTCGGCGGCGGGGGCGAGGACTGCCTGCATGGCGGCAACGGCAATGACAGCCTGGAGGGAGGCGATGGTAACGACACTTTGTGCGGCGACGTGGGTACCGACACCCTGATCAGGCGTGCCGGCAACGACTTcc                                                                                                                             <  6:125871/148‑1 (MQ=255)
                  aCGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCagg                                                                                                                                                <  5:98452/125‑1 (MQ=255)
                    gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACg                                                                                                                       >  2:60886/1‑148 (MQ=255)
                         ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCATGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCg                                                                                                                  >  3:140946/1‑148 (MQ=255)
                          cTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCg                                                                                                                  >  2:5056/1‑147 (MQ=255)
                                    gcgcGGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGGAACGACAGAGTGGTCGGCGACGGCGGTACCGACACCCTGATCCGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGa                                                                                                       <  1:7963/148‑1 (MQ=255)
                                          gCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACAGTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCAGGCAACGACTTCCTCAACGGGGCGGCAGGTAACGAGACggt                                                                                                 >  3:41030/1‑148 (MQ=255)
                                               gACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTGCGGCGACGTGTGTACCGACACCCTGATCAGGGTGGCCGGCAACGACTTCCTCAACGGTGAGGCAGGTACCGATTCGGtggggg                                                                                            >  6:153019/1‑148 (MQ=255)
                                                    ccTGGATGGCGGCAACGGCAATGACACCATGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGtgggtggg                                                                                       <  4:140946/148‑1 (MQ=255)
                                                     cTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTgggg                                                                                      >  7:135075/1‑148 (MQ=255)
                                                             cggcaccgccaaTGCCACCCTGCAGGGCGGCGACGTGACCGACCCGTTGTGCGGCGCCGTGGGTCCCGACACCCTGATCCGGGGGGGCGGCAACGACTTCCTCAACGGGGGGGCAGGTAACGATACGGTGGTGGGGGGGGCCGGTAAc                                                                              <  2:104684/148‑1 (MQ=255)
                                                                    ggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCATCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGTGCCGGTAACGacaaca                                                                       >  8:128610/1‑148 (MQ=255)
                                                                            cacCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGGGGGGACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGTTAACGATACGGTGGTGGGGGGGGCCGGTAACGACACCATGATCGc                                                                <  7:128610/147‑1 (MQ=255)
                                                                                    tgggCGGCGATGGTAACGACACTTTGTTCGGCGACGTGGGTCCCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGAt                                                         <  8:135075/145‑1 (MQ=255)
                                                                                                    cGACACATTGTTCGGCGACGTGGGTACCGACACCCTGACCAGGGGGGCCGGCAACGACTTCCTCAACGGGGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                        <  1:100956/147‑1 (MQ=255)
                                                                                                        acaTTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCg                                   >  6:105063/1‑148 (MQ=255)
                                                                                                                                 acacCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGa           >  4:28164/1‑147 (MQ=255)
                                                                                                                                     ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat       >  1:29822/1‑147 (MQ=255)
                                                                                                                                      cTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGAtcatca     >  1:42750/1‑148 (MQ=255)
                                                                                                                                       tGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCaa    >  2:161493/1‑148 (MQ=255)
                                                                                                                                         aTCAGGGGGGCCAGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAtt  >  1:116851/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                
TGCTGGGGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAATT  >  NC_002947/2921320‑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: