renluyao
/
variant-calling
feito fork de YaqingLiu/variant-calling
Observar 1
Favorito 0
Fork 0
Código Issues 0 Pull requests 0 Versões 0 Wiki Atividade
Germline & Somatic short variant discovery (SNVs + Indels) for WGS & WES.
Você não pode selecionar mais de 25 tópicos Os tópicos devem começar com uma letra ou um número, podem incluir traços ('-') e podem ter até 35 caracteres.
123 Commits
1 Branch
Tag: f02a98c73b
Branches Tags
${ item.name }
Criar branch ${ searchTerm }
de f02a98c73b
${ noResults }
variant-calling/tasks/somaticFilter.wdl

47 linhas
2.6KB
Original Anotar Histórico 

  1. task somaticFilter {

  2.   

  3.   String sample

  4.   File varscan_snp_somatic_hc

  5.   File varscan_snp_germline_hc

  6.   File varscan_snp_loh_hc

  7.   File varscan_indel_somatic_hc

  8.   File varscan_indel_germline_hc

  9.   File varscan_indel_loh_hc

  10.   File varscan_indel

  11. 

  12.   String docker

  13.   String cluster_config

  14.   String disk_size

  15. 

  16.   command <<<

  17.     java -Xmx32g -jar /opt/VarScan.v2.4.3.jar somaticFilter ${varscan_snp_somatic_hc} --min-coverage 10 --min-reads2 2 --min-strands2 1 --min-avg-qual 20 --p-value 0.1 --indel-file ${varscan_indel} --output-file ${sample}.VarScan.SNP.Somatic.filter.vcf

  18. 

  19.     java -Xmx32g -jar /opt/VarScan.v2.4.3.jar somaticFilter ${varscan_snp_germline_hc} --min-coverage 10 --min-reads2 2 --min-strands2 1 --min-avg-qual 20 --p-value 0.1 --indel-file ${varscan_indel} --output-file ${sample}.VarScan.SNP.Germline.filter.vcf

  20. 

  21.     java -Xmx32g -jar /opt/VarScan.v2.4.3.jar somaticFilter ${varscan_snp_loh_hc} --min-coverage 10 --min-reads2 2 --min-strands2 1 --min-avg-qual 20 --p-value 0.1 --indel-file ${varscan_indel} --output-file ${sample}.VarScan.SNP.LOH.filter.vcf

  22. 

  23.     java -Xmx32g -jar /opt/VarScan.v2.4.3.jar somaticFilter ${varscan_indel_somatic_hc} --min-coverage 10 --min-reads2 2 --min-strands2 1 --min-avg-qual 20 --p-value 0.1 --output-file ${sample}.VarScan.INDEL.Somatic.filter.vcf

  24. 

  25.     java -Xmx32g -jar /opt/VarScan.v2.4.3.jar somaticFilter ${varscan_indel_germline_hc} --min-coverage 10 --min-reads2 2 --min-strands2 1 --min-avg-qual 20 --p-value 0.1 --output-file ${sample}.VarScan.INDEL.Germline.filter.vcf

  26. 

  27.     java -Xmx32g -jar /opt/VarScan.v2.4.3.jar somaticFilter ${varscan_indel_loh_hc} --min-coverage 10 --min-reads2 2 --min-strands2 1 --min-avg-qual 20 --p-value 0.1 --output-file ${sample}.VarScan.INDEL.LOH.filter.vcf

  28. 

  29.     # Merge SNP and INDEL

  30.     awk '{if ($1!~/^#/) print}' ${sample}.VarScan.SNP.Somatic.filter.vcf | cat ${sample}.VarScan.INDEL.Somatic.filter.vcf - > ${sample}.VarScan.Somatic.filter.vcf

  31.     awk '{if ($1!~/^#/) print}' ${sample}.VarScan.SNP.Germline.filter.vcf | cat ${sample}.VarScan.INDEL.Germline.filter.vcf - > ${sample}.VarScan.Germline.filter.vcf

  32.     awk '{if ($1!~/^#/) print}' ${sample}.VarScan.SNP.LOH.filter.vcf | cat ${sample}.VarScan.INDEL.LOH.filter.vcf - > ${sample}.VarScan.LOH.filter.vcf

  33.   >>>

  34. 

  35.   runtime {

  36.     docker: docker

  37.     cluster: cluster_config

  38.     systemDisk: "cloud_ssd 40"

  39.     dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"

  40.   }

  41. 

  42.   output {

  43. 		File varscan_somatic_filter = "${sample}.VarScan.Somatic.filter.vcf"

  44.     File varscan_germline_filter = "${sample}.VarScan.Germline.filter.vcf"

  45. 		File varscan_loh_filter = "${sample}.VarScan.LOH.filter.vcf"

  46.   }

  47. }