Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	W3110S.gb	4146667	4146668	2	10 [0]	[0] 16	nirB/tsgA	nitrite reductase, large subunit, NAD(P)H‑binding/predicted transporter

GATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGC  >  W3110S.gb/4146605‑4146666                                                              | gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:1268246/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:1289565/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:1391864/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:197530/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:2048361/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:2483666/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:438415/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:85475/62‑1 (MQ=255) gATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:931390/62‑1 (MQ=255)   tGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGc  <  1:205640/60‑1 (MQ=255)                                                              | GATGCACGTAGCGGGGGAAGAGAGGAAAAGGCCGGACGACGATAGTGCCGTCCGGCATGAGC  >  W3110S.gb/4146605‑4146666

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