Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,474,270 T→C 28.2% A618A (GCT→GCC nrdA → ribonucleoside diphosphate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,474,2700TC28.2% 91.6 / 40.0 92A618A (GCT→GCCnrdAribonucleoside diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (14/52);  new base C (19/7);  total (33/59)
Fisher's exact test for biased strand distribution p-value = 6.12e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

ACGACTGGGAAGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAACCGCC  >  minE/1474206‑1474329
                                                                |                                                           
aCGACTGGGAAGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGa                                                         >  1:2088796/1‑69 (MQ=255)
     tGGGAAGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGctct                                                           >  1:2716776/1‑62 (MQ=255)
     tGGGAAGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGctct                                                           >  1:1491506/1‑62 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:890270/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1834360/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:2735451/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1760106/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1957585/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:3077543/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1740838/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:335529/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1611049/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1582970/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:636405/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1060106/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1141173/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1366440/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1394872/69‑1 (MQ=255)
          aGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:79213/69‑1 (MQ=255)
             tctGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGAc                                            <  1:2122557/69‑1 (MQ=255)
             tctGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGAc                                            <  1:2812995/69‑1 (MQ=255)
             tctGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGAc                                            <  1:1879203/69‑1 (MQ=255)
             tctGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGAc                                            <  1:1126778/69‑1 (MQ=255)
             tctGCGTGAGTCAATCAAAACGCACGGTCAGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGAc                                            <  1:1373330/69‑1 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2491554/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:251219/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2434097/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:1963191/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2330984/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2640794/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2431377/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2116876/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2386763/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2839530/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:1753649/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:2854378/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:1669939/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:1628732/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:404165/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:462645/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:491430/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:758667/1‑47 (MQ=255)
                      gTCAATCAAAACGCACGGCCTGCGTAACTCCACGCTTTCTGCCCTGa                                                         >  1:1304599/1‑47 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTATCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:2741199/69‑1 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGActtctt                                       >  1:3298467/1‑64 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGActtctt                                       >  1:2413198/1‑64 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:2625544/69‑1 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:1037210/69‑1 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:1981933/69‑1 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:1586136/69‑1 (MQ=255)
                       tcaatcaaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:1444281/69‑1 (MQ=255)
                            caaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGAct                                           >  1:2113669/1‑55 (MQ=255)
                            caaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGAtct                              >  1:2558557/1‑68 (MQ=255)
                            caaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGAtct                              >  1:78160/1‑68 (MQ=255)
                            caaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGAtct                              >  1:1850090/1‑68 (MQ=255)
                            caaAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGAtct                              >  1:2934951/1‑68 (MQ=255)
                            caaAAAGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGAtct                              >  1:1010217/1‑68 (MQ=255)
                             aaaaCGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATctct                            <  1:2389219/69‑1 (MQ=255)
                             aaaaCGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATctct                            <  1:410571/69‑1 (MQ=255)
                             aaaaCGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATctct                            <  1:663368/69‑1 (MQ=255)
                               aaCGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAg                                  <  1:3067344/61‑1 (MQ=255)
                                     cGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCa                     >  1:3154077/1‑68 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:629922/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:545371/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:521443/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:711531/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:1431867/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:777958/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:787585/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:81101/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:960723/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:3133039/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:1978789/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:1934157/69‑1 (MQ=255)
                                      ggTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:1757648/69‑1 (MQ=255)
                                       gTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:2882806/68‑1 (MQ=255)
                                        tCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:2960102/67‑1 (MQ=255)
                                        tCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACt                   <  1:1773463/67‑1 (MQ=255)
                                         cTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCga                                               <  1:1752936/38‑1 (MQ=255)
                                                aCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTg         >  1:2250796/1‑69 (MQ=255)
                                                aCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTg         >  1:1412057/1‑69 (MQ=255)
                                                aCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTg         <  1:2041705/69‑1 (MQ=255)
                                                aCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTg         >  1:1756553/1‑69 (MQ=255)
                                                aCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTg         >  1:2933416/1‑69 (MQ=255)
                                                aCTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTg         <  1:1608228/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:1443568/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:1831953/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:1752953/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:1670267/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:833180/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:2107158/69‑1 (MQ=255)
                                                       gCTTTCTGCCCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAAccgcc  <  1:1681329/69‑1 (MQ=255)
                                                         tttCTGCTCTGAAGCCGTCCGAGACTTCTTCGCAg                                  <  1:2460176/35‑1 (MQ=255)
                                                          ttCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATctct                            <  1:78474/40‑1 (MQ=255)
                                                                |                                                           
ACGACTGGGAAGCTCTGCGTGAGTCAATCAAAACGCACGGTCTGCGTAACTCCACGCTTTCTGCTCTGATGCCGTCCGAGACTTCTTCGCAGATCTCTAACGCCACTAACGGTATTGAACCGCC  >  minE/1474206‑1474329

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: