fix header vcf

locate_coding_snps.sh

@@ -10,7 +10,7 @@ intersect_vcf_annotation () {
 
   bedtools intersect -wb -a $2 -b $3 | grep $4 | rev | cut -f10- | rev | uniq > $1_$4.bed
   wait;
-  cat <(head $2) $1_$4.bed > $1_$4.vcf
+  cat <(grep "^#" $2) $1_$4.bed > $1_$4.vcf
 }
 
 # diplodus sargus