code: def hydra_to_vcf_writer(hydra_file, genome_2bit, options, out_handle):
      
      _write_vcf_header(out_handle)
      brends = list(_get_vcf_breakends(hydra_file, genome_2bit, options))
      brends.sort(key=attrgetter("chrom", "pos"))
      for brend in brends:
          _write_vcf_breakend(brend, out_handle)
docs: Write hydra output as sorted VCF file.
  
      Requires loading the hydra file into memory to perform sorting
      on output VCF. Could generalize this to no sorting or by-chromosome
      approach if this proves too memory intensive.