Handle missing chromosomes
Closed this issue · 0 comments
kautto commented
K-mer repeat counting currently fails if a chromosome that is not present in the BAM file is requested, due to Pysam returning an exception. The script should check for available chromosomes before making a request to pysam.AlignmentFile.fetch()