Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,210,179 A→C 29.6% T454T (ACT→ACG yhdP ← conserved membrane protein, predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,210,1790AC29.6% 104.3 / 66.0 98T454T (ACT→ACGyhdPconserved membrane protein, predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (17/52);  new base C (17/12);  total (34/64)
Fisher's exact test for biased strand distribution p-value = 2.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TTGGCTTTAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTG  >  minE/2210114‑2210240
                                                                 |                                                             
ttggctttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGtt                                                            >  1:1146457/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:286350/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:2185524/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:221106/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:2225753/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:1534108/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:2311982/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:2602716/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:948622/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:2673968/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:2921376/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:128381/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:3242212/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:3328134/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:419175/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:466640/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTgg                                                      >  1:547395/1‑69 (MQ=255)
      ttAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTCATTGTTCAGCCAGCTTATCGTTGCCTg                                                       >  1:789975/1‑68 (MQ=255)
           gTCAATATTACGCCCATCCAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGt                                                  >  1:1739727/1‑68 (MQ=255)
           gTCAATATTACGCCCATCCAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTc                                                 >  1:2810146/1‑69 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:2318712/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:2823238/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:1278967/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:483914/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:3019891/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:3023963/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:3196593/69‑1 (MQ=255)
                            ccAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                <  1:1862030/69‑1 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCGGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:2371171/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:51627/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:2063979/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:2983256/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:3223960/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:955685/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcgc                               >  1:756263/1‑69 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                >  1:2524827/1‑68 (MQ=255)
                             cAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGgcg                                >  1:1474400/1‑68 (MQ=255)
                                 tgtaaaCCTTTGTTATTGTTCATCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:1256439/66‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:233546/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:518240/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:3083236/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:3077790/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2001769/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:698585/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:1585351/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:1384828/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2890540/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2868600/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:3242234/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:583971/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2697496/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2319218/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:1234955/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:3331700/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2643743/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2498209/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:589373/69‑1 (MQ=255)
                                 tGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:2580609/69‑1 (MQ=255)
                                                ttGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTGTa                                      >  1:1446326/1‑43 (MQ=255)
                                                ttGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTa                                      >  1:2141374/1‑43 (MQ=255)
                                                ttGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTa                                      >  1:1964678/1‑43 (MQ=255)
                                                ttGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTa                                      >  1:324693/1‑43 (MQ=255)
                                                ttGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTa                                      >  1:3149643/1‑43 (MQ=255)
                                                ttGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACgg                           <  1:1581525/54‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:1010402/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:427261/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:463544/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:310325/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:1918642/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:2109142/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:2260013/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:802641/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:2306608/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:1473693/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:1391871/69‑1 (MQ=255)
                                                 tGTTCAGCCAGCTTATCGTTGCCTGGCCGTCGGCGATTTCCAGTGGCGCACGGAATACCGTTTCGTAAg           <  1:2905816/69‑1 (MQ=255)
                                                  gTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGa                                            <  1:2438228/35‑1 (MQ=255)
                                                  gTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGa                                            <  1:2301125/35‑1 (MQ=255)
                                                  gTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGa                                            <  1:3060812/35‑1 (MQ=255)
                                                  gTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGa                                            <  1:1951873/35‑1 (MQ=255)
                                                  gTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGa                                            <  1:3131600/35‑1 (MQ=255)
                                                     cagccagcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCAt       <  1:3309435/69‑1 (MQ=255)
                                                     cagccagcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCAt       <  1:2448407/69‑1 (MQ=255)
                                                     cagccagcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCAt       <  1:92705/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGTCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:3246141/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:761628/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:990863/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:3307531/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:3007849/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:2662141/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:2658128/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:2535579/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:2305104/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:2156489/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:1470698/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:1187529/69‑1 (MQ=255)
                                                          agcTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTg  <  1:1047671/69‑1 (MQ=255)
                                                             ttATAGTTGCCTGGCCGTAGGCGATTTCTAGTGGCGCACgg                           <  1:455171/41‑1 (MQ=255)
                                                                 |                                                             
TTGGCTTTAACGTCAATATTACGCCCATCCAGCTGGAAACCTTTGTTATTGTTCAGCCAGCTTATAGTTGCCTGGCCGTCGGCGATTTCTAGTGGCGCACGGAATACCGTTTCGTAAGGCATCTTTG  >  minE/2210114‑2210240

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: