Add tests for issue #31

Test for expected results as described in https://github.com/BioJulia/XAM.jl/issues/31.

Project.toml

@@ -20,6 +20,7 @@ BGZFStreams = "0.3"
 BioAlignments = "2"
 BioGenerics = "0.1"
 BioSequences = "2.0.4"
+FormatSpecimens = "1.1"
 GenomicFeatures = "2"
 Indexes = "0.1"
 TranscodingStreams = "0.6, 0.7, 0.8, 0.9"

test/runtests.jl

@@ -25,3 +25,4 @@ end
 
 include("test_sam.jl")
 include("test_bam.jl")
+include("test_issues.jl")

test/test_issues.jl

@@ -0,0 +1,13 @@
+@testset "Issues" begin
+
+    # https://github.com/BioJulia/XAM.jl/issues/31
+    path_bam = joinpath(path_of_format("BAM"), "SRR7993829_1.100K.forward.bam")
+
+    open(BAM.Reader, path_bam, index = path_bam * ".bai") do reader
+
+        @test count(overlap -> true, eachoverlap(reader, "JH584304.1", 51000:51200)) == 0
+        @test count(overlap -> true, eachoverlap(reader, "JH584304.1", 51000:51715)) == 1
+
+    end
+
+end