Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,010,544 A→G 13.6% intergenic (‑165/‑113) yghY ← / → yghZ ECK2994:JW5496:b3000; predicted dienlactone hydrolase/aldo‑keto reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,010,5440AG13.6% 152.9 / 11.0 87intergenic (‑165/‑113)yghY/yghZECK2994:JW5496:b3000; predicted dienlactone hydrolase/aldo‑keto reductase
Reads supporting (aligned to +/- strand):  ref base A (55/20);  new base G (9/3);  total (64/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TTGGCGGTTAAACGCGGCATCGTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCGCTA  >  minE/2010481‑2010611
                                                               |                                                                   
ttGGCAGTTAAACGCGGCATCGTTGCACTCCTCAACTTGTGTTTTTTTAAGCAAAGCATAAGCGCGtgt                                                                <  1:1442992/69‑3 (MQ=255)
ttGGCAGTTAAACGCGGCATCGTTGCACTCCTCAACTGGTGTTTTTTTAAGCAAAGCATAAGCGCGtgt                                                                <  1:1365440/69‑3 (MQ=255)
ttGGCAGTTAAACGCGGCATCGTTGCACTCCTCAACTGGTGTTTTTTTAAGCAAAGCATAAGCGCGtgt                                                                <  1:1041277/69‑3 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGtt                                                    >  1:2005061/1‑60 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCAc                                           >  1:1032941/1‑69 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCAc                                           >  1:2624133/1‑69 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCAc                                           >  1:318105/1‑69 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCAc                                           >  1:1191576/1‑69 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCAc                                           >  1:1113252/1‑69 (MQ=255)
                     gTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCAc                                           >  1:1059863/1‑69 (MQ=255)
                        gCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGcc                                                         <  1:2771990/52‑1 (MQ=255)
                        gCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGcc                                                         <  1:1299101/52‑1 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:285717/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:2097054/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:2890484/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:2648818/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:752831/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:2388559/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:2570338/1‑68 (MQ=255)
                             cctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:1739870/1‑68 (MQ=255)
                              ctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:133105/1‑67 (MQ=255)
                              ctcAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGa                                    >  1:736197/1‑67 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:1432153/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:144079/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:140183/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:275701/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:1108860/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:718682/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:3115640/69‑1 (MQ=255)
                                cAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:1789615/69‑1 (MQ=255)
                                 aGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:732844/68‑1 (MQ=255)
                                 aGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:408755/68‑1 (MQ=255)
                                 aGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:1553056/68‑1 (MQ=255)
                                     ggTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTa                           >  1:1788851/1‑69 (MQ=255)
                                     ggTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTa                           >  1:2692152/1‑69 (MQ=255)
                                     ggTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTa                           >  1:784733/1‑69 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:861276/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:573548/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:531183/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:1163814/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:365141/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:1352844/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:2584610/1‑68 (MQ=255)
                                       tgtTTTTTTAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTAc                          >  1:1335876/1‑68 (MQ=255)
                                       tggttttttAAGCAAAGCATAAGCGCGTGTTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTa                           >  1:1465634/4‑67 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCtg                                       >  1:3154258/1‑51 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGa                      <  1:1892434/68‑1 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGa                      <  1:3138775/68‑1 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGa                      <  1:2437926/68‑1 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGAt                     >  1:825251/1‑69 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGAt                     >  1:1475098/1‑69 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGAt                     >  1:1548489/1‑69 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGAt                     >  1:874497/1‑69 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGAt                     >  1:1033002/1‑69 (MQ=255)
                                           tttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGAt                     >  1:3030949/1‑69 (MQ=255)
                                            ttttAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGa                      <  1:1938478/67‑1 (MQ=255)
                                                            aGCACTTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                >  1:2489245/1‑41 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAg                                   >  1:73302/1‑38 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                >  1:2792268/1‑41 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                >  1:850417/1‑41 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                <  1:484686/41‑1 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCa                                >  1:464580/1‑41 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATAcgc                 >  1:1029738/1‑56 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATc      >  1:1119464/1‑67 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1180861/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1579461/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1317193/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1240304/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1221832/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:820186/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:690163/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1150913/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:934229/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:975893/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1678633/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:2988517/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1953320/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:1972424/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:2182378/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:239248/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:3395449/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:2791431/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:2926723/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:3083244/1‑68 (MQ=255)
                                                            aGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCg     >  1:2984735/1‑68 (MQ=255)
                                                              cACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCGCTa  <  1:827271/69‑1 (MQ=255)
                                                              cACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCGCTa  <  1:3319694/69‑1 (MQ=255)
                                                               |                                                                   
TTGGCGGTTAAACGCGGCATCGTTGCACTCCTCAGTTGGTGTTTTTTTAAGCAAAGCATAAGCACGTATTTTTGCCCAGTTTTTCGTCACTCTGTGAGCCAGACTACGGGATACGCGCTGGCGAATCGCTA  >  minE/2010481‑2010611

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: