Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,885,9611.T71.2% 169.7 / 79.7 73coding (1423/2757 nt)barAhybrid sensory histidine kinase, in two‑component regulatory system with UvrY
Reads supporting (aligned to +/- strand):  ref base . (0/21);  new base T (52/0);  total (52/21)
Fisher's exact test for biased strand distribution p-value = 9.22e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTCCCGCCGTCATGGTGGCACCGGGCT  >  minE/1885905‑1885989
                                                         |                            
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:683289/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:2303325/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:392758/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:2423004/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:584571/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:2024273/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:2019899/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:632920/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:680498/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:1728688/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:1691632/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:1677907/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:158570/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:3072109/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:3013049/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:2809146/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:1545971/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:1574968/61‑1 (MQ=255)
ttCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:700403/61‑1 (MQ=255)
 tCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:268040/60‑1 (MQ=255)
       aCGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTc                          <  1:878266/54‑1 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3184801/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3203371/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3253051/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:585726/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3147853/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3144358/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3260989/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3186176/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:370605/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:431848/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:432082/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:519377/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2731758/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:643334/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:681458/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:705285/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:719237/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:759966/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:762389/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:776973/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:867406/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:957947/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1892950/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1102060/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1156188/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1184251/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1526390/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1571199/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1575652/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1644419/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1716555/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1739042/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1799101/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1816614/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:3088438/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1971379/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2435447/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2962331/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2872265/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:1080815/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2700596/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2584576/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2027262/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2392108/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2269738/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2261192/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2219687/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCgggct  >  1:2056802/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGGTGGCACCggg    >  1:409210/1‑60 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTTCCCGCCGTCATGATGGCACCgggct  >  1:2433706/1‑62 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTAttt                           >  1:678830/1‑36 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGATGCCAGTAttt                           >  1:780549/1‑36 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGCTGAAGCCAGTATTTTCCCGCCGTCATCGTGGCACCg      >  1:1357599/1‑58 (MQ=255)
                        tATTCCAGGCCTTCCGACAGGATGATGCCAGTATTTTCCCGCCGTc                  >  1:1210167/1‑46 (MQ=255)
                                                         |                            
TTCCTGAACGCGATCAATCGCGCTTATTCCAGGCCTTCCGACAGGCTGATGCCAGTATTTCCCGCCGTCATGGTGGCACCGGGCT  >  minE/1885905‑1885989

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: