Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	513,101	1	.	T	64.2%	215.4 / 129.1	95	coding (636/1284 nt)	ybhC	predicted pectinesterase
Reads supporting (aligned to +/- strand):  ref base . (34/0);  new base T (61/0);  total (95/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.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCGCAG  >  minE/513075‑513144                            |                                            tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2689527/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:311694/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2933872/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2903771/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2846867/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2826085/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:280618/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2766057/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2744911/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3192496/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2633340/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2416578/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2394241/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2331387/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2312436/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2304740/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2254997/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2227853/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3617073/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:976001/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:833507/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:789217/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:730482/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:577524/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:526853/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:449576/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3627217/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1042695/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3603510/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3591169/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3547889/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3483682/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3419371/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3389529/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3226514/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:3208892/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1550011/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:203138/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1969397/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1936838/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1830098/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:168828/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:159835/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2159259/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:15502/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1576178/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1495781/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1453735/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1334694/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1309829/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1279902/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1092675/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:1034173/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2038257/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2066236/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2144305/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2117402/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2111961/1‑59 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAc               >  1:367770/1‑58 (MQ=39) tCGCCCAGCGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGACCag            >  1:2440438/1‑61 (MQ=255) tCGCCCAGAGTGTTTTCGATGGTCAGATTTTTGCAGTTGCAGGCCATTGTTTTGTGAcc              >  1:2149506/1‑59 (MQ=38)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTttgcagttg                                  >  1:1414614/1‑36 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:2178452/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:1199960/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:554084/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:1227533/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:385626/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:3615850/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:3468859/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:2708132/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:1981097/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaaga        >  1:2880695/1‑62 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaag         >  1:957310/1‑61 (MQ=255)   gCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCagaag         >  1:2478778/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCATTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:2110677/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:2154735/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:1133197/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:1146163/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:908178/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:1251/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:920417/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:3410917/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:2595734/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:3646965/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:2992732/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:129731/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:1585312/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:3222620/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGAc       >  1:2085181/1‑61 (MQ=255)     ccAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGACATTGTTTTGTGACCAGAAGAc       >  1:3592699/1‑61 (MQ=255)            gtTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCgcag  >  1:2077409/1‑59 (MQ=255)            gtTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCgcag  >  1:2668868/1‑59 (MQ=255)            gtTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCgcag  >  1:414481/1‑59 (MQ=255)            gtTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCgcag  >  1:3177688/1‑59 (MQ=255)            gtTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCgcag  >  1:2345934/1‑59 (MQ=255)                            |                                            TCGCCCAGCGTGTTTTCGATGGTCAGA‑TTTTGCAGTTGCAGGCCATTGTTTTGTGACCAGAAGACCGCAG  >  minE/513075‑513144

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