So I’ve been trying to bgzip around 100 VCF files in parallel, but although the jobs are submitted and files get created there’s definitely something wrong.
So far I’ve been trying:
$for file in *.vcf > do > bsub /foo/bar/bgzip $file >> $file.gz
What is the correct way to do this?
Thanks in advance!
Advertisement
Answer
If, as my websearches suggest, bgzip works the same way as gzip, you don’t specify an output file:
for f in *.vcf; do bsub bgzip "$f" done
You should also double-quote variables in case your filenames have spaces in them.