Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 331,345 A→G 24.3% G121G (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
*minE331,3450AG24.3% 93.5 / 26.2 78G121G (GGT→GGCybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (21/38);  new base G (14/5);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 7.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCAGAG  >  minE/331282‑331390
                                                               |                                             
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAGCCg                                            <  1:168475/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAGCCg                                            <  1:2786580/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAGCCg                                            <  1:2739748/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAGCCg                                            <  1:1775281/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAGCCg                                            <  1:415207/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCg                                            <  1:1444843/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCg                                            <  1:572240/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCg                                            <  1:475439/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCg                                            <  1:1224732/67‑1 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:2354273/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:2124132/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:1122334/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:2698682/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:2372465/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:618343/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGaa                                          >  1:704360/1‑69 (MQ=255)
aattaattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGa                                           >  1:32008/1‑68 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGTAAAGACATAAACCGAacac                                       <  1:779474/69‑1 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac                                       <  1:1680476/69‑1 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac                                       <  1:1837038/69‑1 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac                                       <  1:1902062/69‑1 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac                                       <  1:2769932/69‑1 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac                                       <  1:2910533/69‑1 (MQ=255)
   taattGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAacac                                       <  1:711175/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:3381562/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:2464096/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:3046752/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:3299217/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:3327394/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:958634/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:418646/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:577782/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:1023841/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:1239766/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:1697909/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:2177363/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:2293069/69‑1 (MQ=255)
            ttcttcAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:2026523/69‑1 (MQ=255)
             tcttcAAGTGCGCGCAGTAACACCACGGTACTGGCTGTGGAAAGACATAAACCGAACACGATACCGGt                              <  1:2818637/68‑1 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATAcc                                 >  1:183475/1‑47 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:3330334/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           <  1:1109449/69‑1 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:1406166/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:1635324/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:1733624/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           <  1:2816184/69‑1 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           <  1:2060180/69‑1 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:2669705/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:2060474/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:2183222/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccagc           >  1:267383/1‑69 (MQ=255)
                               aaCACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGAccagcccag            >  1:3276963/1‑68 (MQ=255)
                                acacCACGGTACTGGCGGTGGATAGACATAAGCCGAac                                         >  1:3004519/1‑38 (MQ=38)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGaa                                          >  1:43924/1‑37 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:2379424/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:2057620/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:1327276/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:2517660/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:54902/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:1655486/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:443647/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:360536/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:3248923/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:3322975/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:3345823/1‑38 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCa          >  1:2255450/1‑69 (MQ=255)
                                acacCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCa          >  1:1138221/1‑69 (MQ=255)
                                acacCAAGGTACTGGCGGTGGAAAGACATAAGCCGAac                                         >  1:3358934/1‑38 (MQ=37)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:347963/69‑1 (MQ=255)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:1617793/69‑1 (MQ=255)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:1709435/69‑1 (MQ=255)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:777025/69‑1 (MQ=255)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:956472/69‑1 (MQ=255)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:957443/69‑1 (MQ=255)
                                        gTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:3113407/69‑1 (MQ=255)
                                          aCTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:1818482/67‑1 (MQ=255)
                                             ggCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAAc                       >  1:1992032/1‑43 (MQ=255)
                                                     aaaGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCagag  <  1:93973/56‑1 (MQ=255)
                                                               |                                             
AATTAATTGCCGTTCTTCAAGTGCGCGCAGTAACACCACGGTACTGGCGGTGGAAAGACATAAACCGAACACGATACCGGTCATTAACGACCAGCCCAGCACGGCAGAG  >  minE/331282‑331390

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: