Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 331,123 T→C 56.4% K195K (AAA→AAG 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,1230TC56.4% 24.6 / 66.0 87K195K (AAA→AAGybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base T (13/25);  new base C (5/44);  total (18/69)
Fisher's exact test for biased strand distribution p-value = 7.95e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCCATCAT  >  minE/331069‑331188
                                                      |                                                                 
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:378055/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:1880421/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:187748/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:1632980/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:2638145/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:2983888/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:3013669/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:2881382/68‑1 (MQ=255)
ccACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGa                                                      <  1:1229393/68‑1 (MQ=255)
         cAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATc                                            <  1:1340193/69‑1 (MQ=255)
         cAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATc                                            <  1:2087805/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATGAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1767707/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:2344647/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:204835/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:593805/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:551913/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:2225451/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:414332/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:617335/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:3369465/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:3075091/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:2670410/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:2715984/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:3365585/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1994482/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1321909/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1438982/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1287736/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:156661/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1582551/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1711497/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:1790646/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:721144/69‑1 (MQ=255)
         cAGACGGCGACCTACCAACATCATAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATc                                            <  1:996039/69‑1 (MQ=255)
          aGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATc                                            <  1:3404011/68‑1 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:931784/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:2152242/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:2302986/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:2322369/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:1543855/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:3225150/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:818182/1‑68 (MQ=255)
              ggCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACt                                        >  1:1321650/1‑68 (MQ=255)
                                tAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATCg                                           >  1:1061900/1‑47 (MQ=255)
                                tAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCca                     >  1:2895158/1‑69 (MQ=255)
                                tAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCca                     >  1:700023/1‑69 (MQ=255)
                                tAATGGCGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCca                     >  1:1653105/1‑69 (MQ=255)
                                tAATGGAGATAAATGCGATCACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCca                     >  1:1400822/1‑69 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:824046/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:75132/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:3100918/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:447890/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:2683791/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:2468192/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:2926565/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:2469121/69‑1 (MQ=255)
                                 aaTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:1494948/69‑1 (MQ=255)
                                  aTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:2928831/68‑1 (MQ=255)
                                  aTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:86912/68‑1 (MQ=255)
                                  aTGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCcac                    <  1:1446152/68‑1 (MQ=255)
                                             tGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTc        >  1:485887/1‑69 (MQ=255)
                                             tGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTc        >  1:788194/1‑69 (MQ=255)
                                                 aTCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcat     <  1:3286989/68‑1 (MQ=255)
                                                   cACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  >  1:3254932/1‑69 (MQ=255)
                                                   cACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  >  1:1314196/1‑69 (MQ=255)
                                                   cACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatca   >  1:1218473/1‑68 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2015822/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:651538/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:671473/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:1116554/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2010177/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:182579/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:742875/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2947553/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2142620/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2147329/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:545052/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2364354/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2392137/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:392153/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2514452/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2786443/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:3194646/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:3191902/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:3079970/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:3069133/69‑1 (MQ=255)
                                                   cACCTTGCCGATGGTGATCCCCATATAGACTGAAAGAGTGGCAAAGCCCACATCGCCCTGTTCcatcat  <  1:2910134/69‑1 (MQ=255)
                                                      |                                                                 
CCACGGCACCAGACGGCGACCTACCAGCATCATAATGGCGATAAATGCGATCACTTTGCCGATGGTGATCCCCATATCGACTGCAAGAGTGGCAAAGCCCACATCGCCCTGTTCCATCAT  >  minE/331069‑331188

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