File tree 1 file changed +3
-1
lines changed
1 file changed +3
-1
lines changed Original file line number Diff line number Diff line change @@ -312,7 +312,7 @@ task RegenotypeChunk {
312
312
elif [ ~ {use_delly} -eq 1 ]; then
313
313
${TIME_COMMAND} /delly call --svtype ALL --genome ~ {reference_fa} --vcffile ~ {vcf_to_genotype} --outfile genotypes.bcf ${LOCAL_FILE}
314
314
bcftools view --output-type v --output genotypes.vcf genotypes.bcf
315
- rm -f genotypes.bcf
315
+ # rm -f genotypes.bcf
316
316
fi
317
317
N_LINES=$( grep ' #' genotypes.vcf | wc -l)
318
318
rm -f $( basename ${FILE} ) && echo 0 || echo 1
@@ -340,12 +340,14 @@ task RegenotypeChunk {
340
340
mv tmp.txt genotypes.txt; rm -f tmp.txt
341
341
fi
342
342
head genotypes.txt
343
+ touch genotypes.bcf
343
344
done < ~ {chunk}
344
345
>>>
345
346
346
347
output {
347
348
File format = "format.txt"
348
349
File genotypes = "genotypes.txt"
350
+ File debug = "genotypes.bcf"
349
351
}
350
352
runtime {
351
353
docker : "fcunial/lr-genotyping"
You can’t perform that action at this time.
0 commit comments