Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,053,033 A→G 27.9% G97G (GGT→GGC ydhF ← predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,053,0330AG27.9% 128.9 / 54.6 126G97G (GGT→GGCydhFpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (73/17);  new base G (6/29);  total (80/46)
Fisher's exact test for biased strand distribution p-value = 4.15e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

Only 100 of 126 total aligned reads displayed.

GGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTC  >  minE/1052967‑1053098
                                                                  |                                                                 
ggTCGCGAGATTCATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGAcc                                                                 <  1:896862/69‑1 (MQ=255)
ggTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGAcc                                                                 <  1:663967/69‑1 (MQ=255)
ggTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGAcc                                                                 <  1:695688/69‑1 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:2897227/1‑69 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:1715506/1‑69 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:1711037/1‑69 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:2165251/1‑69 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:23407/1‑69 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:2454732/1‑69 (MQ=255)
                aGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCa                                                 >  1:1153056/1‑69 (MQ=255)
                               ctctTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTcgc                                   >  1:1595544/1‑68 (MQ=255)
                               ctctTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTcgc                                   >  1:452983/1‑68 (MQ=255)
                                  ttAATGATGTGATCGCGGTCAGTGATGTAATGGCCAATGACGTTttct                                                    >  1:1776745/1‑48 (MQ=255)
                                  ttAATGATGTGATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGAt                                >  1:1426009/1‑68 (MQ=255)
                                  ttAATGATGTGATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATa                               >  1:760989/1‑69 (MQ=255)
                                    aaTGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAc                              >  1:387119/1‑68 (MQ=255)
                                    aaTGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAc                              >  1:717263/1‑68 (MQ=255)
                                      tgatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAc                              <  1:449774/66‑1 (MQ=255)
                                      tgatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAc                              <  1:2570627/66‑1 (MQ=255)
                                      tgatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAc                              <  1:2078292/66‑1 (MQ=255)
                                       gatgTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAc                              <  1:1103229/65‑1 (MQ=255)
                                        aagtgATCGCGGTCAGTGATGTAATGCCCAATGACGTTTTCTTCACGCGCGGtcgtcg                                    >  1:1021025/3‑58 (MQ=39)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:931198/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:354230/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:1097768/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:774895/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:2050714/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:728046/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:195949/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:309659/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCAtt                        <  1:2904516/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGGCCAATGACGTTTTCTTAACACGCGGTCGTCGCGATACCGCAtt                        <  1:2378460/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATAGCGCAtt                        <  1:440176/68‑1 (MQ=255)
                                          gtgATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCAtt                        <  1:160077/68‑1 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGa                                 >  1:2331694/1‑56 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:3212323/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:2710832/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:3271501/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:428947/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:1754388/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:23783/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:504207/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:476786/1‑69 (MQ=255)
                                             aTCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACt                    >  1:2248154/1‑69 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:323940/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:2902793/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:1699445/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:2013655/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:3224290/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:2570834/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:2321232/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTg                   >  1:1375709/1‑68 (MQ=255)
                                               cgcgGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTGACTg                   >  1:3348294/1‑68 (MQ=255)
                                                      aGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGa                                 >  1:2440893/1‑47 (MQ=255)
                                                      aGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGAt                                >  1:3422616/1‑48 (MQ=255)
                                                        tGATGTAATGGCCTATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGAAGATTTcc          <  1:636031/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:2652178/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:3035417/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:896631/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:1822734/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:1914299/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:1745948/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:2028035/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:2460304/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:2574854/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:2579952/68‑1 (MQ=255)
                                                        tGATGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:2608942/68‑1 (MQ=255)
                                                          aTGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATc       >  1:511908/1‑69 (MQ=255)
                                                          aTGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATc       >  1:1716253/1‑69 (MQ=255)
                                                          aTGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATc       >  1:341679/1‑69 (MQ=255)
                                                          aTGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATc       >  1:1514041/1‑69 (MQ=255)
                                                           tGTAATGGCCAATGACGTTTTCTTCACACGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:1062668/65‑1 (MQ=255)
                                                             tAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCTCGATACCGCATTTACTGACGATTTcc          <  1:292422/63‑1 (MQ=255)
                                                             tAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:313364/63‑1 (MQ=255)
                                                             tAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:3016962/63‑1 (MQ=255)
                                                             tAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:1640393/63‑1 (MQ=255)
                                                             tAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTcc          <  1:1209053/63‑1 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTcaa         >  1:2140786/1‑60 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCa         >  1:2598755/1‑62 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCa         >  1:3131346/1‑62 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCa         >  1:1465863/1‑62 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCa         >  1:1392881/1‑62 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCa         >  1:1736128/1‑62 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCa         >  1:2995673/1‑62 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGtt   >  1:249215/1‑68 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGtt   >  1:66718/1‑68 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:394883/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:1722424/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:454999/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:219359/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:2427461/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:2843269/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:384814/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:2660021/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:269396/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:276905/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:2792347/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:3004346/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTc  >  1:2839443/1‑69 (MQ=255)
                                                               aTGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGAATTCCATCCGTTc  >  1:538572/1‑69 (MQ=255)
                                                                  |                                                                 
GGTCGCGAGATTAATTAGCGACTGTTCGGCGCTCTTAATGATGTGATCGCGGTCAGTGATGTAATGACCAATGACGTTTTCTTCACGCGCGGTCGTCGCGATACCGCATTTACTGACGATTTCCATCCGTTC  >  minE/1052967‑1053098

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: