vcf header

tasks/mendelian.wdl

@@ -27,7 +27,7 @@ task mendelian {
 
 		zcat ${family_name}.family.anno.mendelian.vcf.gz | grep '#' > header
 		zcat ${family_name}.family.anno.mendelian.vcf.gz | grep -v '#'| grep 'MCU\|MCV' | cat header -  > ${family_name}.violation.vcf
-		zcat ${family_name}.family.anno.mendelian.vcf.gz | grep -v '#'| grep -v 'MCV' | grep -v 'MCU' > ${family_name}.consistent.vcf
+		zcat ${family_name}.family.anno.mendelian.vcf.gz | grep -v '#'| grep -v 'MCV' | grep -v 'MCU' | cat header - > ${family_name}.consistent.vcf
 
 		rtg bgzip ${family_name}.violation.vcf 
 		rtg index -f vcf ${family_name}.violation.vcf.gz