Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,006,3090TC41.7% 9.7 / 43.4 36G64G (GGT→GGC) ydgHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/21);  new base C (0/15);  total (0/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTT  >  minE/1006259‑1006320
                                                  |           
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCATATAAAGAAGGCGCCGCCTCttt  <  1:872644/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1032429/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:3109869/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:363729/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:380476/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:2748054/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:2637993/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:2589226/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:499446/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:682883/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:2258422/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:218767/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1972165/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1937641/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:831588/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1803469/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1496025/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1155921/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1047260/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTCCCGCCTCttt  <  1:1085300/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:895571/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:423679/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:972527/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:3636501/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:3426472/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:323576/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:3187447/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:2920159/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:2689832/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:2506069/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:2343667/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:1631710/62‑1 (MQ=255)
gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:1189987/62‑1 (MQ=255)
  tATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGCGCCGCCTCttt  <  1:308945/60‑1 (MQ=255)
          aaGCTGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGACTCttt  <  1:1013168/52‑1 (MQ=39)
                      ccGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCttt  <  1:1881120/40‑1 (MQ=255)
                                                  |           
GCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTT  >  minE/1006259‑1006320

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: ‑