When I iterate over the compressed fastq file using pysam, I only get 26,107 reads instead of the expected 1,009,470 reads. In [1]: import pysam In [2]: len([x for x in pysam.FastxFile('example.fastq.gz')]) Out[2]: 26107 I've traced the issue back to bgzip not decompress...