code: def hydra_breakpoints(in_bam, pair_stats):
      
      in_bed = convert_bam_to_bed(in_bam)
      if os.path.getsize(in_bed) > 0:
          pair_bed = pair_discordants(in_bed, pair_stats)
          dedup_bed = dedup_discordants(pair_bed)
          return run_hydra(dedup_bed, pair_stats)
      else:
          return None
docs: Detect structural variation breakpoints with hydra.