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

Open Syncmer? #1

Open
jianshu93 opened this issue Jun 26, 2024 · 2 comments
Open

Open Syncmer? #1

jianshu93 opened this issue Jun 26, 2024 · 2 comments

Comments

@jianshu93
Copy link

jianshu93 commented Jun 26, 2024

Hi Team,

Just wonder whether open syncmer could be another option (https://github.com/bluenote-1577/local-kmer-selection-results/blob/db4c5c686433fb09f8571da7ca11b355aaa28e30/src/seeding_methods.rs#L45). It was not bit packed there but only a sample reference implementation.

Thanks,

Jianshu

@imartayan
Copy link
Collaborator

Hi @jianshu93,

From what I understand, this corresponds to selecting windows where the minimizer is at position t, am I correct?
Should we treat differently the case where a single minimizer appears at multiple positions? Or should we only consider the leftmost position?

@jianshu93 jianshu93 changed the title minmer? Open Syncmer? Jun 27, 2024
@jianshu93
Copy link
Author

jianshu93 commented Jun 27, 2024

This is my understanding and the local kmer theory paper explain it very well, for any kmer, the smallest s mer (s<k) by some order equal to some position t in the kmer, this kmer will be selected, otherwise not selected. In practice, t is alway the middle position of the kmer. see figure 1 in this paper(Theory of local k-mer selection with applications to long-read alignment).

For repetitive kmers in a sequence, all will be selected since those kmers are independent, so it is somehow slower than minimizer, but offer more information for sequence with more mutations than minimizer.

No window is needed, but scan all kmers in a sequence.

Thanks,

Jianshu

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