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

bedmap produces incorrect results #43

Closed
sjneph opened this issue Feb 22, 2014 · 1 comment
Closed

bedmap produces incorrect results #43

sjneph opened this issue Feb 22, 2014 · 1 comment
Assignees

Comments

@sjneph
Copy link
Collaborator

sjneph commented Feb 22, 2014

There is a current scenario where bedmap will produce incorrect results.

foo1.bed:
chr1 1 200 a 1
chr1 10 20 b 3
chr1 50 150 c 4

foo2.bed:
chr1 1 200 a 1

Without --faster, bedmap should always produce correct results even with nested elements. The current output from:
bedmap --count --bp-ovr 11 foo1.bed foo2.bed

is:
1
0
0

when it should be:
1
0
1

Row two of the reference cannot have qualifying overlaps as its length is 10, and I used --bp-ovr 11. If --bp-ovr is 10 or less, the output is correct. If the file order is reversed, the output is correct. This needs to be fixed asap.

@sjneph sjneph self-assigned this Feb 22, 2014
@sjneph
Copy link
Collaborator Author

sjneph commented Feb 22, 2014

Fixed in commit 2cda151

@sjneph sjneph closed this as completed Feb 22, 2014
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant