eachoverlap running past end of stream

[rulixxx]

When looping trough a BAM file with eachoverlap I get an error if the selected interval intersects past the last read of the BAM file.

Expected Behavior

samtools view -b chr1:1-10000 bigBAM.bam > small.bam
samtools index small.bam

reader = open(BAM.Reader, "small.bam", index="small.bam.bai")
for record in eachoverlap(reader, "chr1", 1:1000000)
    # ...
    # ...
end
close(reader)

Loops correctly trough all the reads of the small bam file.

Current Behavior

The loop manages to read the last BAM record but then fails with a:

Stacktrace:
 [1] getindex at ./array.jl:728 [inlined]
 [2] virtualoffset(::BGZFStreams.BGZFStream{IOStream}) at /nfs/users/nfs_r/ra11/.julia/packages/BGZFStreams/qApsr/src/bgzfstream.jl:156
 [3] iterate(::BioAlignments.BAM.OverlapIterator{IOStream}, ::BioAlignments.BAM.OverlapIteratorState) at /nfs/users/nfs_r/ra11/.julia/packages/BioAlignments/FOisL/src/bam/overlap.jl:65
 [4] top-level scope at ./REPL[4]:2

Possible Solution / Implementation

Steps to Reproduce (for bugs)

Same as stated above. Take a section of a BAM and run

Context

Your Environment

• Package Version used: 1.0.1 & 2.0.0
• Julia Version used: 1.2.0
• Operating System and version (desktop or mobile): Linux
• Link to your project: