Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 209,373 T→C 50.7% L111L (TTA→TTG yafC ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE209,3730TC50.7% ‑1.3 / 62.8 71L111L (TTA→TTGyafCpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (27/8);  new base C (6/30);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 3.14e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

TTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACAGTCCACGCG  >  minE/209305‑209434
                                                                    |                                                             
ttAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCttaatt                                                               <  1:1163118/69‑1 (MQ=255)
    taatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                            <  1:1152451/68‑1 (MQ=255)
     aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                            <  1:1013566/67‑1 (MQ=255)
     aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                            <  1:2634719/67‑1 (MQ=255)
     aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                            <  1:377306/67‑1 (MQ=255)
     aatCGTTTCGGAGGAGACTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAAc                                                            <  1:654568/67‑1 (MQ=255)
                    gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                           >  1:1187523/1‑69 (MQ=255)
                    gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                           >  1:1229329/1‑69 (MQ=255)
                    gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                           >  1:26496/1‑69 (MQ=255)
                    gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgcag                                           >  1:3212568/1‑69 (MQ=255)
                    gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgca                                            >  1:1333605/1‑68 (MQ=255)
                    gaCTAGCGATAAAGTGACTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTgca                                            >  1:1040679/1‑68 (MQ=255)
                        aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                                       <  1:936930/69‑1 (MQ=255)
                        aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                                       <  1:311030/69‑1 (MQ=255)
                        aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                                       <  1:410285/69‑1 (MQ=255)
                        aGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGcacc                                       <  1:2220082/69‑1 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:2911670/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:2018700/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:3262833/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:3306021/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:33154/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:963912/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:487837/1‑69 (MQ=255)
                                 gTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGtt                              >  1:1121030/1‑69 (MQ=255)
                                    aCTTCCGGATATCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:155184/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:2794134/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:1018506/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:2486163/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:3015857/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:3321795/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:3433369/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:56830/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:681001/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:731244/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:135504/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:1409816/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:157128/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:2412759/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:215074/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:2064912/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:1993744/69‑1 (MQ=255)
                                    aCTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGATGCg                           <  1:353316/69‑1 (MQ=255)
                                     cTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:1679639/68‑1 (MQ=255)
                                     cTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:986309/68‑1 (MQ=255)
                                     cTTCCGGATAGCGTTCACGAAAAGGCTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:1697598/68‑1 (MQ=255)
                                       tCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCa                        <  1:891330/69‑1 (MQ=255)
                                       tCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCa                        <  1:1314802/69‑1 (MQ=255)
                                           gATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGa                    <  1:3369445/69‑1 (MQ=255)
                                           gATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGa                    <  1:2768837/69‑1 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:1457808/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:2486560/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:1023861/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:881218/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:1305014/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:2361209/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:1444865/1‑69 (MQ=255)
                                                    cACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACa           >  1:1487498/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)
                                                             gCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAg                                >  1:328469/1‑39 (MQ=255)
                                                             gCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACAGTCCAcgcg  >  1:2969295/1‑69 (MQ=255)
                                                             gCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACAGTCCAcgcg  >  1:834933/1‑69 (MQ=255)
                                                             gCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACAGTCCAcgcg  >  1:926675/1‑69 (MQ=255)
                                                             gCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACAGTCCAcgcg  >  1:2115018/1‑69 (MQ=255)
                                                              cTTAATCAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:3408584/43‑1 (MQ=255)
                                                                 aattaaCGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGa                    <  1:406807/47‑1 (MQ=255)
                                                                 aatcaaCGGCATCAGAAAGTGCAGCACCACTGGAGTTGCg                           <  1:2309038/40‑1 (MQ=255)
                                                                    |                                                             
TTAATAATCGTTTCGGAGGAGACTAGCGACAAAGTGACTTCCGGATAGCGTTCACGGAAAGGCTTAATTAACGGCATCAGAAAGTGCAGCACCACTGGAGTTGCGGCATCGATCCGTAACAGTCCACGCG  >  minE/209305‑209434

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: