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,4600GA81.2% 32.8 / ‑4.6 17G3494R (GGG→AGG) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/4);  minor base C (0/3);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 6.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921334‑2921603
                                                                                                                              |                                                                                                                                               
ggTACCGACACCCTGACCGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAATGACTTcc                                                                                                                            <  2:207946/148‑1 (MQ=255)
                              gATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACggtgg                                                                                              >  8:121333/1‑148 (MQ=255)
                                tGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACAGTGTTCGGCGACGTGGGTACCGACACCCTGATCCGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGgtggtg                                                                                            <  3:26987/148‑1 (MQ=255)
                                    aGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTcgtgggtg                                                                                        >  7:111167/1‑148 (MQ=255)
                                            tGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGt                                                                                <  7:121333/148‑1 (MQ=255)
                                              gcggcAACGGCAATGCCACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGGAACTCTACGGTGGTGGGTGGGGCAGGTaa                                                                              >  7:98190/1‑148 (MQ=255)
                                                           tGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCTGTAACGACACCATGATc                                                                 <  8:111167/148‑1 (MQ=255)
                                                                                       gACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACAACATGATCGCCACAGATGGAAACGATGTGTTCCAGt                                     >  2:156180/1‑148 (MQ=255)
                                                                                            aTTGTACGGCGCCGTGTGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGCTGGGTGCCGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGAGTTCCAGTTCGCt                                >  6:98567/1‑148 (MQ=255)
                                                                                                ttCGGCGACGTGGGTACCGACACCCTGATCCGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCa                             <  5:82104/147‑1 (MQ=255)
                                                                                                     cGACGTGGGTACCGACACCCTGATCCGGGGGGCGGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCttt                       <  4:123227/148‑1 (MQ=255)
                                                                                                      tccGTGGGGCCCGCCACCCTGATCTGGGGGGGGGGCAACGACTTCCTCTACGGTGCGGCAGGTAACGATACGGTGGGGGGGGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTg                      <  6:62070/146‑1 (MQ=255)
                                                                                                                 cGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACTGTGCGGCACGTAACGATACGGTGGTGGGTGGGTCAGGTAACGCCACCATGATCGCCACCGATGGAAACGa                                                >  6:135081/1‑111 (MQ=255)
                                                                                                                  gACACCCGGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACTAt                                               <  5:135081/111‑1 (MQ=255)
                                                                                                                      cccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat      >  4:13843/1‑148 (MQ=255)
                                                                                                                       ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGACGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat      >  7:112088/1‑147 (MQ=255)
                                                                                                                          gATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAt  >  2:10792/1‑148 (MQ=255)
                                                                                                                              |                                                                                                                                               
GGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATCAAT  >  NC_002947/2921334‑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: