feat(sort): per-dataset minimizer cutoff exponent

packages/nextclade-schemas/input-pathogen-json.schema.json

@@ -464,6 +464,14 @@
           "items": {
             "type": "string"
           }
+        },
+        "cutoff": {
+          "description": "Minimizer acceptance threshold given as an exponent: the retained fraction of k-mers is `2^cutoff / 2^32`. When absent, a default cutoff is used.",
+          "type": [
+            "integer",
+            "null"
+          ],
+          "format": "int64"
         }
       }
     },

packages/nextclade-schemas/input-pathogen-json.schema.yaml

@@ -292,6 +292,12 @@ definitions:
         type: array
         items:
           type: string
+      cutoff:
+        description: 'Minimizer acceptance threshold given as an exponent: the retained fraction of k-mers is `2^cutoff / 2^32`. When absent, a default cutoff is used.'
+        type:
+        - integer
+        - 'null'
+        format: int64
   LabelledMutationsConfig:
     description: Mapping from specific mutations to human-readable labels (e.g. clade-defining or drug-resistance mutations).
     examples:

packages/nextclade-schemas/nextclade-auspice-extensions.schema.json

@@ -798,6 +798,14 @@
           "items": {
             "type": "string"
           }
+        },
+        "cutoff": {
+          "description": "Minimizer acceptance threshold given as an exponent: the retained fraction of k-mers is `2^cutoff / 2^32`. When absent, a default cutoff is used.",
+          "type": [
+            "integer",
+            "null"
+          ],
+          "format": "int64"
         }
       }
     },

packages/nextclade-schemas/nextclade-auspice-extensions.schema.yaml

@@ -511,6 +511,12 @@ definitions:
         type: array
         items:
           type: string
+      cutoff:
+        description: 'Minimizer acceptance threshold given as an exponent: the retained fraction of k-mers is `2^cutoff / 2^32`. When absent, a default cutoff is used.'
+        type:
+        - integer
+        - 'null'
+        format: int64
   LabelledMutationsConfig:
     description: Mapping from specific mutations to human-readable labels (e.g. clade-defining or drug-resistance mutations).
     examples:

packages/nextclade/src/analyze/virus_properties.rs

@@ -150,6 +150,10 @@ pub struct MinimizerIndexConfig {
   #[serde(default, skip_serializing_if = "Vec::is_empty")]
   pub references: Vec<String>,
 
+  /// Minimizer acceptance threshold given as an exponent: the retained fraction of k-mers is `2^cutoff / 2^32`. When absent, a default cutoff is used.
+  #[serde(default, skip_serializing_if = "Option::is_none")]
+  pub cutoff: Option<i64>,
+
   #[serde(flatten)]
   pub other: serde_json::Value,
 }
@@ -482,12 +486,13 @@ mod tests {
 
   #[test]
   fn test_virus_properties_parses_minimizer_index() {
-    let json = r#"{"schemaVersion":"3.0.0","minimizerIndex":{"references":["minimizer_refs/additional_refs.fasta"]}}"#;
+    let json = r#"{"schemaVersion":"3.0.0","minimizerIndex":{"references":["minimizer_refs/additional_refs.fasta"],"cutoff":31}}"#;
     let props = VirusProperties::from_str(&json).unwrap();
     let mi = props
       .minimizer_index
       .expect("minimizerIndex should parse into the typed field");
     assert_eq!(vec_of_owned!["minimizer_refs/additional_refs.fasta"], mi.references);
+    assert_eq!(Some(31), mi.cutoff);
     assert!(props.other.get("minimizerIndex").is_none());
   }