Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 368,298 A→G 47.1% Q196R (CAG→CGG)  fes → enterobactin/ferric enterobactin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE368,2980AG47.1% 4.0 / 74.1 85Q196R (CAG→CGG) fesenterobactin/ferric enterobactin esterase
Reads supporting (aligned to +/- strand):  ref base A (10/35);  new base G (24/16);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 7.57e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

CGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAC  >  minE/368233‑368363
                                                                 |                                                                 
cGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCa                                                                <  1:2534407/69‑1 (MQ=255)
cGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCa                                                                <  1:2388570/69‑1 (MQ=255)
cGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCa                                                                <  1:1515564/69‑1 (MQ=255)
cGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCa                                                                <  1:1627810/69‑1 (MQ=255)
 gATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCa                                                                <  1:1319481/68‑1 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:2832946/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:601425/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:2610701/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:2668660/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:774876/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:59905/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:542407/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:2721580/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:2960252/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:1713084/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:1649051/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:1339338/1‑68 (MQ=255)
       gAATTTTGGGCGCAGAGTATGCCCGTCTGGCCTGCGCTGACTTCGCTGACCCATCGTCGGCAACTtccg                                                         >  1:967674/1‑68 (MQ=255)
                 cgcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCgtg                                                <  1:2121258/68‑1 (MQ=255)
                  gcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTa                                              <  1:52024/69‑1 (MQ=255)
                  gcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTa                                              <  1:2748389/69‑1 (MQ=255)
                  gcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTa                                              <  1:840204/69‑1 (MQ=255)
                  gcAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTa                                              <  1:3273989/69‑1 (MQ=255)
                         aTGCCCGTCTTGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2099350/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:250146/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2695411/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2972953/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:3240526/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:342638/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:3354972/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2443152/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2375731/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2313638/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:2183376/69‑1 (MQ=255)
                         aTGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTg                                       <  1:20400/69‑1 (MQ=255)
                          tGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGa                                      >  1:1195527/1‑69 (MQ=255)
                          tGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGa                                      >  1:3404640/1‑69 (MQ=255)
                          tGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGa                                      >  1:599304/1‑69 (MQ=255)
                                  tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCt                              <  1:668986/69‑1 (MQ=255)
                                  tGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCt                              <  1:133834/69‑1 (MQ=255)
                                   ggCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCt                              <  1:2528138/68‑1 (MQ=255)
                                      cAGTGCTGACTTCGCTGAGCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATc                           <  1:1210909/68‑1 (MQ=255)
                                      cAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATc                           <  1:248475/68‑1 (MQ=255)
                                      cAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATc                           <  1:3044544/68‑1 (MQ=255)
                                      cAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATc                           <  1:1816131/68‑1 (MQ=255)
                                      cAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATc                           <  1:2261867/68‑1 (MQ=255)
                                          gCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGacac                      <  1:2066387/69‑1 (MQ=255)
                                          gCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACCCTATCGacac                      <  1:2055963/69‑1 (MQ=255)
                                                      gccccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:1298579/3‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATc                                    >  1:2601069/1‑43 (MQ=38)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:2739987/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:2637229/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:3253542/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:3069608/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:2971765/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:2774855/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:206325/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgcg           >  1:538785/1‑68 (MQ=255)
                                                      gACCCATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACcgc            >  1:784986/1‑67 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:845311/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:485121/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:833075/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:417460/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:2693437/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:3309335/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:1311128/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:1491254/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:2181869/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:2603216/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:2650567/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:2675791/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:1085036/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:2730244/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:3248105/68‑1 (MQ=255)
                                                         ccATCGTCGGCAACTTCCGCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGccc        <  1:3313382/68‑1 (MQ=255)
                                                            tCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGa    <  1:592988/69‑1 (MQ=255)
                                                            tCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGa    <  1:1970775/69‑1 (MQ=255)
                                                            tCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGa    <  1:1230733/69‑1 (MQ=255)
                                                              gTCAGCAACTTCCTCCCGCCGTGTATGTGTTGGTCGACGCTATCGACACCACGCCCCGCGCCCCCGAAc  >  1:2431871/1‑69 (MQ=255)
                                                              gTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGaa   >  1:535894/1‑68 (MQ=255)
                                                              gTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAc  >  1:2753359/1‑69 (MQ=255)
                                                              gTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAc  >  1:3266201/1‑69 (MQ=255)
                                                              gTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAc  >  1:3431746/1‑69 (MQ=255)
                                                              gTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAc  >  1:980503/1‑69 (MQ=255)
                                                              gTCAGCAACTTCCTACCGCCGTGTATGTGTTGATCGACGc                               >  1:2201133/1‑40 (MQ=255)
                                                                 |                                                                 
CGATGGCGAATTTTGGGCGCAAAGTATGCCCGTCTGGCCAGTGCTGACTTCGCTGACCCATCGTCAGCAACTTCCTCCCGCCGTGTATGTGTTGATCGACGCTATCGACACCACGCACCGCGCCCACGAAC  >  minE/368233‑368363

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: