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

Turn off reverse complements #56

Closed
Turakar opened this issue Dec 2, 2024 · 3 comments · May be fixed by #57
Closed

Turn off reverse complements #56

Turakar opened this issue Dec 2, 2024 · 3 comments · May be fixed by #57

Comments

@Turakar
Copy link

Turakar commented Dec 2, 2024

It would be useful if a user could pass revcomp=False to disable reverse complementing for matching. In my use case, reverse complements should not match each other, but be considered separate. Could this be added to the library?

It seems like updating _align_patterns() should be enough?

@jmschrei
Copy link
Owner

jmschrei commented Dec 2, 2024

Conceptually, it could be done. However, I think this is built-in to many methods, the most important being the seqlet similarity scores, and I don't have the bandwidth to add this feature. Sorry about that.

@jmschrei jmschrei closed this as completed Dec 2, 2024
@Turakar
Copy link
Author

Turakar commented Dec 3, 2024

Not considering possible performance improvements - is it wrong to only change _align_patterns()?

@jmschrei
Copy link
Owner

jmschrei commented Dec 3, 2024

Yeah, you'd need too also change the seqlet similarity score in at least https://github.com/jmschrei/tfmodisco-lite/blob/main/modiscolite/affinitymat.py#L67. I don't remember all the places it would matter -- it's been a while since I've looked at the code.

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

Successfully merging a pull request may close this issue.

2 participants