Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,889,221 A→G 9.5% L365P (CTG→CCG)  PP_1695 ← sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,889,2210AG 9.5% 109.9 / 3.1 42L365P (CTG→CCG) PP_1695sodium‑solute symporter/two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base A (20/18);  new base G (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGCCGATGGTCGCCAGGCTGGCGGTGGAACCCAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAGGGGCAGGCTGAT  >  NC_002947/1889075‑1889366
                                                                                                                                                  |                                                                                                                                                 
ggCCGATGGTCGCCAGGCTGGCGGTGGAACCCAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCgcagca                                                                                                                                                <  1:698186/150‑1 (MQ=255)
            ccAGGCTGGCGGTGGAACCCAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGgcagca                                                                                                                                    <  2:346077/150‑1 (MQ=255)
              aGGCTGGCGGTGGAACCCAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGgcagcagc                                                                                                                                  <  2:1609008/150‑1 (MQ=255)
                gCTGGCGGTGGAACCCAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCAt                                                                                                                                <  1:1430632/150‑1 (MQ=255)
                             cccAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACAcc                                                                                                                   >  2:1060682/1‑150 (MQ=255)
                             cccAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACAcc                                                                                                                   >  1:504460/1‑150 (MQ=255)
                             cccAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACAcc                                                                                                                   >  1:797541/1‑150 (MQ=255)
                              ccAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCTTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCa                                                                                                                  >  2:1323054/1‑150 (MQ=255)
                                agcagGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCGGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATg                                                                                                                >  1:886862/1‑150 (MQ=255)
                                agcagGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCGGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATg                                                                                                                >  1:886659/1‑150 (MQ=255)
                                 gcagGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATgg                                                                                                               <  1:1060682/150‑1 (MQ=255)
                                  cagGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGt                                                                                                              <  2:1489002/150‑1 (MQ=255)
                                     gCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCGGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGa                                                                                                           <  2:886862/150‑1 (MQ=255)
                                     gCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCGGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGa                                                                                                           <  2:886659/150‑1 (MQ=255)
                                              gACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCAc                                                                                                  <  2:797541/150‑1 (MQ=255)
                                                 aTAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGc                                                                                               <  1:1153711/150‑1 (MQ=255)
                                                 aTAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGc                                                                                               >  1:1638042/1‑150 (MQ=255)
                                                 aTAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGACAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGc                                                                                               >  1:1172336/1‑150 (MQ=255)
                                                   aGGCCAGCAGCAGGCAAACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCt                                                                                             >  2:1471077/1‑150 (MQ=255)
                                                       cagcagcagGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTcgac                                                                                         >  2:1554932/1‑150 (MQ=255)
                                                       cagcagcagGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAAAGCAACGGCCTcgac                                                                                         >  2:689715/1‑150 (MQ=255)
                                                                   caccacGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGAGAACGCCACGGCCTCGACGATCACGATGcc                                                                             >  2:777831/1‑150 (MQ=255)
                                                                     ccacGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCgg                                                                           >  1:1428261/1‑150 (MQ=255)
                                                                       acGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTa                                                                         <  2:1286338/150‑1 (MQ=255)
                                                                           tGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAgccg                                                                     >  2:1288017/1‑150 (MQ=255)
                                                                               tACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCTTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGa                                                                 <  1:1323054/150‑1 (MQ=255)
                                                                                  ccGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGc                                                              >  2:162921/1‑150 (MQ=255)
                                                                                          ccGAAAGCATCCAGTGGCGGCGCGCTGGGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCTTGTCGTTGGACCCCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGa                                                      <  2:1313353/150‑1 (MQ=255)
                                                                                               cgCATCCAGTGGCGGACCGATACACAAGGCCGCGCGGCTTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAAc                                                 <  2:1172336/149‑1 (MQ=255)
                                                                                                  aTCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTc                                                                                                             >  2:532182/1‑87 (MQ=255)
                                                                                                  aTCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTc                                                                                                             <  1:532182/87‑1 (MQ=255)
                                                                                                  aTCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGa                                                                                                           <  1:1647920/89‑1 (MQ=255)
                                                                                                  aTCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGa                                                                                                           >  2:1647920/1‑89 (MQ=255)
                                                                                                          gCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGc                                      >  2:361920/1‑150 (MQ=255)
                                                                                                                         aGGCCGCTCGGCGTTGTTGCGGAGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTgcgc                       <  1:285038/150‑1 (MQ=255)
                                                                                                                           gCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGAGCAAGGCTCGGGTGCGCCt                     >  2:584959/1‑150 (MQ=255)
                                                                                                                              gCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCgg                  >  2:1458601/1‑150 (MQ=255)
                                                                                                                              gCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCgg                  >  2:790003/1‑150 (MQ=255)
                                                                                                                                  ggCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAg              <  2:504460/150‑1 (MQ=255)
                                                                                                                                  ggCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCACCGAGGCACCGCAAGCGAGCGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAg              >  1:429535/1‑150 (MQ=255)
                                                                                                                                   gCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAgg             <  1:1288017/150‑1 (MQ=255)
                                                                                                                                        gttgCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAGGGGCAg        <  1:229728/150‑1 (MQ=255)
                                                                                                                                              gcgcAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAGGGGCAGGCtgat  <  2:1428261/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGCCGATGGTCGCCAGGCTGGCGGTGGAACCCAGCAGGCGGTAGCTGACATAGGCCAGCAGCAGGATCACCACGATGGTTACCCGGCGTACCGAAAGCATCCAGTGGCGGAACGCTTCGAAAGGCCGCTCGGCGTTGTTGCGGCGCAGCAGCCAGGGCAGCAGCATGTCGTTGGACACCATGGTCGACAACGCCACGGCCTCGACGATCACCATGCCGGTAGCCGCCGAGGCACCGCCGATGAACGCCAGCAGGGCAAGGCTCGGGTGCGCCTCGGCCAGGGGCAGGCTGAT  >  NC_002947/1889075‑1889366

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