millironx/XAM.jl
1
0
Fork 0
mirror of https://github.com/MillironX/XAM.jl.git synced 2024-11-14 22:33:14 +00:00
Code Issues Releases Packages Activity

Add tryread! check

Browse source
This commit is contained in:
Ciarán O’Mara 2021-04-02 01:00:24 +11:00
parent 0e6ef62e3e
commit 742face886
1 changed files with 4 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

5
src/bam/overlap.jl
View file

@ -63,7 +63,10 @@ function Base.iterate(iter::OverlapIterator, state)
while state.chunkid ≤ lastindex(state.chunks) while state.chunkid ≤ lastindex(state.chunks)
chunk = state.chunks[state.chunkid] chunk = state.chunks[state.chunkid]
while BGZFStreams.virtualoffset(iter.reader.stream) < chunk.stop while BGZFStreams.virtualoffset(iter.reader.stream) < chunk.stop
read!(iter.reader, state.record) if BioGenerics.IO.tryread!(iter.reader, state.record) === nothing
# no more overlapping records in this chunk since the end of the chunk was reached.
break
end
c = compare_intervals(state.record, (state.refindex, iter.interval)) c = compare_intervals(state.record, (state.refindex, iter.interval))
if c == 0 if c == 0
return copy(state.record), state return copy(state.record), state