Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	389336	389371	36	2 [1]	[1] 2	hemB	porphobilinogen synthase

GCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTCGACGCCATGCTCGCACAGCACAC  >  NC_000913/389207‑389336                                                                                                                                 |  ccccgTACCCGGCCGTCCATCGCGGCCGAAGGGGCGATGAAAGTTGCACCTGCAGCAGCTTCAACCCAGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTCGACGCCATGCTCgcacagcaca   <  1:287715/125‑1 (MQ=255) cccTGTCCCGGGCCGCCCACCG‑GGCGGAAAGGGCGAAGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTCGACGCCATGCTCGCACAGCacac  <  1:285286/128‑1 (MQ=255)                                                                                                                                 |  GCCTGTACCTGGCCGTCCATCGCGGCGGAAGGGGCGATGAAGTCTGCACCTGCAGCAGCTGCAACCACGGCTTGCTTGCCTAAATTTTCCAGAGTCGCGTCGTTGTCGACGCCATGCTCGCACAGCACAC  >  NC_000913/389207‑389336

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