Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4295429	4296042	614	16 [0]	[0] 42	nrfE	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

TGGCGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTAA  >  W3110S.gb/4296043‑4296103 |                                                             tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:683635/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1078679/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:3195/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:365851/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:379181/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:469497/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:539849/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:586630/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:602018/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:60739/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:649677/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:267448/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:683882/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:69094/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:698428/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:788534/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:788762/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:848250/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:886470/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:972603/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:977922/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1397981/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1194670/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1201737/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:120338/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1262408/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1333167/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1339190/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1385714/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1387526/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1390336/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1390914/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:221467/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1404584/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1426115/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1434247/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1487334/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1519202/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1597169/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:1643887/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTaa  <  1:174864/61‑1 (MQ=255) tggcGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGGTATCGCGGGAAATGTTaa  <  1:528978/61‑1 (MQ=255) |                                                             TGGCGAGCGCGGGACGGTGGCCCGGCGGTGCATTTTTCGGGGTTATCGCGGGAAATGTTAA  >  W3110S.gb/4296043‑4296103

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