Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,674,145 A→G 25.5% S131S (TCT→TCC gntT ← gluconate transporter, high‑affinity GNT I system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,674,1450AG25.5% 99.6 / 31.1 86S131S (TCT→TCCgntTgluconate transporter, high‑affinity GNT I system
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base G (13/9);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AAGCCGTGGGTCACAGACAGTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGCAAA  >  minE/2674080‑2674210
                                                                 |                                                                 
aaGCCGTGGGTCACAGACAGTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGaa                                                                >  1:119893/1‑69 (MQ=255)
            acagacagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTa                                                    <  1:3057828/69‑1 (MQ=255)
            acagacagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTa                                                    <  1:2747646/69‑1 (MQ=255)
            acagacagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTa                                                    <  1:3067393/69‑1 (MQ=255)
            acagacagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTa                                                    <  1:1580859/69‑1 (MQ=255)
            acagacagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTa                                                    <  1:2609266/69‑1 (MQ=255)
            acagacagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTa                                                    <  1:1253674/69‑1 (MQ=255)
                 cagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAacac                                               <  1:611220/69‑1 (MQ=255)
                 cagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAacac                                               <  1:460079/69‑1 (MQ=255)
                 cagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAacac                                               <  1:1874922/69‑1 (MQ=255)
                 cagTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAacac                                               <  1:2637905/69‑1 (MQ=255)
                 cagAGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAacac                                               <  1:2282276/69‑1 (MQ=255)
                  agTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGATTATTCGCAGaa                                                                >  1:3246604/1‑51 (MQ=255)
                  agTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACAcc                                              >  1:1645205/1‑69 (MQ=255)
                  agTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACAcc                                              >  1:2282726/1‑69 (MQ=255)
                     gCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAg                                            >  1:959646/1‑68 (MQ=255)
                     gCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:1136927/1‑69 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCcgcg                                                          >  1:2831612/1‑53 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:2058612/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:2142432/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:2597656/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:2326988/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:2768318/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:3065220/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:181266/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:145989/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:3251202/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:3302393/1‑68 (MQ=255)
                      cAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGc                                           >  1:1104042/1‑68 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:2958237/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:2962787/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:1905334/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:1580200/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:3146488/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:1169075/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:2283112/68‑1 (MQ=255)
                        gccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:982687/68‑1 (MQ=255)
                         ccgccATTGGTACACCAACATACAGCAGTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCgg                                         <  1:2323989/67‑1 (MQ=255)
                           gccATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGgcagc                                     <  1:1078629/69‑1 (MQ=255)
                           gccATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGgcagc                                     <  1:318008/69‑1 (MQ=255)
                            ccATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGgcagc                                     <  1:1703837/68‑1 (MQ=255)
                                 ggTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAac                                                 >  1:1764698/1‑51 (MQ=255)
                                 ggTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAac                                                 >  1:469924/1‑51 (MQ=255)
                                 ggTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAac                                                 >  1:3404458/1‑51 (MQ=255)
                                 ggTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAac                                                 >  1:126009/1‑51 (MQ=255)
                                 ggTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGAGATGGTAAac                                                 >  1:1512578/1‑51 (MQ=255)
                                        cAACATAGAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGc                        <  1:1283357/69‑1 (MQ=255)
                                        cAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGc                        <  1:47417/69‑1 (MQ=255)
                                        cAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGc                        <  1:77555/69‑1 (MQ=255)
                                        cAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGc                        <  1:2813999/69‑1 (MQ=255)
                                        cAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGc                        <  1:1372846/69‑1 (MQ=255)
                                        cAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGc                        <  1:1063633/69‑1 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGcacga                           >  1:2456745/1‑55 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:1392446/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:3248395/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:1440676/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:3088840/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:583923/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:2364754/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:2928520/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:2883186/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:934138/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:2070587/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACtt                  >  1:2816770/1‑66 (MQ=255)
                                                 gcagTGGAATATTCGCGGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACGAAACCCACt                   >  1:1868655/1‑65 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:1657954/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:480525/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:2146106/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:3331743/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:939517/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:3268097/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:754393/69‑1 (MQ=255)
                                                   agTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             <  1:892942/69‑1 (MQ=255)
                                                    gTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCagag             >  1:669231/1‑68 (MQ=255)
                                                    gTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             >  1:25965/1‑68 (MQ=255)
                                                    gTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAg             >  1:300876/1‑68 (MQ=255)
                                                       gAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACa         >  1:2427420/1‑69 (MQ=255)
                                                       gAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACa         >  1:1835073/1‑69 (MQ=255)
                                                       gAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACa         >  1:2010490/1‑69 (MQ=255)
                                                       gAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACa         >  1:2241618/1‑69 (MQ=255)
                                                            ttCGCAGAAGCCGCGATGGTAAACACCAGCGgcagc                                     >  1:1450100/1‑36 (MQ=255)
                                                            ttCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGc     >  1:1269474/1‑68 (MQ=255)
                                                            ttCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGc     >  1:3130643/1‑68 (MQ=255)
                                                            ttCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGc     >  1:1786587/1‑68 (MQ=255)
                                                              cGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGCaaa  <  1:1419577/69‑1 (MQ=255)
                                                              cGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGCaaa  <  1:272200/69‑1 (MQ=255)
                                                                 |                                                                 
AAGCCGTGGGTCACAGACAGTGCAGCCGCCATTGGTACACCAACATACAGCAGTGGAATATTCGCAGAAGCCGCGATGGTAAACACCAGCGGCAGCATCAGCACAAAGCCCACTTCATAGAACAGGGCAAA  >  minE/2674080‑2674210

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: