Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,370,3250CA50.0% 23.2 / 71.5 50T93N (ACC→AAC) yohKpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (22/3);  new base A (0/25);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 5.18e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGT  >  minE/1370276‑1370381
                                                 |                                                        
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1597524/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:2696585/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:3341226/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:2462466/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:3376342/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:2190204/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:34967/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:2040562/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1960934/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1918111/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1907798/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:2819649/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1682961/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1564173/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:2847036/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1458974/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1339592/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:416433/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1297608/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1246372/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:614771/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:1186843/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:115663/62‑1 (MQ=255)
ttATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:3198429/62‑1 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3088996/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3215698/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:2472920/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3288188/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:350272/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3617047/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3648383/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:66870/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:966747/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:1924435/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:1024324/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:1325434/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:1860650/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3077697/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:2182853/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:2233007/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:1007230/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:2797419/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:2812882/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCa                                         >  1:3049545/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTCATCGGCa                                         >  1:29120/1‑62 (MQ=255)
     tGAGCAGCTACACCAGATCCGCGCGCGCAGGAAATCGATCATCACCATCTGTTTTTTCGGCa                                         >  1:390230/1‑62 (MQ=255)
       agcagcTACACCAGATCCGCGCGCGCTGGAAATCGATCATCAACATCTGTTTTAt                                              <  1:751951/55‑1 (MQ=255)
                                            catcaCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCAttgt  <  1:1501724/62‑1 (MQ=255)
                                            catcaCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCAttgt  <  1:1240347/62‑1 (MQ=255)
                                            catcaCCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCAttgt  <  1:1096210/62‑1 (MQ=255)
                                                 |                                                        
TTATATGAGCAGCTACACCAGATCCGCGCGCGCTGGAAATCGATCATCACCATCTGTTTTATCGGCAGCGTGGTTGCGATGGTGACGGGGACTTCCGTGGCATTGT  >  minE/1370276‑1370381

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: