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.

53 lines
1.1KB

  1. task charger{
  2. String sample_id
  3. File vep_vcf
  4. File inheritanceGeneList
  5. File clinvar
  6. File pathogenic
  7. File HotSpot3D
  8. File PP2GeneList
  9. String disk_size
  10. String docker
  11. String cluster
  12. command <<<
  13. set -o pipefail
  14. set -e
  15. vep_input=`echo ${vep_vcf}|awk -F "/" '{print $NF}'`
  16. vep_index=`echo ${vep_vcf}|awk -F "/" '{print $NF}'|sed 's/.vcf//'`
  17. echo ${sample_id}
  18. mkdir ./input/
  19. mkdir ./output/
  20. grep "#" ${vep_vcf} > ./input/$vep_input
  21. grep -v "#" ${vep_vcf}|sed 's/\.0//' >> ./input/$vep_input
  22. charger -f ./input/$vep_input \
  23. -o ./output/$vep_index.charger.tsv \
  24. -D ${inheritanceGeneList} \
  25. -l --mac-clinvar-tsv ${clinvar} \
  26. -z ${pathogenic} \
  27. -H ${HotSpot3D} \
  28. --PP2GeneList ${PP2GeneList}
  29. >>>
  30. runtime {
  31. docker: docker
  32. cluster: cluster
  33. systemDisk: "cloud_ssd 40"
  34. dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  35. }
  36. output {
  37. Array[File] filter_result=glob("./output/${sample_id}*")
  38. }
  39. }