Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,916,227 C→T 13.3% N1749N (AAC→AAT)  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,916,2270CT13.3% 71.2 / 3.6 30N1749N (AAC→AAT) PP_2561hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base C (17/9);  new base T (2/2);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GACCTGGTCACCAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGTTGATCACCA  >  NC_002947/2916084‑2916375
                                                                                                                                               |                                                                                                                                                    
gACCTGGTCACCAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCcga                                                                                                                                                 >  3:246106/1‑149 (MQ=255)
      gTCACCAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCcga                                                                                                                                                 >  8:241756/1‑143 (MQ=255)
      gTCACCAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCcga                                                                                                                                                 <  7:241756/143‑1 (MQ=255)
         aCCAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAATCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCcg                                                                                                                                        <  4:314408/149‑1 (MQ=255)
          ccAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCcgc                                                                                                                                       <  4:246106/149‑1 (MQ=255)
             aGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGaa                                                                                                                                    <  2:141221/149‑1 (MQ=255)
                          tACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACacaccaac                                                                                                                           >  4:302062/1‑145 (MQ=255)
                          tACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACacaccaac                                                                                                                           <  3:302062/145‑1 (MQ=255)
                             cATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGACCCCACCAAcaccaccc                                                                                                                    >  4:42149/1‑148 (MQ=255)
                                  aCATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACCCACCAACACCACCACGCCg                                                                                                               >  6:169565/1‑149 (MQ=255)
                                              aaCCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAATGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGac                                                                                                      >  3:71177/1‑146 (MQ=255)
                                              aaCCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACTAACACCACCACGCCGTTCGTCGac                                                                                                      >  8:258880/1‑146 (MQ=255)
                                                                 cGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCg                                                                                 <  1:61193/148‑1 (MQ=255)
                                                                                 aCCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGgt                                                                >  7:45352/1‑149 (MQ=255)
                                                                                     aCTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACTAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTc                                                            <  7:258880/149‑1 (MQ=255)
                                                                                                                aCCTGCCTGGGCCGGACGGCATCCTCGGCAACGCTGACGACATCCGCGAACACACCAACACCACCACGCCGATCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGAACAAGGTGTTCCTGCGCGCCTATGTCATGACTgacgac                                 >  7:228730/1‑149 (MQ=255)
                                                                                                                     cctgGGCCGGACGGCATCCTCGGCAATGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGcc                             <  4:71177/148‑1 (MQ=255)
                                                                                                                                cGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACgg                 >  8:224748/1‑149 (MQ=255)
                                                                                                                                 ggCATCCTCGGCAATGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCAt                                           >  4:263404/1‑122 (MQ=255)
                                                                                                                                 ggCATCCTCGGCAATGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCAt                                           <  3:263404/122‑1 (MQ=255)
                                                                                                                                  gCATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGc               >  4:169597/1‑149 (MQ=255)
                                                                                                                                   cATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGcc              >  4:230772/1‑149 (MQ=255)
                                                                                                                                    aTCCTCGGCAACGCCGACGACATACGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCAGGACTGACGACGGCCCGGTCGCCACGGGCCg             >  7:205768/1‑149 (MQ=255)
                                                                                                                                      ccTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGt           >  4:250723/1‑149 (MQ=255)
                                                                                                                                       cTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGtt          <  7:25627/149‑1 (MQ=255)
                                                                                                                                       cTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGtt          >  7:83418/1‑149 (MQ=255)
                                                                                                                                        tCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGtt          >  1:26573/1‑148 (MQ=255)
                                                                                                                                          ggCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACGAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTc                                             <  2:150328/111‑1 (MQ=255)
                                                                                                                                          ggCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTc                                             >  1:150328/1‑111 (MQ=255)
                                                                                                                                              aCGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGTTGATCAc    >  3:116880/1‑148 (MQ=255)
                                                                                                                                               cGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGTTGATCACCa  <  1:125515/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GACCTGGTCACCAAGGGCGGCTCCGGTACCATCTACATCCCGCTGCAACCGGATGACCCGCTGTACGTGGAGGGTGGCTTCACCAACTTCATGGTGGTAACTCGTGCCACCAACCTGCCTGGGCCGGACGGCATCCTCGGCAACGCCGACGACATCCGCGAACACACCAACACCACCACGCCGTTCGTCGACCAGAACCAGACCTACAGCTCGCATCCCTCGCACCAGGTGTTCCTGCGCGCCTATGTCATGACTGACGACGGCCCGGTCGCCACGGGCCGGTTGATCACCA  >  NC_002947/2916084‑2916375

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑