@@ -185,9 +185,25 @@ task SortBed {
185
185
186
186
command <<<
187
187
set -eo pipefail
188
- echo " Sort bed file"
188
+
189
+ echo " [DEBUG] Current working directory:"
190
+ pwd
191
+ echo " [DEBUG] Directory contents:"
192
+ ls -la
193
+ echo " [DEBUG] Environment variables:"
194
+ env
195
+
196
+ echo " [DEBUG] Checking GATK installation:"
197
+ which gatk || echo " GATK not found"
198
+ gatk --help || echo " GATK command failed"
199
+
200
+ echo " [DEBUG] Sort bed file"
189
201
sort -k1,1V -k2,2n -k3,3n " ~{unsorted_bed}" > sorted.bed
190
- echo " Transform bed file to intervals list with Picard----------------------------------------"
202
+
203
+ echo " [DEBUG] Sorted bed contents:"
204
+ cat sorted.bed
205
+
206
+ echo " [DEBUG] Running BedToIntervalList"
191
207
gatk --java-options " -Xms4g" \
192
208
BedToIntervalList \
193
209
-I sorted.bed \
@@ -224,12 +240,28 @@ task SamToFastq {
224
240
225
241
command <<<
226
242
set -eo pipefail
243
+
244
+ echo " [DEBUG] Current working directory:"
245
+ pwd
246
+ echo " [DEBUG] Directory contents:"
247
+ ls -la
248
+ echo " [DEBUG] Environment variables:"
249
+ env
250
+
251
+ echo " [DEBUG] Checking GATK installation:"
252
+ which gatk || echo " GATK not found"
253
+ gatk --help || echo " GATK command failed"
254
+
255
+ echo " [DEBUG] Input BAM header:"
256
+ samtools view -H " ~{input_bam}"
257
+
258
+ echo " [DEBUG] Running SamToFastq"
227
259
gatk --java-options " -Dsamjdk.compression_level=5 -Xms4g" \
228
260
SamToFastq \
229
261
--INPUT " ~{input_bam}" \
230
262
--FASTQ " ~{base_file_name}.fastq" \
231
263
--INTERLEAVE true \
232
- --INCLUDE_NON_PF_READS true
264
+ --INCLUDE_NON_PF_READS true
233
265
>>>
234
266
235
267
output {
0 commit comments