|
|
@@ -36,8 +36,8 @@ task mendelian { |
|
|
|
rtg index -f vcf ${family_name}.consistent.vcf.gz |
|
|
|
|
|
|
|
zcat ${child_vcf_gz} | grep '##' > child_header |
|
|
|
zcat ${family_name}.violation.vcf.gz | grep -v '##' | cut -f1-10 | awk '$10!~/\./' | cat header - > ${child_name}.violation.vcf |
|
|
|
zcat ${family_name}.consistent.vcf.gz | grep -v '##' | cut -f1-10 | awk '$10!~/\./' | cat header - > ${child_name}.consistent.vcf |
|
|
|
zcat ${family_name}.violation.vcf.gz | grep -v '##' | cut -f1-10 | awk '$10!~/\./' | cat child_header - > ${child_name}.violation.vcf |
|
|
|
zcat ${family_name}.consistent.vcf.gz | grep -v '##' | cut -f1-10 | awk '$10!~/\./' | cat child_header - > ${child_name}.consistent.vcf |
|
|
|
|
|
|
|
rtg bgzip ${child_name}.violation.vcf |
|
|
|
rtg index -f vcf ${child_name}.violation.vcf.gz |