Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,910,534 G→A 9.3% intergenic (‑54/‑490) aprDA ← / → PP_2561 alkaline protease secretion ABC transporter ATP‑binding protein/hemolysin‑type calcium‑binding bacteriocin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,910,5340GA 9.3% 112.3 / 2.4 43intergenic (‑54/‑490)aprDA/PP_2561alkaline protease secretion ABC transporter ATP‑binding protein/hemolysin‑type calcium‑binding bacteriocin
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CGCTGAAGGCGGCTAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATACTTTGAAC  >  NC_002947/2910393‑2910676
                                                                                                                                             |                                                                                                                                              
cgcTGAAGGCGGCTAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCtgttgt                                                                                                                                        <  1:566106/150‑1 (MQ=255)
cgcTGAAGGCGGCTAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCtgttgt                                                                                                                                        >  1:442630/1‑150 (MQ=255)
     aaGGCGGCTAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCtgtt                                                                                                                                          >  2:866900/1‑143 (MQ=255)
     aaGGCGGCTAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCtgtt                                                                                                                                          <  1:866900/143‑1 (MQ=255)
             tAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCg                                                                                                                           >  2:191668/1‑150 (MQ=255)
                   gTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGc                                                                                                                     >  2:864596/1‑150 (MQ=255)
                   gTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACAGCcgttgtacgg                                                                                                                                    >  2:1205936/1‑131 (MQ=255)
                   gTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACAGCcgttgtacgg                                                                                                                                    <  1:1205936/135‑5 (MQ=255)
                   gTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACAGCcgttgtacgg                                                                                                                                    <  1:1206003/135‑5 (MQ=255)
                   gTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACAGCcgttgtacgg                                                                                                                                    >  2:1206003/1‑131 (MQ=255)
                      gAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAgg                                                                                                                  >  2:64162/1‑150 (MQ=255)
                         gCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCATTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCgg                                                                                                               >  1:1087139/1‑150 (MQ=255)
                               gcgcAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCa                                                                                                         <  1:1099520/150‑1 (MQ=255)
                                 gcAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGAGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCCCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACt                                                                                                       >  1:1223662/1‑150 (MQ=255)
                                    aGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCATTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGAt                                                                                                    <  2:1087139/150‑1 (MQ=255)
                                      tGAACCAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATgg                                                                                                  <  1:191668/150‑1 (MQ=255)
                                      tGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGGGGCCCAGTCGAGAGGCGCCCCGTACAACGGATGTTCTACCACTTTCACCGCCGCGCAGGCGTGCCGCACTGATgg                                                                                                  >  1:988307/1‑150 (MQ=255)
                                               aCGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGc                                                                                         <  1:64162/150‑1 (MQ=255)
                                                 gTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGccc                                                                                       >  2:1126316/1‑150 (MQ=255)
                                                    ggCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCACTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCaaa                                                                                    <  2:684582/150‑1 (MQ=255)
                                                          ttCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGc                                                                              <  2:442630/150‑1 (MQ=255)
                                                             ggAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTGGTACTGATGGCGATATTGCCCAAACGCTGCGtt                                                                           <  1:864596/150‑1 (MQ=255)
                                                                           gtggCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCt                                                                                                                                             >  1:1184256/1‑70 (MQ=255)
                                                                           gtggCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCt                                                                                                                                             <  2:1184256/70‑1 (MQ=255)
                                                                              gCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGAcagcag                                                          >  1:521899/1‑150 (MQ=255)
                                                                                ttAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTc                                                                       <  1:1058357/135‑1 (MQ=255)
                                                                                ttAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTc                                                                       >  2:1058357/1‑135 (MQ=255)
                                                                                ttAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCg                                                        >  1:1027691/1‑150 (MQ=255)
                                                                                    gCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAgg                                                    <  1:923793/150‑1 (MQ=255)
                                                                                          gAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCaa                                              <  1:464363/150‑1 (MQ=255)
                                                                                               cTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGaaa                                         <  2:521899/150‑1 (MQ=255)
                                                                                                      tATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAg                                  <  1:1202053/150‑1 (MQ=255)
                                                                                                        tCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGc                                                                              >  2:1206184/1‑104 (MQ=255)
                                                                                                        tCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGc                                                                              <  1:1206184/104‑1 (MQ=255)
                                                                                                          acaAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAat                              >  2:8337/1‑150 (MQ=255)
                                                                                                                    aaGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGc                    <  1:1100364/150‑1 (MQ=255)
                                                                                                                               gcgcCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATAc         <  2:1027691/150‑1 (MQ=255)
                                                                                                                                  ccccGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATACttt      >  2:987680/1‑150 (MQ=255)
                                                                                                                                  ccccGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATACttt      <  1:8337/150‑1 (MQ=255)
                                                                                                                                    ccGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATACTTTGa    >  2:416419/1‑150 (MQ=255)
                                                                                                                                      gTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATACTTTGAAc  >  1:646589/1‑150 (MQ=255)
                                                                                                                                        aCAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTa                                                      >  1:239241/1‑96 (MQ=255)
                                                                                                                                        aCAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTa                                                      <  2:239241/96‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
CGCTGAAGGCGGCTAGGGCGTAGAAGCTGCGGCGCAAGTGAAACAGCACGTCGGCCAGTTCGGAGCGGGTGCGTGGTGGCTTAGGCATCGGAGGCCTCCCTTTATCACAAAGCGGGAAGCCCTTGCGGCGCCCCGTACAACGGCTGTTGTGGATTTTTGGTCGGCTCGCAGGCGGGTCGCACTGATGGCGATATTGCCCAAACGCTGCGTTGCTCCCACAGACAGCAGCGTAGGTTGCAACGAAACATCCAGTAATATCATTTCGCTCTATAAATACTTTGAAC  >  NC_002947/2910393‑2910676

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