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

Two modes --fast (default) and --sensitive in strobealign? #135

Open
ksahlin opened this issue Nov 1, 2022 · 0 comments
Open

Two modes --fast (default) and --sensitive in strobealign? #135

ksahlin opened this issue Nov 1, 2022 · 0 comments
Labels
feature New feature or request

Comments

@ksahlin
Copy link
Owner

ksahlin commented Nov 1, 2022

Seeding
At the seeding step, the current algorithm selects, for each syncmer (start strobe), a downstream strobe from a set of downstream syncmers in a window. It is obvious that there will be some redundancies in downstream strobe sampling (see attached figure for an example).

We could instead for each syncmer start strobe, select a k-mer from the window of k-mers downstream in "nucleotide space" (i.e. not syncmer-space).

The envisioned outcome is improved randomness in seeds, hence number of matches, hence accuracy.. The tradeoff is a slowdown in the seeding step as the windows have to be larger to have the same range in nucleotide space (in expectation 5x more candidate second strobes when syncmers are sampled with a density of 1/5).

Candidate site testing

We could in sensitive mode also increase the number of candidate sites tested with -M 40 instead of default -M 20.

Current strobe creation

See some redundancies in the second strobe sampling below.
image

@ksahlin ksahlin added the feature New feature or request label Nov 1, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
feature New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant