Skip to content
Snippets Groups Projects
Commit fe773ff8 authored by Brandi Cantarel's avatar Brandi Cantarel
Browse files

add Callset to vcfheader union.sh

parent 0eb2c5f6
No related merge requests found
......@@ -48,6 +48,8 @@ for i in ${dir}/*.vcf.gz; do
done
echo "##fileformat=VCFv4.2" > header.vcf
echo "##INFO=<ID=CallSet,Number='.',Type=String,Description=\"Caller\">" >> header.vcf
echo "##INFO=<ID=CallSetInconsistent,Number=1,Type=Integer,Description=\"Caller Inconsistency Manual Inspection\">" >> header.vcf
zcat ${dir}/*.vcf.gz |grep "##" |grep -v '#fileformat' |sort -u |grep 'ALT\|FILTER\|FORMAT\|INFO' >> header.vcf
perl $baseDir/unionvcf.pl header.vcf $list2
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment