Conversation
MicahR-Y
left a comment
There was a problem hiding this comment.
@yueyaog I made some suggestions for this code. The two major points are changing the inputs to an array to allow for more than 2 samples to be cross-checked at once and including an option for not failing the workflow if the samples don't match.
| String output_name | ||
| String picard_docker = "us.gcr.io/broad-gotc-prod/genomes-in-the-cloud:2.4.1-1540490856" | ||
| } | ||
|
|
There was a problem hiding this comment.
| Int mismatch_rc = if fail_on_mismatch then 1 else 2 |
|
|
||
| workflow CrosscheckFingerprintsWorkflow { | ||
| input { | ||
| File input_bamA |
There was a problem hiding this comment.
| File input_bamA |
| workflow CrosscheckFingerprintsWorkflow { | ||
| input { | ||
| File input_bamA | ||
| File input_bamB |
There was a problem hiding this comment.
| File input_bamB | |
| Array[File] input_bams |
| File input_bamB | ||
| File hapmap | ||
| File reference_fasta | ||
| String output_name |
There was a problem hiding this comment.
| String output_name | |
| String output_name | |
| Boolean fail_on_mismatch = false |
|
|
||
| call CrosscheckFingerprints { | ||
| input: | ||
| input_bamA = input_bamA, |
There was a problem hiding this comment.
| input_bamA = input_bamA, |
| I=~{input_bamB} \ | ||
| HAPLOTYPE_MAP=~{hapmap} \ | ||
| CROSSCHECK_BY=SAMPLE \ | ||
| O=~{output_name}.crosscheck_metrics \ |
There was a problem hiding this comment.
| O=~{output_name}.crosscheck_metrics \ | |
| O=~{output_name}.crosscheck_metrics.txt \ |
| HAPLOTYPE_MAP=~{hapmap} \ | ||
| CROSSCHECK_BY=SAMPLE \ | ||
| O=~{output_name}.crosscheck_metrics \ | ||
| REFERENCE_SEQUENCE=~{reference} \ |
There was a problem hiding this comment.
| REFERENCE_SEQUENCE=~{reference} \ | |
| REFERENCE_SEQUENCE=~{reference} \ | |
| EXIT_CODE_WHEN_MISMATCH=~{mismatch_rc} \ | |
| ~{extra_args} \ |
| } | ||
|
|
||
| output { | ||
| File crosscheckmetrics = "~{output_name}.crosscheck_metrics" |
There was a problem hiding this comment.
| File crosscheckmetrics = "~{output_name}.crosscheck_metrics" | |
| File crosscheckmetrics = "~{output_name}.crosscheck_metrics.txt" | |
| String crosscheck_results = if read_int("rc") == 0 then "PASS" else "FAIL" |
| File crosscheckmetrics = "~{output_name}.crosscheck_metrics" | ||
| } | ||
| } | ||
| # |
There was a problem hiding this comment.
Can you remove the commented out task here?
| output_name = output_name, | ||
| picard_docker = picard_docker | ||
| } | ||
| # call FingerprintsResult { |
There was a problem hiding this comment.
Can you remove the commented out call here?
|
|
||
| runtime { | ||
| docker: picard_docker | ||
| disks: "local-disk " + disk_gb + " HDD" |
There was a problem hiding this comment.
| disks: "local-disk " + disk_gb + " HDD" | |
| disks: "local-disk " + disk_gb + " HDD" | |
| continueOnReturnCode: [0, 2] |
@MicahR-Y Here is the original crosscheck FP WDL. Feel free to suggest changes or send me the updated version you have.