室间质评
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

33 lines
936B

  1. task JI {
  2. File vcf_a
  3. File vcf_b
  4. File dir_name
  5. File sdf
  6. String name_a = basename(vcf_a,".vcf")
  7. String name_b = basename(vcf_b,".vcf")
  8. String docker
  9. String cluster_config
  10. String disk_size
  11. command <<<
  12. set -o pipefail
  13. set -e
  14. /opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg bgzip ${vcf_a} -c > ${name_a}.rtg.vcf.gz
  15. /opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg index -f vcf ${name_a}.rtg.vcf.gz
  16. /opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg bgzip ${vcf_a} -c > ${name_b}.rtg.vcf.gz
  17. /opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg index -f vcf ${name_b}.rtg.vcf.gz
  18. /opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg vcfeval -b ${name_a}.rtg.vcf.gz -c ${name_b}.rtg.vcf.gz -o ${dir_name} -t ${sdf}
  19. >>>
  20. runtime {
  21. docker:docker
  22. cluster: cluster_config
  23. systemDisk: "cloud_ssd 40"
  24. dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  25. }
  26. output {
  27. File JI_summary = "${dir_name}/summary.txt"
  28. }
  29. }