瀏覽代碼

fix(CNVkit): output

master
YaqingLiu 2 年之前
父節點
當前提交
28cf9d337f
共有 1 個檔案被更改,包括 4 行新增8 行删除
  1. +4
    -8
      tasks/CNVkit.wdl

+ 4
- 8
tasks/CNVkit.wdl 查看文件

set -e set -e
nt=$(nproc) nt=$(nproc)
CNV_ANALYSIS_PATH=/cromwell_root/tmp/cnvkit
mkdir -p $CNV_ANALYSIS_PATH
cd $CNV_ANALYSIS_PATH
cnvkit.py access ${ref_dir}/${fasta} -o access.bed cnvkit.py access ${ref_dir}/${fasta} -o access.bed
# Prepare the target bed # Prepare the target bed
cnvkit.py target ${regions} --annotate ${ref_flat} --split --short-names -o my_baits.bed cnvkit.py target ${regions} --annotate ${ref_flat} --split --short-names -o my_baits.bed
cnvkit.py genemetrics ${sample}.cnr -s ${sample}.call.cns -t 0.2 -m 3 -o ${sample}.segment-genes.txt cnvkit.py genemetrics ${sample}.cnr -s ${sample}.call.cns -t 0.2 -m 3 -o ${sample}.segment-genes.txt
cat ${sample}.ratio-genes.txt | tail -n+2 | cut -f1 | sort | uniq > ratio-genes.txt cat ${sample}.ratio-genes.txt | tail -n+2 | cut -f1 | sort | uniq > ratio-genes.txt
cat ${sample}.segment-genes.txt | tail -n+2 | cut -f1 | sort | uniq > segment-genes.txt cat ${sample}.segment-genes.txt | tail -n+2 | cut -f1 | sort | uniq > segment-genes.txt
comm -12 ratio-genes.txt segment-genes.txt > trusted_cnv_genes.txt
for gene in `cat trusted_cnv_genes.txt`
comm -12 ratio-genes.txt segment-genes.txt > ${sample}.trusted_cnv_genes.txt
for gene in `cat ${sample}.trusted_cnv_genes.txt`
do do
cnvkit.py scatter ${sample}.cnr -s ${sample}.call.cns -g $gene -o ./gainloss/${sample}.$gene.scatter.pdf cnvkit.py scatter ${sample}.cnr -s ${sample}.call.cns -g $gene -o ./gainloss/${sample}.$gene.scatter.pdf
done done
File stats = "${sample}.stats" File stats = "${sample}.stats"
File call_cns = "${sample}.call.cns" File call_cns = "${sample}.call.cns"
File genemetrics = "${sample}.cnv.txt" File genemetrics = "${sample}.cnv.txt"
File gainloss_genes = "/cromwell_root/tmp/cnvkit/trusted_cnv_genes.txt"
Array[File] gainloss = glob("/cromwell_root/tmp/cnvkit/gainloss/*")
File gainloss_genes = "${sample}.trusted_cnv_genes.txt"
Array[File] gainloss = glob("./gainloss/*")
} }
} }

Loading…
取消
儲存