Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	minE	687344	687358	15	49 [0]	[0] 37	yccK	predicted sulfite reductase subunit

CACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAT  >  minE/687295‑687343                                                 | cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTGTCCGTTTCTAt  <  1:2134090/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3175574/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2305051/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2428199/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2641465/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2730407/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2815452/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2916560/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2955159/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3026282/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3048535/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3073456/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1092756/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3355396/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3445196/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3450964/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3593286/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3631163/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:383997/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:468277/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:595900/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:703275/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:777258/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:847418/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1866223/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1269190/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1320141/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1597361/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1616806/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1621223/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1708320/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1721938/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1785544/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1242864/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1859319/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:189233/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1943971/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2006803/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1077011/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2060845/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2137618/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2164805/49‑1 (MQ=255) cACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:218904/49‑1 (MQ=255)  aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:3191313/48‑1 (MQ=255)  aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1224388/48‑1 (MQ=255)  aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2246509/48‑1 (MQ=255)  aCTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:2818543/48‑1 (MQ=255)    tCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:1927092/46‑1 (MQ=255)            ctgctTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAt  <  1:866134/38‑1 (MQ=255)                                                 | CACTCCACTGGCTGCTTTCTTTGAGATAGCCTTCGGTATCCGTTTCTAT  >  minE/687295‑687343

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