Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,065,0490CT41.7% 83.2 / 151.1 120Y430Y (TAC→TATrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base C (70/0);  new base T (0/50);  total (70/50)
Fisher's exact test for biased strand distribution p-value = 5.45e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 120 total aligned reads displayed.

AGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCA  >  minE/2065001‑2065107
                                                |                                                          
ggCGGTTGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2560383/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2519057/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1962716/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3599425/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2213685/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2252727/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2256178/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3572035/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3534422/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:23306/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2342826/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2444982/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2467463/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3522409/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:907328/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2535117/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3473958/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3327492/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2741361/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3312722/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2768813/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2875146/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2875193/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3211259/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:2932611/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:3208731/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1060166/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1317036/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1299511/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1390847/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1168523/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1481656/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1157873/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:152487/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:668823/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1628038/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:590326/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:17835/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:589133/53‑1 (MQ=255)
ggCGGTGGACAAGTTCGAATACCGTCGCGGCTATAAATTCTCGACCTATGCGAc                                                       <  1:1862745/53‑1 (MQ=255)
                                             cTACGCAACCTGGTTGATCCGTCAGGCGATCACCCGCTCTAt                      >  1:3531941/1‑42 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGCTCAGGCGCGCACCa  >  1:3507203/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATc              >  1:1210462/1‑50 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCCCCa  >  1:2939660/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:843431/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3161469/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:78992/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2693464/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3181639/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3177617/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:431438/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:337592/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3418957/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:74250/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:730990/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3507366/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:657627/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:630706/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:618521/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3578057/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:510619/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3594360/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:376632/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:427515/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:170214/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3109852/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2150912/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2044868/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2017622/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1938951/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1855412/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1816720/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1713081/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2181235/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1644734/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1605311/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1579589/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1444669/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1439675/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1020864/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1011266/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2281432/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3099088/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:3021587/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:300832/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2933503/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2928801/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2755944/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2698022/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:1007561/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2635871/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2623396/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2488432/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2473913/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:247227/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2310134/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2292368/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCa  >  1:2267706/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCAGATCAGGCGCGCACCa  >  1:614301/1‑62 (MQ=255)
                                             cTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATAGCGGATCAGGCGCGCACCa  >  1:1076315/1‑62 (MQ=255)
                                                |                                                          
AGCGGTTGATAAATTCGAATACCGCCGTGGTTACAAGTTCTCCACCTACGCAACCTGGTGGATCCGTCAGGCGATCACCCGCTCTATCGCGGATCAGGCGCGCACCA  >  minE/2065001‑2065107

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: