Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	AE015451	2388272	2388337	66	4 [2]	[3] 5	pyrD/rmf	Dihydroorotate dehydrogenase (quinone)/ribosome modulation factor

TCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGC  >  AE015451/2388199‑2388272                                                                         |  tCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccg   >  1:405449/1‑73 (MQ=255) tCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccg   <  2:405447/73‑1 (MQ=255)           aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccgc  <  1:93047/64‑1 (MQ=255)           aGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAgcccgc  >  2:93046/1‑64 (MQ=255)                                                                         |  TCATCTACAAAGGGCCGGCGCTGATTCGCGAAGCGGTGGATGCTATCGCGGCAATGCCGCGGTGACCAGCCCGC  >  AE015451/2388199‑2388272

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