let seq2hla fq =
let fastq = get_fastq fq in
let r1 = fastq#product#r1 in
let r2 =
match fastq#product#r2 with
| Some r -> r
| None ->
failf "Seq2HLA doesn't support Single_end_sample(s)."
in
let work_dir =
Config.work_dir //
sprintf "%s-%s_seq2hla-workdir"
fastq#product#escaped_sample_name
fastq#product#fragment_id_forced
in
Seq2hla_result (
Tools.Seq2HLA.hla_type
~work_dir ~run_with ~run_name:fastq#product#escaped_sample_name ~r1 ~r2
)