Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,587,208 A→G 22.6% G23G (GGT→GGC mdtF ← multidrug transporter, RpoS‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,587,2080AG22.6% 157.0 / 38.0 119G23G (GGT→GGCmdtFmultidrug transporter, RpoS‑dependent
Reads supporting (aligned to +/- strand):  ref base A (16/76);  new base G (10/17);  total (26/93)
Fisher's exact test for biased strand distribution p-value = 3.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

GTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATA  >  minE/2587146‑2587267
                                                              |                                                           
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:1052495/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:816955/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:2154457/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:1699072/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:2396339/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:163359/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:1756781/69‑1 (MQ=255)
gTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCa                                                       <  1:434532/69‑1 (MQ=255)
  aaTGTTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:1692398/69‑1 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCaa                                                      >  1:3169824/1‑68 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:2552770/69‑1 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:3338952/1‑69 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:1756412/1‑69 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:2562701/1‑69 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:3336128/69‑1 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:1832801/69‑1 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:2452526/69‑1 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:2110578/1‑69 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:2116682/1‑69 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:2376575/1‑69 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:1437098/69‑1 (MQ=255)
  aaTGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     >  1:513460/1‑69 (MQ=255)
     ggTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGcat                                                  <  1:1481198/69‑1 (MQ=255)
      gTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGcat                                                  >  1:2526830/1‑68 (MQ=255)
      gTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGcat                                                  <  1:732773/68‑1 (MQ=255)
      gTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGcat                                                  >  1:1188125/1‑68 (MQ=255)
      gTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGcat                                                  >  1:880318/1‑68 (MQ=255)
           tGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:2986101/69‑1 (MQ=255)
           tGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:1634151/69‑1 (MQ=255)
           tGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:1649179/69‑1 (MQ=255)
           tGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:580150/69‑1 (MQ=255)
           tGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:1825029/69‑1 (MQ=255)
           tGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:334464/69‑1 (MQ=255)
            ggCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            <  1:1800215/68‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCATACCACCTGCAAg                                                     <  1:755613/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:641362/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:2564471/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:511297/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:2880703/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:3409925/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:2902838/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:3409359/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:1027396/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:1373008/50‑1 (MQ=255)
                     tGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:1460293/50‑1 (MQ=255)
                         gATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATCATAATGGCAAGTAc                                <  1:3195246/67‑1 (MQ=255)
                         gATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTa                                 <  1:2852178/66‑1 (MQ=255)
                         gATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTAc                                <  1:2278354/67‑1 (MQ=255)
                         gATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTAc                                <  1:1826818/67‑1 (MQ=255)
                             cTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAgg                           >  1:1852440/1‑68 (MQ=255)
                              tGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAg                                                     <  1:262651/41‑1 (MQ=255)
                               gcgcAACCGGTAAGTTCATGATCGCCAGACCACCTGCAACATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:2154769/69‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaaa                      <  1:2754178/68‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaaa                      <  1:2992573/68‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaaa                      <  1:254436/68‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaaa                      <  1:855295/68‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:430933/67‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:147311/67‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:1770969/67‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:923774/67‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:3176617/67‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:20979/67‑1 (MQ=255)
                                  cAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:2591188/67‑1 (MQ=255)
                                       ggTAAGTTCATGATCGCCAGACCACCTGCAAGCATCataat                                            >  1:502742/1‑41 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:345685/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:1733845/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:2459207/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:1008256/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:1638263/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:3106149/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:1708900/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:1214469/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCg              >  1:1437078/1‑68 (MQ=255)
                                          aaGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:2725146/59‑1 (MQ=255)
                                          aaGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:1805873/59‑1 (MQ=255)
                                          aaGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCaaa                       <  1:2735490/59‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:2105576/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:1379362/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:58615/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:79025/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:1099690/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:89422/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:2287017/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:458992/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:1820436/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:3336251/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:1686299/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:3181973/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:2996399/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:1729332/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:2960998/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAAcc                   <  1:2706523/56‑1 (MQ=255)
                                                 tGATCGCCAGACCGCCTGCAAGCATCATAATAATGGCAAGTA‑‑CAGGCAAAAAcc                   <  1:1475236/54‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:1530378/68‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:1551633/68‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:1660912/68‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:661810/68‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:3283537/68‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:2593759/68‑1 (MQ=255)
                                                  gATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATaa      <  1:2296039/68‑1 (MQ=255)
                                                              |                                                           
GTAATGGTCGGTGGCGCAATCTGCGGATACTGCGCAACCGGTAAGTTCATGATCGCCAGACCACCTGCAAGCATCATAATAATGGCAAGTACCCAGGCAAAAACCGGGCGATCAATAAAATA  >  minE/2587146‑2587267

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