Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 330,220 A→G 20.8% G496G (GGT→GGC ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE330,2200AG20.8% 116.6 / 23.0 86G496G (GGT→GGCybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (21/47);  new base G (14/4);  total (35/51)
Fisher's exact test for biased strand distribution p-value = 7.85e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

CATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGC  >  minE/330156‑330277
                                                                |                                                         
cATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCg                                                       >  1:1196223/1‑69 (MQ=255)
cATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCg                                                       >  1:503313/1‑69 (MQ=255)
cATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCg                                                       >  1:1361152/1‑69 (MQ=255)
 aTAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGc                                                      <  1:2365651/69‑1 (MQ=255)
 aTAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGc                                                      <  1:2805311/69‑1 (MQ=255)
 aTAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGc                                                      <  1:1602325/69‑1 (MQ=255)
 aTAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGc                                                      <  1:2053425/69‑1 (MQ=255)
 aTAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGc                                                      <  1:648538/69‑1 (MQ=255)
    aTGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:1255383/69‑1 (MQ=255)
    aTGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:1508399/69‑1 (MQ=255)
    aTGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:488657/69‑1 (MQ=255)
    aTGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:965247/69‑1 (MQ=255)
    aTGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:1040204/69‑1 (MQ=255)
     tGGGCGCGGTCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:1586949/68‑1 (MQ=255)
     tGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:790020/68‑1 (MQ=255)
                aataatCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTc                                                   <  1:2217930/57‑1 (MQ=255)
                   aatCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGaa                                     >  1:1148485/1‑68 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             <  1:1747147/68‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             >  1:2143186/1‑68 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             <  1:1489540/68‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             >  1:2521851/1‑68 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             <  1:2983711/68‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             <  1:3038917/68‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAg                              <  1:2373875/67‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAg                              <  1:311248/67‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAg                              <  1:690625/67‑1 (MQ=255)
                           tatCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAg                              <  1:2376900/67‑1 (MQ=255)
                            atCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAgg                             <  1:1607717/67‑1 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCATAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:2352162/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCGGGATc                          >  1:3112691/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:869849/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:3093926/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:1410553/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:3098009/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:2723385/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:36298/1‑69 (MQ=255)
                             tCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATc                          >  1:1154569/1‑69 (MQ=255)
                               cGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCAt                                                 >  1:1039703/1‑44 (MQ=255)
                               cGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCAt                                                 >  1:1567620/1‑44 (MQ=255)
                                  aTTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCa                     <  1:457031/69‑1 (MQ=255)
                                  aTTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCa                     <  1:1637264/69‑1 (MQ=255)
                                  aTTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCa                     <  1:226254/69‑1 (MQ=255)
                                  aTTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCa                     <  1:233599/69‑1 (MQ=255)
                                  aTTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCa                     <  1:2013303/69‑1 (MQ=255)
                                   tttttCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGTAATCGTCAGGATCAGCCa                     <  1:385182/68‑1 (MQ=255)
                                          cgcgGGCAGATGCCACAATCTCGCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:633339/68‑1 (MQ=255)
                                          cgcgGGCAGATGCCACAATCTCGCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:2054182/68‑1 (MQ=255)
                                          cgcgGGCAGATGCCACAATCTCGCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:3115821/68‑1 (MQ=255)
                                          cgcgGGCAGATGCCACAATCTCGCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:448501/68‑1 (MQ=255)
                                             gggCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACAtt          <  1:2769078/69‑1 (MQ=255)
                                              ggCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACAtt          <  1:1901057/68‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAAGCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:3206430/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:1615198/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:623779/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:2634192/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:1229154/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:2425146/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:1178547/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:2906927/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:2936613/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:2969880/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:1441647/69‑1 (MQ=255)
                                                cAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:2288959/69‑1 (MQ=255)
                                                 aGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:1768507/68‑1 (MQ=255)
                                                 aGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCa       <  1:516111/68‑1 (MQ=255)
                                                    tGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATg   <  1:3153765/69‑1 (MQ=255)
                                                    tGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATg   <  1:2839203/69‑1 (MQ=255)
                                                    tGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATg   <  1:135739/69‑1 (MQ=255)
                                                     gCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATgc  >  1:684127/1‑69 (MQ=255)
                                                     gCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATgc  >  1:1900391/1‑69 (MQ=255)
                                                     gCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATgc  >  1:3430835/1‑69 (MQ=255)
                                                     gCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATgc  >  1:798090/1‑69 (MQ=255)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAg                        >  1:2074364/1‑45 (MQ=38)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCa                     >  1:1396967/1‑48 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCa                     >  1:2561553/1‑48 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1830547/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1174910/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1911938/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1979345/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:558982/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:864090/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1468360/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:253029/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1092363/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:2747900/1‑49 (MQ=39)
                                                       cacaATCTCGCCCGCTTCATAGCCGTTGGGAATCGTCAGGATCAGCCAt                    >  1:1468000/1‑49 (MQ=39)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:1649570/46‑1 (MQ=255)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:997614/46‑1 (MQ=255)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:1504572/46‑1 (MQ=255)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:2510276/46‑1 (MQ=255)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:3320825/46‑1 (MQ=255)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:267945/46‑1 (MQ=255)
                                                                aCCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca              <  1:3182324/46‑1 (MQ=255)
                                                                |                                                         
CATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCACATTCCAGATGC  >  minE/330156‑330277

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