Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,862,627 T→C 59.6% R163R (CGT→CGC ygcN → predicted oxidoreductase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,862,6270TC59.6% 38.6 / 57.1 84R163R (CGT→CGCygcNpredicted oxidoreductase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (6/28);  new base C (11/39);  total (17/67)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGA  >  minE/1862562‑1862695
                                                                 |                                                                    
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:1130323/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:518745/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:47596/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:411797/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:34689/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:3180566/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:300835/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:2832311/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:2025698/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:1851630/1‑69 (MQ=255)
gACGATATTCTTCGCGCCCGTTATGTGGTGCTGGCAGAAGGTCCCAACAGCGTCCTGGCGGAACGCCAc                                                                   >  1:3306068/1‑69 (MQ=255)
 aCGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCAc                                                                   >  1:2020098/1‑68 (MQ=255)
 aCGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACg                                                                  >  1:1260418/1‑69 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTTCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:2251109/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:3185280/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:3308624/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:422081/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:668484/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:202399/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:3227736/69‑1 (MQ=255)
     tatTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGtt                                                              <  1:890137/69‑1 (MQ=255)
                      aTGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCctgctg                                             >  1:3038520/1‑69 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1938974/69‑1 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1828195/69‑1 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:3212437/69‑1 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:2929580/69‑1 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:2161446/69‑1 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:2859072/69‑1 (MQ=255)
                             gCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1126123/69‑1 (MQ=255)
                                       ggTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGaa                                        >  1:1709442/1‑57 (MQ=255)
                                       ggTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTg                            >  1:2462037/1‑69 (MQ=255)
                                       ggTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTg                            >  1:943649/1‑69 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGTATCa                     <  1:2622979/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:1702242/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:291265/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:2717702/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:1345056/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:2685995/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:1687294/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:1560701/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:3001737/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:2485553/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:1267467/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:490025/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:518251/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:318690/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:1784578/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:78140/69‑1 (MQ=255)
                                              aCAGCGTCCTGGCGGAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCa                     <  1:328103/69‑1 (MQ=255)
                                                   gTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:392250/47‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:1227467/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:3329199/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:123904/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:1611257/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:648411/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:667835/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:690400/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:787222/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:874700/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:971396/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:300028/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:1815512/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:1995556/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2102695/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2129870/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2142763/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2178821/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2240795/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2515880/56‑1 (MQ=255)
                                                           ggAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:286277/56‑1 (MQ=255)
                                                            gAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2432717/56‑1 (MQ=255)
                                                            gAACGCCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCATTGGGGATCaa                    <  1:2935360/56‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGTCTCGTCCTGCTGTCGAAGc                                      <  1:435709/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1584100/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:175124/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1801561/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:3307882/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:3311925/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1560221/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:1463777/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:554842/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:2548939/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:2370858/35‑1 (MQ=255)
                                                               cGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGc                                      <  1:2429395/35‑1 (MQ=255)
                                                                 tCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGa  <  1:3371881/69‑1 (MQ=255)
                                                                 tCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGa  <  1:796410/69‑1 (MQ=255)
                                                                 |                                                                    
GACGATATTCTCCGCGCCCGTTATGTGGTGCTGGCAGAAGGTGCCAACAGCGTCCTGGCTGAACGTCACGGGTTAGTGACTCGTCCTGCTGGCGAAGCGATGGCGTTGGGGATCAAAGAAGTGCTGTCGCTGGA  >  minE/1862562‑1862695

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: