code: def run_radia_with_merge(job, rna_bam, tumor_bam, normal_bam, univ_options, radia_options):
      
      spawn = job.wrapJobFn(run_radia, rna_bam['rna_genome'], tumor_bam,
                            normal_bam, univ_options, radia_options, disk='100M',
                            memory='100M').encapsulate()
      merge = job.wrapJobFn(merge_perchrom_vcfs, spawn.rv(), univ_options, disk='100M', memory='100M')
      job.addChild(spawn)
      spawn.addChild(merge)
      return merge.rv()
docs: A wrapper for the the entire RADIA sub-graph.
  
      :param dict rna_bam: Dict dicts of bam and bai for tumor RNA-Seq obtained by running STAR within
             ProTECT.
      :param dict tumor_bam: Dict of bam and bai for tumor DNA-Seq
      :param dict normal_bam: Dict of bam and bai for normal DNA-Seq
      :param dict univ_options: Dict of universal options used by almost all tools
      :param dict radia_options: Options specific to RADIA
      :return: fsID to the merged RADIA calls
      :rtype: toil.fileStore.FileID