|
|
|
|
|
|
|
|
import "./tasks/fastp.wdl" as fastp |
|
|
|
|
|
|
|
|
#import "./tasks/fastp.wdl" as fastp |
|
|
import "./tasks/fastqTobam.wdl" as fastqTobam |
|
|
import "./tasks/fastqTobam.wdl" as fastqTobam |
|
|
import "./tasks/qualimap.wdl" as qualimap |
|
|
import "./tasks/qualimap.wdl" as qualimap |
|
|
import "./tasks/BQSR.wdl" as Sentieon_BQSR |
|
|
import "./tasks/BQSR.wdl" as Sentieon_BQSR |
|
|
|
|
|
|
|
|
#config |
|
|
#config |
|
|
String cluster_config |
|
|
String cluster_config |
|
|
String disk_size |
|
|
String disk_size |
|
|
Boolean fastqc |
|
|
|
|
|
|
|
|
Boolean? fastqc |
|
|
Boolean bamqc |
|
|
Boolean bamqc |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
if (normal_in1 != ""){ |
|
|
if (normal_in1 != ""){ |
|
|
|
|
|
|
|
|
if (fastqc){ |
|
|
|
|
|
call fastp.fastp as fastp_normal{ |
|
|
|
|
|
input: |
|
|
|
|
|
in1=normal_in1, |
|
|
|
|
|
in2=normal_in2, |
|
|
|
|
|
sample_id=sample_id+'_N', |
|
|
|
|
|
docker=docker_fastp, |
|
|
|
|
|
cluster_config=cluster_config, |
|
|
|
|
|
disk_size=disk_size |
|
|
|
|
|
|
|
|
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
if (!fastqc){ |
|
|
|
|
|
call fastqTobam.SentieonFastqToBam as SentieonFastqToBam_normal{ |
|
|
call fastqTobam.SentieonFastqToBam as SentieonFastqToBam_normal{ |
|
|
input: |
|
|
input: |
|
|
fastq1=normal_in1, |
|
|
fastq1=normal_in1, |
|
|
|
|
|
|
|
|
} |
|
|
} |
|
|
|
|
|
|
|
|
} |
|
|
} |
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
if (tumor_in1 != ""){ |
|
|
if (tumor_in1 != ""){ |
|
|
|
|
|
|
|
|
if (fastqc){ |
|
|
|
|
|
call fastp.fastp as fastp_tumor{ |
|
|
|
|
|
input: |
|
|
|
|
|
in1=tumor_in1, |
|
|
|
|
|
in2=tumor_in2, |
|
|
|
|
|
sample_id=sample_id+'_T', |
|
|
|
|
|
docker=docker_fastp, |
|
|
|
|
|
cluster_config=cluster_config, |
|
|
|
|
|
disk_size=disk_size |
|
|
|
|
|
|
|
|
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
if (!fastqc){ |
|
|
|
|
|
call fastqTobam.SentieonFastqToBam as SentieonFastqToBam_tumor{ |
|
|
call fastqTobam.SentieonFastqToBam as SentieonFastqToBam_tumor{ |
|
|
input: |
|
|
input: |
|
|
fastq1=tumor_in1, |
|
|
fastq1=tumor_in1, |
|
|
|
|
|
|
|
|
} |
|
|
} |
|
|
|
|
|
|
|
|
} |
|
|
} |
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
} |
|
|
|
|
|
|
|
|
|
|
|
call TNseq.sentieon_TNseq as sentieon_TNseq{ |
|
|
call TNseq.sentieon_TNseq as sentieon_TNseq{ |
|
|
input: |
|
|
input: |