Skip to content

Commit

Permalink
temporary fix for GTF file error
Browse files Browse the repository at this point in the history
  • Loading branch information
dmcmanam committed Apr 3, 2023
1 parent dbad5a2 commit 2fef272
Showing 1 changed file with 2 additions and 1 deletion.
3 changes: 2 additions & 1 deletion scripts/alignment_and_picard.py
Original file line number Diff line number Diff line change
Expand Up @@ -258,7 +258,8 @@ def rna_alignment_and_metrics(sample, run, sample_parameters, rna_directory, wor
metric_file = "{}___P{}___{}___{}".format(run, prjct, sample.sample_id, sample_parameters["GTAG"])
fastq_list = "/igo/staging/FASTQ/{}/Reports/fastq_list.csv ".format(run)

launch_dragen_rna = "/opt/edico/bin/dragen -f -r {} --fastq-list {} --fastq-list-sample-id {} -a {} --intermediate-results-dir /staging/temp --enable-map-align true --enable-sort true --enable-bam-indexing true --enable-map-align-output true --output-format BAM --enable-rna true --enable-duplicate-marking true --enable-rna-quantification true --output-file-prefix {} --output-directory {}".format(rna_path, fastq_list, sample.sample_id, sample_parameters["GTF"], metric_file, rna_directory)
gtf_file = "/igo/work/genomes/H.sapiens/hg38/gencode.v40.chr_patch_hapl_scaff.annotation.gtf"
launch_dragen_rna = "/opt/edico/bin/dragen -f -r {} --fastq-list {} --fastq-list-sample-id {} -a {} --intermediate-results-dir /staging/temp --enable-map-align true --enable-sort true --enable-bam-indexing true --enable-map-align-output true --output-format BAM --enable-rna true --enable-duplicate-marking true --enable-rna-quantification true --output-file-prefix {} --output-directory {}".format(rna_path, fastq_list, sample.sample_id, gtf_file, metric_file, rna_directory)
bsub_launch_dragen_rna = "bsub -J {0}{1} -o {0}{1}.out -cwd \"{2}\" -m \"id01 id02 id03\" -q dragen -n 48 -M 4 {3}".format(rna_dragen_job_name_header, sample.sample_id, rna_directory, launch_dragen_rna)
print(bsub_launch_dragen_rna)
call(bsub_launch_dragen_rna, shell = True)
Expand Down

0 comments on commit 2fef272

Please sign in to comment.