Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 128,172 C→T 19.6% T218T (ACC→ACT yadG → predicted transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE128,1720CT19.6% 122.6 / 18.0 81T218T (ACC→ACTyadGpredicted transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (45/20);  new base T (2/14);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 9.22e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GGAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAG  >  minE/128106‑128237
                                                                  |                                                                 
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:3055235/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:1828545/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:1885498/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:414436/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:414752/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:2088082/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:221412/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTc                                                                 >  1:2714142/1‑69 (MQ=255)
ggAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAACACCTc                                                                 >  1:1736099/1‑69 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTTGTGGAAAATACCTCGa                                                               >  1:2077246/1‑67 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGa                                                               >  1:142846/1‑67 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGa                                                               >  1:318255/1‑67 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGa                                                               >  1:1649947/1‑67 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGa                                                               >  1:406159/1‑67 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGa                                                               >  1:2059054/1‑67 (MQ=255)
    gaagCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTAGa                                                               >  1:2909526/1‑67 (MQ=255)
                   tgCCGCAATATCGGCATTATTCAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTgg                                                >  1:125587/1‑67 (MQ=255)
                   tgCCGCAATATCGGCATTATTCAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTgg                                                >  1:2122301/1‑67 (MQ=255)
                     ccGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGaa                                            >  1:1913327/1‑69 (MQ=255)
                     ccGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGa                                             >  1:2251460/1‑68 (MQ=255)
                             tCGTCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:2759122/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:3139058/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:3246839/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:2760059/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:2060638/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:2318210/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:1045330/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:1052498/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:1070861/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:869875/51‑1 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAAt                                    >  1:2844235/1‑69 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAAt                                    >  1:764377/1‑69 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAAt                                    >  1:1364471/1‑69 (MQ=255)
                             tCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAAt                                    >  1:728674/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGgcgc                                                      <  1:1734689/43‑1 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:1807008/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:293752/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:2920734/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:2487213/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:2468273/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:3287747/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:3296019/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:2324854/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:699926/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:3323880/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:1420860/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:1544280/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:489398/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:1620346/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:3990/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAcc                            >  1:3430169/1‑69 (MQ=255)
                                     attCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAAc                             >  1:2807377/1‑68 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:846961/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:101725/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:3401598/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:941407/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:2747614/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:2410885/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1105935/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:2098946/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1261121/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1756806/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1404722/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1709034/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1702309/68‑1 (MQ=255)
                                        cAACACGGTGAGTTGGTGGAAAATACTTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCAGAAACCtt                          <  1:1597458/68‑1 (MQ=255)
                                            acGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTAtt                      >  1:2464895/1‑68 (MQ=255)
                                            acGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATtc                     >  1:906238/1‑69 (MQ=255)
                                            acGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATtc                     >  1:2176725/1‑69 (MQ=255)
                                                          gAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCAcc        >  1:1224440/1‑68 (MQ=255)
                                                              aTACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGaaaa   <  1:2704457/69‑1 (MQ=255)
                                                              aTACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGaaaa   <  1:402195/69‑1 (MQ=255)
                                                              aTACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGaaaa   <  1:3433279/69‑1 (MQ=255)
                                                              aTACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGaaaa   <  1:3389982/69‑1 (MQ=255)
                                                              aTACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTATATTCTCGATCTCGCACCGaaaa   <  1:1127799/69‑1 (MQ=255)
                                                               tACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAg  <  1:1635889/69‑1 (MQ=255)
                                                               tACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAg  >  1:524073/1‑69 (MQ=255)
                                                               tACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAg  <  1:687492/69‑1 (MQ=255)
                                                               tACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAg  <  1:1831615/69‑1 (MQ=255)
                                                               tACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAg  >  1:3296547/1‑69 (MQ=255)
                                                               tACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAg  <  1:1105283/69‑1 (MQ=255)
                                                                  |                                                                 
GGAAGAAGCAGAAATGCTGTGCCGCAATATCGGCATTATTCAACACGGTGAGCTGGTGGAAAATACCTCGATGAAGGCGCTGCTGGCGAAGCTGAAATCGGAAACCTTTATTCTCGATCTCGCACCGAAAAG  >  minE/128106‑128237

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: