-
Notifications
You must be signed in to change notification settings - Fork 6
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
Issue with coordinates in modkit pileup #135
Comments
Hello @ccastignani, Could you tell me how you're observing the off-by-one error? It may be an artifact of the viewer. For example, the BED specification is zero-based whereas many viewers (such as IGV) are one-based. This can make everything appear off-by-one. Also, with the |
Hi,
I'm running modkit pileup (v0.2.4) as follows:
Basecalling has been done using dorado basecaller (modified-bases 5mCG_5hmCG and [email protected]), and alignment has been done using minimap2.
Unfortunately, when looking at the output bed file from modkit pileup, there coordinates seem to be off consistently by one base pair upstream. The read counts regarding modified/canonical base pairs seem correct, it's just that they are all off by one position.
Is this something others have experienced before? (couldn't find any other issue about this). Any suggestions or recommendations to fix this would be highly appreciated!
Thanks :)
The text was updated successfully, but these errors were encountered: