|
|
@@ -29,7 +29,7 @@ task cnvkit{ |
|
|
|
/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 |
|
|
|
#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 call ${sample}.segmetrics.cns --drop-low-coverage --filter ci -m threshold --purity $PURITY -o ${sample}.call.cns |