@@ -179,15 +179,15 @@ def dragen(sample, run, sample_parameters, work_directory, dragen_directory, fas
179
179
180
180
# get the correct path for the reference
181
181
if (sample_parameters ["GTAG" ] == "GRCh38" ) or (sample .genome == "Synthetic" ):
182
- dragen_path = "/igo/work/igo/dragen_hash_tables/4.2 /hg38-alt_masked.cnv.graph.hla.rna-9-r3 .0-1"
183
- vcfFileOption = "--qc-cross-cont-vcf /opt/edico/config /sample_cross_contamination_resource_hg38.vcf.gz"
182
+ dragen_path = "/igo/work/igo/dragen_hash_tables/4.3.6 /hg38-alt_masked.cnv.graph.hla.rna-10-r4 .0-1"
183
+ vcfFileOption = "--qc-cross-cont-vcf /opt/dragen/4.3.6/resources/qc /sample_cross_contamination_resource_hg38.vcf.gz"
184
184
else :
185
- dragen_path = "/igo/work/igo/dragen_hash_tables/4.2 /{}" .format (sample_parameters ["GTAG" ])
185
+ dragen_path = "/igo/work/igo/dragen_hash_tables/4.3.6 /{}" .format (sample_parameters ["GTAG" ])
186
186
vcfFileOption = ""
187
187
188
188
metric_file_prefix = "{}___P{}___{}___{}" .format (run , sample .project [8 :], sample .sample_id , sample_parameters ["GTAG" ])
189
- launch_dragen = "/opt/edico /bin/dragen --force --ref-dir {} --fastq-list {} --fastq-list-sample-id {} --intermediate-results-dir /staging/temp --output-directory {} --output-file-prefix {} {} --enable-duplicate-marking true --enable-sort true --enable-map-align true --enable-map-align-output true --output-format cram --bin_memory 70000000000 " .format (dragen_path , fastq_list , sample .sample_id , dragen_directory , sample .sample_id , vcfFileOption )
190
- bsub_launch_dragen = "bsub -J {0}{1} -o {0}{1}.out -cwd \" {2}\" -m \" id02 \" -q dragen -n48 -M4 {3}" .format (dragen_job_name_header , sample .sample_id , dragen_directory , launch_dragen )
189
+ launch_dragen = "/opt/dragen/4.3.6 /bin/dragen --force --ref-dir {} --fastq-list {} --fastq-list-sample-id {} --intermediate-results-dir /staging/temp --output-directory {} --output-file-prefix {} {} --enable-duplicate-marking true --enable-sort true --enable-map-align true --enable-map-align-output true --output-format cram --bin_memory 70000000000 " .format (dragen_path , fastq_list , sample .sample_id , dragen_directory , sample .sample_id , vcfFileOption )
190
+ bsub_launch_dragen = "bsub -J {0}{1} -o {0}{1}.out -cwd \" {2}\" -m \" id01 \" -q dragen -n48 -M4 {3}" .format (dragen_job_name_header , sample .sample_id , dragen_directory , launch_dragen )
191
191
print (bsub_launch_dragen )
192
192
call (bsub_launch_dragen , shell = True )
193
193
0 commit comments