code: def run_and_save(data):
      
      run(None, data)
      stats_file, idxstats_file = _get_stats_files(data)
      data = tz.update_in(data, ["depth", "samtools", "stats"], lambda x: stats_file)
      data = tz.update_in(data, ["depth", "samtools", "idxstats"], lambda x: idxstats_file)
      return data
docs: Run QC, saving file outputs in data dictionary.