let vaxrank
?(configuration = Tools.Vaxrank.Configuration.default)
vcfs bam predictor alleles =
fun ~(var_count : int) ->
let vcfs_compiled = List.map vcfs ~f:(fun v -> v ~var_count) in
let bam_compiled = bam ~var_count in
function_call "vaxrank" ([
"configuration", string Tools.Vaxrank.Configuration.(name configuration);
"alleles", alleles ~var_count;
"predictor", string Hla_utils.(predictor_to_string predictor);
"bam", bam_compiled;
] @ (List.mapi ~f:(fun i v -> (sprintf "vcf%d" i, v)) vcfs_compiled))