Add eachoverlap(reader::Reader, refname::String)
Opened this issue · 0 comments
apraga commented
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,