Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ USA300TCH1516_ALE 1573790 1573830 41 2 [1] [1] 2 USA300HOU_RS07705 Glycyl‑glycine endopeptidase ALE‑1

TTGTCGATTTAACAGATTGCCAAGCTTTTCCAAACCATTTACCAGTACTATTAGCTATAGCCTTAGTGTGGTATCCTACAGAACTTTTGGCTGAGCTCCAACCTGAACTTAATTTGCTTGGAATCCCTTTGATTCCGCTCCACATTTTTTTCATTTCGCCGCCAAAATG  >  USA300TCH1516_ALE/1573649‑1573817
                                                                                                                                            |                            
ttgtCGATTTAACAGATTGCCAAGCTTTTCCAAACCATTTACCAGTACTATTAGCTATAGCCTTAGTGTGGTATCCTACAGAACTTTTGGCTGAGCTCCAACCTGAACTTAATTTGCTTGGAATCCCTTTGATTCCGCTcc                              <  1:27990/141‑1 (MQ=255)
                            tCCAAACCATTTACCAGTACTATTAGCTATAGCCTTAGTGTGGTATCCTACAGAACTTTTGGCTGAGCTCCAACCTGAACTTAATTTGCTTGGAATCCCTTTGATTCCGCTCCACATTTTTTTCATTTCGCCGCCAAAATg  >  1:25511/1‑141 (MQ=255)
                                                                                                                                            |                            
TTGTCGATTTAACAGATTGCCAAGCTTTTCCAAACCATTTACCAGTACTATTAGCTATAGCCTTAGTGTGGTATCCTACAGAACTTTTGGCTGAGCTCCAACCTGAACTTAATTTGCTTGGAATCCCTTTGATTCCGCTCCACATTTTTTTCATTTCGCCGCCAAAATG  >  USA300TCH1516_ALE/1573649‑1573817

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