+ my $ignore_indels = $$opts{snps_only} ? 1 : 0;
+ my $ignore_snps = $$opts{indels_only} ? 1 : 0;
+ my $keep_ref = $$opts{keep_ref} ? 1 : 0;
+ my $title = exists($$opts{title}) ? $$opts{title} : 'data';
+
+ print $fh_out
+ qq[##fileformat=VCFv3.3\n],
+ qq[##INFO=DP,1,Integer,"Total Depth"\n],
+ qq[##FORMAT=GT,1,String,"Genotype"\n],
+ qq[##FORMAT=GQ,1,Integer,"Genotype Quality"\n],
+ qq[##FORMAT=DP,1,Integer,"Read Depth"\n],
+ qq[#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\t$title\n]
+ ;