Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 209,361 G→A 39.3% F115F (TTC→TTT yafC ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE209,3610GA39.3% 45.3 / 58.9 89F115F (TTC→TTTyafCpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base A (6/29);  total (32/57)
Fisher's exact test for biased strand distribution p-value = 3.34e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CAAATTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACA  >  minE/209301‑209425
                                                            |                                                                
cAAATTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1563531/69‑1 (MQ=255)
   aTTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCttaa                                                        <  1:1702679/68‑1 (MQ=255)
   aTTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCttaa                                                        <  1:940938/68‑1 (MQ=255)
    ttAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCttaatt                                                      <  1:1163118/69‑1 (MQ=255)
        taatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                   <  1:1152451/68‑1 (MQ=255)
         aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                   <  1:654568/67‑1 (MQ=255)
         aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                   <  1:2634719/67‑1 (MQ=255)
         aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                   <  1:1013566/67‑1 (MQ=255)
         aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                   <  1:377306/67‑1 (MQ=255)
                  ggaggaGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCttaa                                                        <  1:343179/53‑1 (MQ=255)
                        gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                  >  1:26496/1‑69 (MQ=255)
                        gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                  >  1:1187523/1‑69 (MQ=255)
                        gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                  >  1:1229329/1‑69 (MQ=255)
                        gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                  >  1:3212568/1‑69 (MQ=255)
                        gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgca                                   >  1:1040679/1‑68 (MQ=255)
                        gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgca                                   >  1:1333605/1‑68 (MQ=255)
                            aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                              <  1:311030/69‑1 (MQ=255)
                            aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                              <  1:2220082/69‑1 (MQ=255)
                            aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                              <  1:936930/69‑1 (MQ=255)
                            aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                              <  1:410285/69‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:2450218/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          >  1:537862/1‑35 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:2856235/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:533635/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          >  1:1953591/1‑35 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:263149/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:2992418/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1750129/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          >  1:2749515/1‑35 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:2785065/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1574295/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:3260558/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1137079/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1240736/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1302736/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1307861/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:1023099/35‑1 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          >  1:792743/1‑35 (MQ=255)
                                  aaaGTGACTTCCGGATAGCGTTCACGGAAAGGCtt                                                          <  1:3277536/35‑1 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:487837/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:1121030/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:3306021/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:963912/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:3262833/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:2018700/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:2911670/1‑69 (MQ=255)
                                     gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                     >  1:33154/1‑69 (MQ=255)
                                        aCTTCCGGATATCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:155184/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:2794134/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:3321795/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:3433369/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:3015857/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:2412759/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1018506/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:135504/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1409816/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:731244/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:681001/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:157128/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:56830/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1749720/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1803311/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1863007/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1993744/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:2064912/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:215074/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:2486163/69‑1 (MQ=255)
                                        aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGATGCg                  <  1:353316/69‑1 (MQ=255)
                                         cTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:986309/68‑1 (MQ=255)
                                         cTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1697598/68‑1 (MQ=255)
                                         cTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                  <  1:1679639/68‑1 (MQ=255)
                                           tCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCa               <  1:1314802/69‑1 (MQ=255)
                                           tCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCa               <  1:891330/69‑1 (MQ=255)
                                               gATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGa           <  1:3369445/69‑1 (MQ=255)
                                               gATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGa           <  1:2768837/69‑1 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:1023861/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:1305014/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:881218/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:1444865/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:1457808/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:1487498/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:2361209/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:2486560/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:2690593/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:2786711/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:2804840/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:2818146/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa  >  1:3260718/1‑69 (MQ=255)
                                                        cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATACGTAACa  >  1:10766/1‑69 (MQ=255)
                                                            |                                                                
CAAATTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACA  >  minE/209301‑209425

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: