"{{ project_name }}.ballgown_docker": "{{ ballgown_docker }}", | "{{ project_name }}.ballgown_docker": "{{ ballgown_docker }}", | ||||
"{{ project_name }}.ballgown_cluster": "{{ ballgown_cluster }}", | "{{ project_name }}.ballgown_cluster": "{{ ballgown_cluster }}", | ||||
"{{ project_name }}.count_docker": "{{ count_docker }}", | "{{ project_name }}.count_docker": "{{ count_docker }}", | ||||
"{{ project_name }}.count_cluster": "{{ count_cluster }}", | |||||
"{{ project_name }}.count_length": "{{ count_length }}", | |||||
"{{ project_name }}.count_cluster": "{{ count_cluster }}", | |||||
"{{ project_name }}.count_length": "{{ count_length }}", | |||||
"{{ project_name }}.pre_alignment_qc": "{{ pre_alignment_qc }}" | "{{ project_name }}.pre_alignment_qc": "{{ pre_alignment_qc }}" | ||||
} | } |
runtime { | runtime { | ||||
docker:docker | docker:docker | ||||
cluster: cluster | |||||
systemDisk: "cloud_ssd 40" | |||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | |||||
cluster: cluster | |||||
systemDisk: "cloud_ssd 40" | |||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | |||||
} | } | ||||
output { | output { | ||||
File read1_html = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_fastqc.html") | File read1_html = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_fastqc.html") |
nt=$(nproc) | nt=$(nproc) | ||||
# mkdir -p /cromwell_root/tmp | # mkdir -p /cromwell_root/tmp | ||||
# cp -r ${screen_ref_dir} /cromwell_root/tmp/ | # cp -r ${screen_ref_dir} /cromwell_root/tmp/ | ||||
sed -i "s#/cromwell_root/fastq_screen_reference#${screen_ref_dir}#g" ${fastq_screen_conf} | |||||
sed -i "s#/cromwell_root/fastq_screen_reference#${screen_ref_dir}#g" ${fastq_screen_conf} | |||||
fastq_screen --aligner bowtie2 --conf ${fastq_screen_conf} --top 100000 --threads $nt ${read1} | fastq_screen --aligner bowtie2 --conf ${fastq_screen_conf} --top 100000 --threads $nt ${read1} | ||||
fastq_screen --aligner bowtie2 --conf ${fastq_screen_conf} --top 100000 --threads $nt ${read2} | fastq_screen --aligner bowtie2 --conf ${fastq_screen_conf} --top 100000 --threads $nt ${read2} | ||||
>>> | >>> | ||||
runtime { | runtime { | ||||
docker:docker | docker:docker | ||||
cluster: cluster | |||||
systemDisk: "cloud_ssd 40" | |||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | |||||
cluster: cluster | |||||
systemDisk: "cloud_ssd 40" | |||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | |||||
} | } | ||||
output { | output { | ||||
File png1 = "${read1name}_screen.png" | File png1 = "${read1name}_screen.png" |
command { | command { | ||||
nt=$(nproc) | |||||
hisat2 -t -p $nt -x ${idx}/${idx_prefix} -1 ${read_1P} -2 ${read_2P} -S ${base}.sam --un-conc-gz ${base}_un.fq.gz | |||||
nt=$(nproc) | |||||
hisat2 -t -p $nt -x ${idx}/${idx_prefix} -1 ${read_1P} -2 ${read_2P} -S ${base}.sam --un-conc-gz ${base}_un.fq.gz | |||||
} | } | ||||
runtime { | runtime { | ||||
docker: docker | |||||
cluster: cluster | |||||
systemDisk: "cloud_ssd 40" | |||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | |||||
docker: docker | |||||
cluster: cluster | |||||
systemDisk: "cloud_ssd 40" | |||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | |||||
} | } | ||||
output { | output { |
nt=$(nproc) | nt=$(nproc) | ||||
/opt/qualimap/qualimap rnaseq -bam ${bam} -outformat HTML -outdir ${bamname}_rnaseq -gtf ${gtf} -pe --java-mem-size=10G | /opt/qualimap/qualimap rnaseq -bam ${bam} -outformat HTML -outdir ${bamname}_rnaseq -gtf ${gtf} -pe --java-mem-size=10G | ||||
tar -zcvf ${bamname}_rnaseq_qualimap.zip ${bamname}_rnaseq | tar -zcvf ${bamname}_rnaseq_qualimap.zip ${bamname}_rnaseq | ||||
/opt/qualimap/qualimap bamqc -bam ${bam} -outformat PDF:HTML -nt $nt -outdir ${bamname}_bamqc --java-mem-size=32G | |||||
/opt/qualimap/qualimap bamqc -bam ${bam} -outformat PDF:HTML -nt $nt -outdir ${bamname}_bamqc --java-mem-size=32G | |||||
tar -zcvf ${bamname}_bamqc_qualimap.zip ${bamname}_bamqc | tar -zcvf ${bamname}_bamqc_qualimap.zip ${bamname}_bamqc | ||||
>>> | >>> | ||||
output { | output { | ||||
File rnaseq_zip = "${bamname}_rnaseq_qualimap.zip" | File rnaseq_zip = "${bamname}_rnaseq_qualimap.zip" | ||||
File bamqc_zip = "${bamname}_bamqc_qualimap.zip" | |||||
File bamqc_zip = "${bamname}_bamqc_qualimap.zip" | |||||
} | } | ||||
} | } |
File read1 | File read1 | ||||
File read2 | File read2 | ||||
File idx | |||||
File idx | |||||
File screen_ref_dir | File screen_ref_dir | ||||
File fastq_screen_conf | File fastq_screen_conf | ||||
String idx_prefix | String idx_prefix | ||||
File gtf | |||||
File gtf | |||||
String disk_size | String disk_size | ||||
String fastqc_docker | |||||
String fastqc_cluster | |||||
String fastqc_docker | |||||
String fastqc_cluster | |||||
String fastqscreen_docker | String fastqscreen_docker | ||||
String fastqscreen_cluster | |||||
String fastqscreen_cluster | |||||
String hisat2_docker | String hisat2_docker | ||||
String hisat2_cluster | String hisat2_cluster | ||||
String stringtie_docker | String stringtie_docker | ||||
disk_size= disk_size | disk_size= disk_size | ||||
} | } | ||||
call samtools.samtools as samtools { | |||||
input: | |||||
call samtools.samtools as samtools { | |||||
input: | |||||
docker = samtools_docker, | docker = samtools_docker, | ||||
cluster = samtools_cluster, | cluster = samtools_cluster, | ||||
sam = hisat2.sam, | sam = hisat2.sam, | ||||
disk_size = disk_size | disk_size = disk_size | ||||
} | } | ||||
call stringtie.stringtie as stringtie { | |||||
input: | |||||
call stringtie.stringtie as stringtie { | |||||
input: | |||||
docker = stringtie_docker, | docker = stringtie_docker, | ||||
cluster = stringtie_cluster, | cluster = stringtie_cluster, | ||||
gtf = gtf, | gtf = gtf, |