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

Reset seqlet coordinates with respect to "--window" #20

Open
jvierstra opened this issue May 3, 2023 · 1 comment
Open

Reset seqlet coordinates with respect to "--window" #20

jvierstra opened this issue May 3, 2023 · 1 comment

Comments

@jvierstra
Copy link

The seqlet coordinates (per example) are returned relative to the "trimmed" contributions For example, if the user specificies a window of 100 with the "--window" argument, on a full sequence input width of 100 then the return seqlet coordinates, then the seqlet coordinates should be shifted by 450 bp when saving to the file MODISCO output file.

Alternatively, this could be just documented. I was trying to use the H5 file (after backwards conversion) to call seqlet hits on the examples and this messed me up for a couple of hours.

@snaqvi1990
Copy link

Could you clarify how you get a 450 bp shift when setting --window to 100? Doesn't it mean that there is a 100bp window centered on the peak center, and then all seqlet coordinates per example are with respect to that 100bp window? I am also trying to call seqlet hits and having issues. Thanks

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants