Browse Source

更新 'tasks/CNVkit.wdl'

master
meng 2 years ago
parent
commit
97a83e54b1
1 changed files with 3 additions and 2 deletions
  1. +3
    -2
      tasks/CNVkit.wdl

+ 3
- 2
tasks/CNVkit.wdl View File

File ref_dir File ref_dir
String? fasta String? fasta
File ref_flat File ref_flat
File hrd
#File hrd
File tumor_bam File tumor_bam
File tumor_bam_index File tumor_bam_index
File? normal_bam File? normal_bam
/root/miniconda2/bin/cnvkit.py metrics ${sample}.cnr -s ${sample}.cns > ${sample}.stats /root/miniconda2/bin/cnvkit.py metrics ${sample}.cnr -s ${sample}.cns > ${sample}.stats
# Derive each segment's absolute integer copy number, ploidy must be int value # Derive each segment's absolute integer copy number, ploidy must be int value
PURITY=`awk -F'\t' '{print $6}' ${hrd} | sed -n '2p'`
#PURITY=`awk -F'\t' '{print $6}' ${hrd} | sed -n '2p'`
PURITY=3.3
/root/miniconda2/bin/cnvkit.py segmetrics ${sample}.cnr -s ${sample}.cns --ci -o ${sample}.segmetrics.cns /root/miniconda2/bin/cnvkit.py segmetrics ${sample}.cnr -s ${sample}.cns --ci -o ${sample}.segmetrics.cns
/root/miniconda2/bin/cnvkit.py call ${sample}.segmetrics.cns --drop-low-coverage --filter ci -m threshold --purity $PURITY -o ${sample}.call.cns /root/miniconda2/bin/cnvkit.py call ${sample}.segmetrics.cns --drop-low-coverage --filter ci -m threshold --purity $PURITY -o ${sample}.call.cns

Loading…
Cancel
Save