code: def _print_header(data):
      
      print("##fileformat=VCFv4.2", file=STDOUT, end="")
      print("##source=seqbuster2.3", file=STDOUT, end="")
      print("##reference=mirbase", file=STDOUT, end="")
      for pos in data:
          print("##contig=<ID=%s>" % pos["chrom"], file=STDOUT, end="")
      print('##INFO=<ID=ID,Number=1,Type=String,Description="miRNA name">', file=STDOUT, end="")
      print('##FORMAT=<ID=GT,Number=1,Type=Integer,Description="Genotype">', file=STDOUT, end="")
      print('##FORMAT=<ID=NR,Number=A,Type=Integer,Description="Total reads supporting the variant">', file=STDOUT, end="")
      print('##FORMAT=<ID=NS,Number=A,Type=Float,Description="Total number of different sequences supporting the variant">', file=STDOUT, end="")
      print("#CHROM\tPOS\tID\tREF\tALT\tQUAL\tFILTER\tINFO\tFORMAT\tSAMP001", file=STDOUT, end="")
docs: Create vcf header to make
      a valid vcf.