code: def build_vcf_parts(feature, genome_2bit, info=None):
      
      base1 = genome_2bit[feature.chrom1].get(
          feature.start1, feature.start1 + 1).upper()
      id1 = "hydra{0}a".format(feature.name)
      base2 = genome_2bit[feature.chrom2].get(
          feature.start2, feature.start2 + 1).upper()
      id2 = "hydra{0}b".format(feature.name)
      orientation = _breakend_orientation(feature.strand1, feature.strand2)
      return (VcfLine(feature.chrom1, feature.start1, id1, base1,
                      _vcf_alt(base1, feature.chrom2, feature.start2,
                               orientation.is_rc1, orientation.is_first1),
                      _vcf_info(feature.start1, feature.end1, id2, info)),
              VcfLine(feature.chrom2, feature.start2, id2, base2,
                      _vcf_alt(base2, feature.chrom1, feature.start1,
                               orientation.is_rc2, orientation.is_first2),
                      _vcf_info(feature.start2, feature.end2, id1, info)))
docs: Convert BedPe feature information into VCF part representation.
  
      Each feature will have two VCF lines for each side of the breakpoint.