Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,216,188 A→G 15.3% F97L (TTT→CTT)  yehW ← glycine betaine ABC transporter membrane subunit YehW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,216,1880AG15.3% 64.1 / 4.2 26F97L (TTT→CTT) yehWglycine betaine ABC transporter membrane subunit YehW
Reads supporting (aligned to +/- strand):  ref base A (11/11);  new base G (2/2);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

TTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGCAGTTGTGGGAACAAGGCAGCAAACAGCGGCTGGCTGTAAGGCAGCCAGAAAATCAGCGCCACAAACAGAGCA  >  NC_000913/2216020‑2216441
                                                                                                                                                                        |                                                                                                                                                                                                                                                             
ttCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGCCGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAACTTCCGCGCCCCCCCGGCGAGTGACAGACAATCACACACCAGTCACAATGacacacccaaaaaa                                                                                                                                                      >  2:291324/1‑261 (MQ=255)
                                                      cAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAAt                                                                                                                                                                                                                    >  2:217296/1‑158 (MQ=255)
                                                      cAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTACCACCAGCCAGAAATGAACCAGTGCCAGAGCTGCAAAACTTTc                                                                                                >  1:106572/1‑274 (MQ=255)
                                                         gtgtCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTcc                                                                                                                                                                                                                 <  1:217296/158‑1 (MQ=255)
                                                                                                      tAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGCAg                                                                        >  1:99817/1‑250 (MQ=255)
                                                                                                          gtcggggtGCAGACCAAATCCGTTCACCGGACTGGCGTTCCCAAGCACTGGGACGGACGGAAAAGCCGGACCACCGGCGGCAATAGTTTCCACCAGTGGACGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGCAGTTGTGGGAACAAGGCAGCAAACAGCggc                                            <  2:106572/267‑1 (MQ=255)
                                                                                                          cGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGCAGTTGTGGGAACAAGGCAGCAAACAGCgg                                             <  1:258138/273‑1 (MQ=255)
                                                                                                           gCTGGTGGCCGACCAAGGCCGATCCCCGGAACGGCGTGCGCCAGCACGGCGAGGGGCGGAAAAGTCTGTCCAACGGGGGCAAGAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATCCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGGGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGCAGTTGTGGGAACAAGGCAGCAAACAGCggc                                            <  1:291324/273‑1 (MQ=255)
                                                                                                                            gCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAGAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATg                                                                                                                           >  1:281540/1‑177 (MQ=255)
                                                                                                                            gCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCg                                                                                                                                                                                   >  2:99099/1‑121 (MQ=255)
                                                                                                                               gATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAGAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGc                                                                                                                        <  2:281540/177‑1 (MQ=255)
                                                                                                                               gATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCa                                                                                                                                                                                <  1:99099/121‑1 (MQ=255)
                                                                                                                                 tCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCtt                                                                                              >  1:79847/1‑201 (MQ=255)
                                                                                                                                      ggAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCg                                                                                                                                                                                   >  2:112186/1‑111 (MQ=255)
                                                                                                                                         aCGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCa                                                                                                                                                                                <  1:112186/111‑1 (MQ=255)
                                                                                                                                              gATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCa                                                                                                                                                                                                         >  1:248845/1‑81 (MQ=255)
                                                                                                                                                 cGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACgg                                                                                                                                                                                                      <  2:248845/81‑1 (MQ=255)
                                                                                                                                                     aGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGGGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAGCGGGTCGTGGCAGTTGTGGGAACAAGGCAGCAAACAGCGGCTGGCTGTAAGGCAGCCAGAAAATCAGCGCCACAAACAGAGCa  <  1:158101/273‑1 (MQ=255)
                                                                                                                                                             gACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCcgcg                                                                                                                                                                                                              >  1:25574/1‑61 (MQ=255)
                                                                                                                                                                gggCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGccc                                                                                                                                                                                                           <  2:25574/61‑1 (MQ=255)
                                                                                                                                                                  gCGGAAGAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAAt                                                                                                                            >  2:101316/1‑138 (MQ=255)
                                                                                                                                                                  gCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATg                                                                                                                           >  1:221939/1‑139 (MQ=255)
                                                                                                                                                                  gCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGt                                                                             >  1:283273/1‑185 (MQ=255)
                                                                                                                                                                     gAAGAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAg                                                                                                                         <  1:101316/138‑1 (MQ=255)
                                                                                                                                                                     gAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGc                                                                                                                        <  2:221939/139‑1 (MQ=255)
                                                                                                                                                                     gAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGc                                                                          <  2:283273/185‑1 (MQ=255)
                                                                                                                                                                        |                                                                                                                                                                                                                                                             
TTCCCATACCTTTCGCAACTTCTGTCACGCTGGCATCAATCGCTCCCAGCCCGGCAAGTGTCGCCTGCAGGACGGGCAGCACACCGTAAAGGATCAAGGCGATAATCGCTGGTTGCAGACCAAAGCCGATCACCGGAACGGCGATCGCCAGCACTGCGACGGGCGGAAAAGTCTGTCCAACGGCGGCAATAGTTTCCACCAGTGGGCGAAATTCCGCGCCCCACGGGCGAGTGACAGCAATTCCGGCACCAGTGCCAATGATCACCGCAAACAAACTCGAAATTCCCACCAGCCAGAAATGAGCCAGTGCCAGAGCTGCAAAACTTTCTTGCTGATAAACGGGTCGTGGCAGTTGTGGGAACAAGGCAGCAAACAGCGGCTGGCTGTAAGGCAGCCAGAAAATCAGCGCCACAAACAGAGCA  >  NC_000913/2216020‑2216441

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑