@@ -30,7 +30,6 @@ task cnvkit{ | |||
# Derive each segment's absolute integer copy number, ploidy must be int value | |||
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 call ${sample}.segmetrics.cns --drop-low-coverage --filter ci -m threshold --purity $PURITY -o ${sample}.call.cns | |||