Browse Source

Update 'tasks/fastqscreen.wdl'

output file names from fastq.gz and fq.gz
tags/v0.2.0
renluyao 3 years ago
parent
commit
dbacb202c7
1 changed files with 12 additions and 13 deletions
  1. +12
    -13
      tasks/fastqscreen.wdl

+ 12
- 13
tasks/fastqscreen.wdl View File

File read2 File read2
File screen_ref_dir File screen_ref_dir
File fastq_screen_conf File fastq_screen_conf
String read1name = basename(read1,".fastq.gz")
String read2name = basename(read2,".fastq.gz")
String docker String docker
String cluster_config String cluster_config
String disk_size String disk_size
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/
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} --subset 1000000 --threads $nt ${read1}
fastq_screen --aligner bowtie2 --conf ${fastq_screen_conf} --subset 1000000 --threads $nt ${read2}
>>> >>>


runtime { runtime {
docker:docker docker:docker
cluster: cluster_config
systemDisk: "cloud_ssd 40"
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
cluster: cluster_config
systemDisk: "cloud_ssd 40"
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
} }
output { output {
File png1 = "${read1name}_screen.png"
File txt1 = "${read1name}_screen.txt"
File html1 = "${read1name}_screen.html"
File png2 = "${read2name}_screen.png"
File txt2 = "${read2name}_screen.txt"
File html2 = "${read2name}_screen.html"
File png1 = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_screen.png")
File txt1 = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_screen.txt")
File html1 = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_screen.html")
File png2 = sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_screen.png")
File txt2 = sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_screen.txt")
File html2 = sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_screen.html")
} }
} }

Loading…
Cancel
Save