Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	USA300TCH1516_ALE	2082828	2082849	22	2 [1]	[0] 2	USA300HOU_RS10520	65 kDa membrane protein

CATTTGCCATTGCTGATTTAAGTTTTATACTTAAATCATTAATTGTAATTTTGCTTGCATTATAGAACACAAATTCATTTGAAACTATGTTTTCAAATTTATTTATTGTTACTGGGAATTTCACATTATTAGATTTATTTATTTTTTT  >  USA300TCH1516_ALE/2082850‑2082997 |                                                                                                                                                    cATTTGCCATTGCTGATTTAAGTTTTATACTTAAATCATTAATTGTAATTTTGCTTGCATTATAGAACACAAATTCATTTGAAACTATGTTTTCAAATTTATTTATTGTTACTGGGAATTTCACATTATTAGATTTATTTAttttttt  <  1:236088/148‑1 (MQ=255) cATTTGCCATTGCTGATTTAAGTTTTATACTTAAATCATTAATTGTAATTTTGCTTGCATTATAGAACACAAATTCATTTGAAACTATGTTTTCAAATTTATTTATTGTTACTGGGAATTTCACATTATTAGATTTATTTAttttttt  <  1:236124/148‑1 (MQ=255) |                                                                                                                                                    CATTTGCCATTGCTGATTTAAGTTTTATACTTAAATCATTAATTGTAATTTTGCTTGCATTATAGAACACAAATTCATTTGAAACTATGTTTTCAAATTTATTTATTGTTACTGGGAATTTCACATTATTAGATTTATTTATTTTTTT  >  USA300TCH1516_ALE/2082850‑2082997

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