소스 검색

first commit

master
YaqingLiu 5 년 전
부모
커밋
574caef601
3개의 변경된 파일28개의 추가작업 그리고 25개의 파일을 삭제
  1. +1
    -1
      tasks/BQSR.wdl
  2. +9
    -4
      tasks/PON.wdl
  3. +18
    -20
      workflow.wdl

+ 1
- 1
tasks/BQSR.wdl 파일 보기

@@ -30,7 +30,7 @@ task BQSR {
${SENTIEON_INSTALL_DIR}/bin/sentieon driver -t $nt --algo QualCal --plot --before ${sample}_recal_data.table --after ${sample}_recal_data.table.post ${sample}_recal_data.csv
${SENTIEON_INSTALL_DIR}/bin/sentieon plot bqsr -o ${sample}_bqsrreport.pdf ${sample}_recal_data.csv
${SENTIEON_INSTALL_DIR}/bin/sentieon plot QualCal -o ${sample}_bqsrreport.pdf ${sample}_recal_data.csv
>>>

+ 9
- 4
tasks/PON.wdl 파일 보기

@@ -3,6 +3,7 @@ task PON {
String SENTIEON_INSTALL_DIR
String SENTIEON_LICENSE
String sample
Boolean PONmode

File ref_dir
String fasta
@@ -22,10 +23,14 @@ task PON {
set -e
export SENTIEON_LICENSE=${SENTIEON_LICENSE}
nt=$(nproc)
mkdir -p /cromwell_root/tmp/cosmic/
cp ${cosmic_dir}/${cosmic_vcf} /cromwell_root/tmp/cosmic/
${SENTIEON_INSTALL_DIR}/bin/sentieon util vcfindex /cromwell_root/tmp/cosmic/${cosmic_vcf}
${SENTIEON_INSTALL_DIR}/bin/sentieon driver -t $nt -r ${ref_dir}/${fasta} -i ${normal_recaled_bam} --algo TNhaplotyper --detect_pon --cosmic /cromwell_root/tmp/cosmic/${cosmic_vcf} --dbsnp ${dbsnp_dir}/${dbsnp} ${sample}_pon.vcf
if [ ${PONmode} ]; then
mkdir -p /cromwell_root/tmp/cosmic/
cp ${cosmic_dir}/${cosmic_vcf} /cromwell_root/tmp/cosmic/
${SENTIEON_INSTALL_DIR}/bin/sentieon util vcfindex /cromwell_root/tmp/cosmic/${cosmic_vcf}
${SENTIEON_INSTALL_DIR}/bin/sentieon driver -t $nt -r ${ref_dir}/${fasta} -i ${normal_recaled_bam} --algo TNhaplotyper --detect_pon --cosmic /cromwell_root/tmp/cosmic/${cosmic_vcf} --dbsnp ${dbsnp_dir}/${dbsnp} ${sample}_pon.vcf
else
touch ${sample}_pon.vcf
fi
>>>
runtime {

+ 18
- 20
workflow.wdl 파일 보기

@@ -232,32 +232,30 @@ workflow {{ project_name }} {
cluster_config=cluster_config
}

if (PONmode) {
call PON.PON as PON {
input:
SENTIEON_INSTALL_DIR=SENTIEON_INSTALL_DIR,
SENTIEON_LICENSE=SENTIEON_LICENSE,
sample=sample[0],
fasta=fasta,
ref_dir=ref_dir,
normal_recaled_bam=normal_BQSR.recaled_bam,
normal_recaled_bam_index=normal_BQSR.recaled_bam_index,
cosmic_vcf=cosmic_vcf,
cosmic_dir=cosmic_dir,
dbsnp=dbsnp,
dbsnp_dir=dbsnp_dir,
docker=sentieon_docker,
disk_size=disk_size,
cluster_config=cluster_config
}
call PON.PON as PON {
input:
SENTIEON_INSTALL_DIR=SENTIEON_INSTALL_DIR,
SENTIEON_LICENSE=SENTIEON_LICENSE,
sample=sample[0],
PONmode=PONmode,
fasta=fasta,
ref_dir=ref_dir,
normal_recaled_bam=normal_BQSR.recaled_bam,
normal_recaled_bam_index=normal_BQSR.recaled_bam_index,
cosmic_vcf=cosmic_vcf,
cosmic_dir=cosmic_dir,
dbsnp=dbsnp,
dbsnp_dir=dbsnp_dir,
docker=sentieon_docker,
disk_size=disk_size,
cluster_config=cluster_config
}
File pon_vcf=PON.pon_vcf
}

call bcftools.bcftools as bcftools {
input:
PONmode=PONmode,
pon_vcf=pon_vcf,
pon_vcf=PON.pon_vcf,
docker=bcftools_docker,
disk_size=disk_size,
cluster_config=cluster_config

Loading…
취소
저장