Skip to content

Commit

Permalink
Update LaunchMetrics.py
Browse files Browse the repository at this point in the history
pointing script to new methylated tables for hg38 and grcm39 for dragen 4.2 added id02 and id03 dragen servers back to production.
  • Loading branch information
darrelln32 committed Apr 19, 2024
1 parent 766ee25 commit 13a544a
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions scripts/LaunchMetrics.py
Original file line number Diff line number Diff line change
Expand Up @@ -16,7 +16,7 @@
# Global Variable : we do not want to process these experiments in this script
DO_NOT_PROCESS = ["DLP"]
# These recipes will be evaluated using DRAGEN because of their larger size of fastqs
RUN_ON_DRAGEN = ["MissionBio", "SingleCellCNV", "MouseWholeGenome", "HumanWholeGenome", "PombeWholeGenome", "ChIPSeq", "AmpliconSeq"]
RUN_ON_DRAGEN = ["MissionBio", "SingleCellCNV", "MouseWholeGenome", "HumanWholeGenome", "PombeWholeGenome", "ChIPSeq", "AmpliconSeq", "MethylCaptureSeq"]
# these projects willl only need demux stats
DEMUX_ONLY = ["SMARTSeq", "10X_Genomics"]

Expand Down Expand Up @@ -218,13 +218,13 @@ def dragen_methylation(sample, run, sample_parameters, work_directory, dragen_di

# get the correct path for the reference
if (sample_parameters["GTAG"] == "GRCh38"):
dragen_path = "/igo/work/igo/dragen_hash_tables/hg38_methylated"
dragen_path = "/igo/work/igo/dragen_hash_tables/4.2/hg38_methylated"
else:
dragen_path = "/igo/work/igo/dragen_hash_tables/grcm39_methylated"
dragen_path = "/igo/work/igo/dragen_hash_tables/4.2/grcm39_methylated"

metric_file_prefix = "{}___P{}___{}___{}".format(run, sample.project[8:], sample.sample_id, sample_parameters["GTAG"])
launch_dragen_methylation = "/opt/edico/bin/dragen --enable-methylation-calling true --methylation-protocol directional --ref-dir {} --fastq-list {} --fastq-list-sample-id {} --intermediate-results-dir /staging/temp --output-directory {} --output-file-prefix {} --enable-sort true --enable-duplicate-marking true".format(dragen_path, fastq_list, sample.sample_id, dragen_directory, sample.sample_id)
bsub_launch_dragen = "bsub -J {0}{1} -o {0}{1}.out -cwd \"{2}\" -m \"id03\" -q dragen -n 48 -M 4 {3}".format(dragen_methylation_job_name_header, sample.sample_id, dragen_directory, launch_dragen_methylation)
bsub_launch_dragen = "bsub -J {0}{1} -o {0}{1}.out -cwd \"{2}\" -m \"id01 id02 id03\" -q dragen -n 48 -M 4 {3}".format(dragen_methylation_job_name_header, sample.sample_id, dragen_directory, launch_dragen_methylation)
print(bsub_launch_dragen)
call(bsub_launch_dragen, shell = True)

Expand Down

0 comments on commit 13a544a

Please sign in to comment.