Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,240,9091.A61.7% 195.8 / 144.2 94coding (342/870 nt)yebKpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (23/13);  new base A (0/58);  total (23/71)
Fisher's exact test for biased strand distribution p-value = 4.67e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAA  >  minE/1240868‑1240956
                                           |                                               
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:961243/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:2399082/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:2174750/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:2933493/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:3029869/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:3246232/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:3324995/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:3328949/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:1602971/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:362784/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:949193/61‑1 (MQ=255)
tCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACcgcg                              <  1:1200789/61‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGTTTAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1643030/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTTCTCACTCAgg         <  1:337697/62‑1 (MQ=38)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:919223/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3004419/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2811056/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2799554/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2764102/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:275614/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2678417/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2559405/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2509080/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2489832/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2453801/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2401119/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:248401/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2369286/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2950077/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:492854/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:80907/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3105877/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:780513/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:648955/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3359120/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3402841/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3412687/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:637032/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3414900/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:567031/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:38339/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:478578/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1878695/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1517908/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1646351/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1728154/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1480299/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1775669/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1807941/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1809569/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1834286/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:147301/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1860246/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1376020/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1308048/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1987658/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2061494/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2078035/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1206808/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:113614/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2242401/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1112220/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2331676/62‑1 (MQ=255)
                     cATGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCCCCGTCGACTTGCTCACTCAgg         <  1:36838/62‑1 (MQ=255)
                      aTGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:1250572/61‑1 (MQ=255)
                      aTGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:557590/61‑1 (MQ=255)
                      aTGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:894158/61‑1 (MQ=255)
                      aTGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:3287561/61‑1 (MQ=255)
                      aTGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:2273656/61‑1 (MQ=255)
                      aTGTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGGCGACTTGCTCACTCAgg         <  1:2303493/61‑1 (MQ=255)
                        gTCCGTCATTCACTGGATAAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAgg         <  1:103312/59‑1 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCg                        >  1:448412/1‑40 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:185545/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:2430861/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:1157801/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:2168352/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:827537/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:817197/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:1185509/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:2058286/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:161637/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:1412757/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:46841/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:1728406/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:381168/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:3390790/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:1629831/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:3473670/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaaa  >  1:1537268/1‑62 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaa   >  1:183762/1‑61 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaa   >  1:1896468/1‑61 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaa   >  1:2768766/1‑61 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaaa   >  1:2078478/1‑61 (MQ=255)
                           cGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCaaaa    >  1:389514/1‑60 (MQ=255)
                                           |                                               
TCCGCAATGGCAACGCTTGATCATGTCCGTCATTCACTGGATAAATCTGCCATCAACCGCGCCGTCGACTTGCTCACTCAGGCAAAAAA  >  minE/1240868‑1240956

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: