Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,893,7940TA37.7% 46.4 / 58.2 53Y141F (TAC→TTC) mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base T (33/0);  new base A (20/0);  total (53/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGC  >  minE/1893730‑1893821
                                                                |                           
aCTCAATGCCCCCGCGTAGATCTCCGCCCGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:976607/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:204113/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1097423/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:254369/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:279174/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:363559/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:369019/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:375573/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:394231/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:406599/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:426909/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:470135/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:572399/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:647690/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:668845/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:679209/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:758605/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:788793/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:925832/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:954549/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:972276/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:158002/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1045129/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1146130/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1177940/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1337957/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1343699/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1347403/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1364148/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1392514/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:169253/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1504499/1‑68 (MQ=255)
aCTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTaaa                          >  1:1529895/1‑68 (MQ=255)
                     ctcCGCACGAGACGGCAGCCGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTg    >  1:421811/1‑69 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCCCGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:1549868/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCg         >  1:1098330/1‑64 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTgg   >  1:556798/1‑70 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:3055/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:967958/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:114140/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:266813/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:899463/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:116457/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:767715/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:270693/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:751119/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:14262/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:344268/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:5133/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:1367547/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:1375817/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:1401152/1‑71 (MQ=255)
                     ctcCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGAAAATAGGATACTGGAACTCGCCCTGGc  >  1:376920/1‑71 (MQ=255)
                                                                |                           
ACTCAATGCCCCCGCGTAGATCTCCGCACGAGACGGCAGACGACCACGTTTTGGCGGCATACGGTAAATAGGATACTGGAACTCGCCCTGGC  >  minE/1893730‑1893821

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