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

GGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATC  >  NC_002947/2921326‑2921600
                                                                                                                                      |                                                                                                                                            
gggAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAAc                                                                                                                                <  4:306484/149‑1 (MQ=255)
         gTAACGACAACCTGAACGGGGGCGCTGGCGATGAAAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTACGGCGACCTGGGAACCGACACCATGAACAGGGGGGCCGCCAACGACtaact                                                                                                                        >  1:210813/1‑144 (MQ=255)
            aCGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGTGCCGGCACCGACTtcatcacc                                                                                                                    >  4:245785/1‑147 (MQ=255)
              gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGTGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACgg                                                                                                                  >  7:43693/1‑149 (MQ=255)
              gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACgg                                                                                                                  >  7:95063/1‑149 (MQ=255)
              gACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACgg                                                                                                                  >  3:56121/1‑149 (MQ=255)
                   ccTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCgg                                                                                                             >  4:111190/1‑149 (MQ=255)
                        aagggggCGCTGGAGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGAAGTGCGTACCGACACCCTGATCAGGGCGGCCGGCAACGACTTCATCAACGGTGCGTCAgc                                                                                                         >  8:161619/3‑147 (MQ=255)
                          gggtgCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTGCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTaa                                                                                                      <  2:240766/149‑1 (MQ=255)
                                 cTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGc                                                                                                               <  6:55140/133‑1 (MQ=255)
                                 cTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGc                                                                                                               >  5:55140/1‑133 (MQ=255)
                                        tGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGtggtgg                                                                                        >  4:249159/1‑149 (MQ=255)
                                         gACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTggtggg                                                                                       >  7:214466/1‑149 (MQ=255)
                                             gCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGtgggtggg                                                                                   <  3:111190/149‑1 (MQ=255)
                                              ccTGGATGGCGGCAACGGCAATGACACCCTGCCGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCCGCAGGTAACGATACGGTGGTGGGTgggg                                                                                  <  1:181844/149‑1 (MQ=255)
                                                   aTGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTc                                                                                                                       <  1:187259/107‑1 (MQ=255)
                                                   aTGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTc                                                                                                                       >  2:187259/1‑107 (MQ=255)
                                                        ggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCg                                                                                                              >  5:19137/1‑111 (MQ=255)
                                                        ggcaacggcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCg                                                                                                              <  6:19137/111‑1 (MQ=255)
                                                               gcaaTGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATg                                                                 <  4:56121/149‑1 (MQ=255)
                                                                     acacCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGc                                                            <  3:245785/148‑1 (MQ=255)
                                                                              aGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGa                                                  <  3:249159/149‑1 (MQ=255)
                                                                                   ggcgATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGa                                             >  4:135180/1‑149 (MQ=255)
                                                                                   ggcgATGGTAACGACACATTGTTCGGCGACGTGGGCACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGa                                             >  8:255637/1‑149 (MQ=255)
                                                                                    gcgATGGTACCGACACATTGTTCGGCGACGTGGGGACCGACACCCGGATCCTGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGAt                                            <  7:299838/149‑1 (MQ=255)
                                                                                          gTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTc                                      <  8:214466/149‑1 (MQ=255)
                                                                                             aCGACACATTGTTCGGTGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                    <  8:43693/148‑1 (MQ=255)
                                                                                             aCGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCa                                    <  8:95063/148‑1 (MQ=255)
                                                                                              cGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGc                                                            >  7:256616/1‑123 (MQ=255)
                                                                                              cGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGc                                                            <  8:256616/123‑1 (MQ=255)
                                                                                                       gTTCGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAg                         >  7:205056/1‑149 (MQ=255)
                                                                                                          cGGCGACGTGGGTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAAc                                                                          <  8:189998/97‑1 (MQ=255)
                                                                                                                     gTACCGACACCCTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGAc           >  7:107570/1‑149 (MQ=255)
                                                                                                                              cccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcat   >  5:119632/1‑148 (MQ=255)
                                                                                                                               ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc  >  4:147752/1‑148 (MQ=255)
                                                                                                                               ccTGATCAGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGatcatc  >  8:610/1‑148 (MQ=255)
                                                                                                                                      |                                                                                                                                            
GGGAGGCCGGTAACGACAACCTGAACGGGGGCGCTGGCGATGACAGCCTGGATGGCGGCAACGGCAATGACACCCTGCAGGGCGGCGATGGTAACGACACATTGTTCGGCGACGTGGGTACCGACACCCTGATCGGGGGGGCCGGCAACGACTTCCTCAACGGTGCGGCAGGTAACGATACGGTGGTGGGTGGGGCCGGTAACGACACCATGATCGCCACCGATGGAAACGATGTGTTCCAGTTCGCTGCAGGCTTTGGCAACGACCTGATCATC  >  NC_002947/2921326‑2921600

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: