Fix edge case where no reads map for a species

biomedicalinformaticsgroup / Sargasso

Sargasso disambiguates mixed-species high-throughput sequencing data.

Other

8 stars 4 forks source link

Closed lweasel closed 5 years ago

lweasel commented 6 years ago

Causes an exception in filter_sample_reads.py. This used to work (by fluke I think), but something has changed in htslib.

hxin commented 5 years ago

Did we fix this already when I hit the bug where one of the species bam files is empty? That is the function hits_generator in samutils.