diff --git a/src/subcommand/giraffe_main.cpp b/src/subcommand/giraffe_main.cpp index ce2e7f40825..b440cb73dfe 100644 --- a/src/subcommand/giraffe_main.cpp +++ b/src/subcommand/giraffe_main.cpp @@ -859,8 +859,8 @@ int main_giraffe(int argc, char** argv) { // Use downsampling instead of max unique minimizer count .add_entry("max-min", 0) .add_entry("num-bp-per-min", 1000) - .add_entry("downsample-window-count", 125) - .add_entry("downsample-window-length", 120) + .add_entry("downsample-window-count", 120) + .add_entry("downsample-window-length", 125) // Don't use the hit-cap||score-fraction filter because it doesn't do anything after downsampling .add_entry("hit-cap", 0) .add_entry("score-fraction", 1.0) @@ -909,14 +909,14 @@ int main_giraffe(int argc, char** argv) { .add_entry("watchdog-timeout", 30) .add_entry("batch-size", 10) // Use downsampling instead of max unique minimizer count - .add_entry("max-min", 100) - .add_entry("num-bp-per-min", 500) - .add_entry("downsample-window-count", 500) - .add_entry("downsample-window-length", 20) + .add_entry("max-min", 79) + .add_entry("num-bp-per-min", 152) + .add_entry("downsample-window-count", 15) + .add_entry("downsample-window-length", 227) // Don't use the hit-cap||score-fraction filter because it doesn't do anything after downsampling .add_entry("hit-cap", 0) .add_entry("score-fraction", 1.0) - .add_entry("hard-hit-cap", 20000) + .add_entry("hard-hit-cap", 13614) .add_entry("mapq-score-scale", 1) .add_entry("mapq-score-window", 150) .add_entry("zipcode-tree-score-threshold", 100.0) @@ -940,10 +940,10 @@ int main_giraffe(int argc, char** argv) { .add_entry("min-chaining-problems", 6) .add_entry("max-chaining-problems", std::numeric_limits::max()) .add_entry("max-lookback-bases", 20000) - .add_entry("max-lookback-bases-per-base", 0.15) + .add_entry("max-lookback-bases-per-base", 0.10501002120802233) .add_entry("item-bonus", 20) .add_entry("item-scale", 1) - .add_entry("gap-scale", 2.75) + .add_entry("gap-scale", 0.06759721757973396) .add_entry("max-indel-bases", 5000) .add_entry("max-indel-bases-per-base", 2.45) .add_entry("chain-score-threshold", 100.0)