Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	276,639	1	.	T	61.5%	200.6 / 146.4	96	coding (491/615 nt)	cyoC	cytochrome o ubiquinol oxidase subunit III
Reads supporting (aligned to +/- strand):  ref base . (37/0);  new base T (59/0);  total (96/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.75e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGACGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACCGTG  >  minE/276630‑276694           |                                                        cgACTGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2569447/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCTGATAAGACCAGAAGTGACGTGCAGAc      >  1:2076537/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGAcc                       >  1:137197/1‑45 (MQ=38) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGAc                        >  1:728797/1‑44 (MQ=38) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGaa                   >  1:1829199/1‑49 (MQ=38) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTg                >  1:496313/1‑52 (MQ=39) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:197892/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2059276/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2074051/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1961818/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:234874/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2358912/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2392761/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:240858/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2420938/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2462642/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2470758/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1932742/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2516168/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2566175/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2626433/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:3033980/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:93604/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:815102/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:672931/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:48000/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:434087/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:386075/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:3157806/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:314362/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2590388/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2914943/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2796118/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2757936/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:27315/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2725856/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2700432/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2164129/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:162392/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1335908/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1361579/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1372152/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:12601/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1486418/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1226230/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1522347/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1192991/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1609268/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1705408/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1191898/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1185531/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1153891/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1028188/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1783329/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1780492/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1031884/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1744170/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:1137898/1‑62 (MQ=255) cgACGGGCGATTTTGCACCATCAGCACCGCCAACCAGATAAGACCAGAAGTGACGTGCAGAc      >  1:2386331/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCATCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1220065/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCTCCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1705747/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGCCcgtg  >  1:1156985/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2986525/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:922224/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:898526/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1284112/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:275963/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1280530/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2925427/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:602154/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:3055216/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1129720/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2577307/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:3197150/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2227840/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1920446/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:180740/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1958272/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1777941/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:176283/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2125645/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1026314/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2181290/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:222114/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2639314/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1580593/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1535579/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2498956/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2499670/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:2503538/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1456491/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1440198/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1839621/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1431396/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgtg  >  1:1397054/1‑62 (MQ=255)    cGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcgt   >  1:2224261/1‑61 (MQ=255)           |                                                        CGACGGGCGA‑TTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACCGTG  >  minE/276630‑276694

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