-
Notifications
You must be signed in to change notification settings - Fork 61
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
feat: deprecate old alignment params #1270
Conversation
The latest updates on your projects. Learn more about Vercel for Git ↗︎
|
(o!("--min-seeds (minSeeds)"), &self.min_seeds), | ||
(o!("--seed-length (seedLength)"), &self.seed_length), | ||
(o!("--seed-spacing (seedSpacing)"), &self.seed_spacing), | ||
(o!("--min-match-rate (minMatchRate)"), &self.min_match_rate), | ||
(o!("--mismatches-allowed (mismatchesAllowed)"), &self.mismatches_allowed), | ||
(o!("--max-indel (maxIndel)"), &self.max_indel), |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
These are the params being deprecated.
And the error message is constructed just below this fragment.
Let's use `kmer_length` instead of removed `seed_length` and adjust the coefficient to roughly match the old value.
@@ -46,7 +46,7 @@ pub fn align_nuc( | |||
); | |||
} | |||
|
|||
if ref_len + qry_len < (10 * params.seed_length) { | |||
if ref_len + qry_len < (20 * params.kmer_length) { | |||
// for very short sequences, use full square |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
According to Richard's recommendation, for short sequence check let's use kmer_length
instead of removed seed_length
and adjust the coefficient to roughly match the old value.
Here I deprecate some of the alignment params which lost their use after the alignment algo was changed.
I emit an error if any of these parameters are still present in CLI args and/or in dataset's
pathogen.json
. We could also make it a warning.