瀏覽代碼

revise fastqscreen

tags/v0.1.0
stead99 4 年之前
父節點
當前提交
0339137a26
共有 5 個文件被更改,包括 57 次插入57 次删除
  1. +2
    -2
      tasks/ballgown.wdl
  2. +4
    -4
      tasks/fastqc.wdl
  3. +9
    -9
      tasks/fastqscreen.wdl
  4. +3
    -3
      tasks/hisat2.wdl
  5. +39
    -39
      workflow.wdl

+ 2
- 2
tasks/ballgown.wdl 查看文件

task ballgown { task ballgown {
File gene_abundance File gene_abundance
String base = basename(gene_abundance, ".gene.abundance.txt")
String base=basename(gene_abundance, ".gene.abundance.txt")
Array[File] ballgown Array[File] ballgown
String docker String docker
String cluster String cluster
} }
output { output {
File mat_expression = "${base}.txt"
File mat_expression="${base}.txt"
} }
} }

+ 4
- 4
tasks/fastqc.wdl 查看文件

dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
} }
output { output {
File read1_html = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_fastqc.html")
File read1_zip = sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_fastqc.zip")
File read2_html = sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_fastqc.html")
File read2_zip = sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_fastqc.zip")
File read1_html=sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_fastqc.html")
File read1_zip=sub(basename(read1), "\\.(fastq|fq)\\.gz$", "_fastqc.zip")
File read2_html=sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_fastqc.html")
File read2_zip=sub(basename(read2), "\\.(fastq|fq)\\.gz$", "_fastqc.zip")
} }
} }

+ 9
- 9
tasks/fastqscreen.wdl 查看文件

set -o pipefail set -o pipefail
set -e set -e
nt=$(nproc) nt=$(nproc)
# mkdir -p /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}
mkdir -p /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}
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}
>>> >>>
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" 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="${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"
} }
} }

+ 3
- 3
tasks/hisat2.wdl 查看文件

} }


output { output {
File sam = base + ".sam"
File unmapread_1p = base + "_un.fq.1.gz"
File unmapread_2p = base + "_un.fq.2.gz"
File sam=base + ".sam"
File unmapread_1p=base + "_un.fq.1.gz"
File unmapread_2p=base + "_un.fq.2.gz"
} }
} }

+ 39
- 39
workflow.wdl 查看文件

input: input:
read1=read1, read1=read1,
read2=read2, read2=read2,
docker = fastqc_docker,
cluster = fastqc_cluster,
disk_size = disk_size
docker=fastqc_docker,
cluster=fastqc_cluster,
disk_size=disk_size
} }
call fastqscreen.fastqscreen as fastqscreen { call fastqscreen.fastqscreen as fastqscreen {
input: input:
read1 = read1,
read2 = read2,
docker = fastqscreen_docker,
cluster = fastqscreen_cluster,
screen_ref_dir = screen_ref_dir,
fastq_screen_conf = fastq_screen_conf,
disk_size = disk_size
read1=read1,
read2=read2,
docker=fastqscreen_docker,
cluster=fastqscreen_cluster,
screen_ref_dir=screen_ref_dir,
fastq_screen_conf=fastq_screen_conf,
disk_size=disk_size
} }
} }
call hisat2.hisat2 as hisat2 { call hisat2.hisat2 as hisat2 {
input: input:
docker = hisat2_docker,
cluster = hisat2_cluster,
docker=hisat2_docker,
cluster=hisat2_cluster,
idx=idx, idx=idx,
idx_prefix=idx_prefix, idx_prefix=idx_prefix,
read_1P=read1, read_1P=read1,
read_2P=read2, read_2P=read2,
disk_size= disk_size
disk_size=disk_size
} }


call samtools.samtools as samtools { call samtools.samtools as samtools {
input: input:
docker = samtools_docker,
cluster = samtools_cluster,
sam = hisat2.sam,
disk_size= disk_size
docker=samtools_docker,
cluster=samtools_cluster,
sam=hisat2.sam,
disk_size=disk_size
} }


call qualimap.qualimap as qualimap { call qualimap.qualimap as qualimap {
input: input:
bam = samtools.out_bam,
bai = samtools.out_bam_index,
gtf = gtf,
docker = qualimap_docker,
cluster = qualimap_cluster,
disk_size = disk_size
bam=samtools.out_bam,
bai=samtools.out_bam_index,
gtf=gtf,
docker=qualimap_docker,
cluster=qualimap_cluster,
disk_size=disk_size
} }


call stringtie.stringtie as stringtie { call stringtie.stringtie as stringtie {
input: input:
docker = stringtie_docker,
cluster = stringtie_cluster,
gtf = gtf,
bam = samtools.out_bam,
disk_size = disk_size
docker=stringtie_docker,
cluster=stringtie_cluster,
gtf=gtf,
bam=samtools.out_bam,
disk_size=disk_size
} }


call ballgown.ballgown as ballgown { call ballgown.ballgown as ballgown {
input: input:
docker = ballgown_docker,
cluster = ballgown_cluster,
ballgown = stringtie.ballgown,
gene_abundance = stringtie.gene_abundance,
disk_size = disk_size
docker=ballgown_docker,
cluster=ballgown_cluster,
ballgown=stringtie.ballgown,
gene_abundance=stringtie.gene_abundance,
disk_size=disk_size
} }
call count.count as count { call count.count as count {
input: input:
sample_id = sample_id,
docker = count_docker,
cluster = count_cluster,
ballgown = stringtie.ballgown,
disk_size = disk_size,
count_length = count_length
sample_id=sample_id,
docker=count_docker,
cluster=count_cluster,
ballgown=stringtie.ballgown,
disk_size=disk_size,
count_length=count_length
} }
} }

Loading…
取消
儲存