Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
apraga opened this issue Feb 26, 2023 · 0 comments
Open

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

apraga opened this issue Feb 26, 2023 · 0 comments
Labels
enhancement New feature or request

Comments

@apraga
Copy link

apraga commented Feb 26, 2023

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,
@kescobo kescobo added the enhancement New feature or request label Feb 26, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants