From 7b5cae51ea79e01f2411bcdc4c6a8c297ef96b14 Mon Sep 17 00:00:00 2001 From: ekiernan Date: Fri, 10 Jan 2025 15:41:01 -0500 Subject: [PATCH] testing gtf path --- pipelines/skylab/atac/atac.wdl | 18 +++++++++--------- 1 file changed, 9 insertions(+), 9 deletions(-) diff --git a/pipelines/skylab/atac/atac.wdl b/pipelines/skylab/atac/atac.wdl index e147a24a75..437a87fa5f 100644 --- a/pipelines/skylab/atac/atac.wdl +++ b/pipelines/skylab/atac/atac.wdl @@ -605,16 +605,16 @@ task CreateFragmentFile { # Add GTF to uns field # Original path from args.annotation_file - annotation_gtf = str(atac_gtf) # e.g., '/cromwell_root/gcp-public-data--broad-references/hg38/v0/star/v2_7_10a/modified_v43.annotation.gtf' - - # Transform the path - if annotation_gtf.startswith('/cromwell_root/'): - stripped_path = annotation_gtf[len('/cromwell_root/'):] # Remove '/cromwell_root/' - updated_path = f'gs://{stripped_path}' # Add 'gs://' prefix - else: - updated_path = str(atac_gtf) + gtf_path = str(~{annotations_gtf}) # e.g., '/cromwell_root/gcp-public-data--broad-references/hg38/v0/star/v2_7_10a/modified_v43.annotation.gtf' + + # # Transform the path + # if annotation_gtf.startswith('/cromwell_root/'): + # stripped_path = annotation_gtf[len('/cromwell_root/'):] # Remove '/cromwell_root/' + # updated_path = f'gs://{stripped_path}' # Add 'gs://' prefix + # else: + # updated_path = str(atac_gtf) - atac_data.uns["reference_gtf_file"] = updated_path + atac_data.uns["reference_gtf_file"] = gtf_path # calculate tsse metrics snap.metrics.tsse(atac_data, atac_gtf) # Write new atac file