BioJulia/XAM.jl

Add eachoverlap(reader::Reader, refname::String)

Opened this issue · 0 comments

Hi,

Would it be possible to add eachoverlap such that it can be browsed by chromosome when the BAM is indexed ?

Expected Behavior

chrom = "chr1"
itr = eachoverlap(reader, chrom) 

Current Behavior

chrom = "chr1"
idx = findfirst(isequal(chrom), reader.refseqnames)
refseqlen = reader.refseqlens[idx]

itr = eachoverlap(reader, chrom, 1:refseqlen))

Thanks,