|
|
@@ -23,13 +23,13 @@ task benchmark { |
|
|
|
/opt/rtg-tools/dist/rtg-tools-3.10.1-4d58ead/rtg index -f vcf ${sample}.rtg.vcf.gz |
|
|
|
|
|
|
|
if [ ${sample_mark} == "LCL5" ];then |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL5.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL5.consensus.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
elif [ ${sample_mark} == "LCL6" ]; then |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL6.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL6.consensus.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
elif [ ${sample_mark} == "LCL7" ]; then |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL7.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL7.consensus.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
elif [ ${sample_mark} == "LCL8" ]; then |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL8.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
/opt/hap.py/bin/hap.py ${benchmarking_dir}/LCL8.consensus.vcf.gz ${sample}.rtg.vcf.gz --threads $nt -o ${sample} |
|
|
|
else |
|
|
|
echo "only for quartet samples" |
|
|
|
fi |