#!/bin/bash
#
#SBATCH --output=/home/jvierstra/proj/footprinting.2019/processed/fHeart-DS24812/logs/%J.out
#SBATCH --mem=32G
#SBATCH --cpus-per-task=1
#SBATCH --partition=queue0

set -e -o pipefail

TMPDIR=/tmp/$SLURM_JOB_ID
mkdir -p ${TMPDIR}

sort-bed --tmpdir ${TMPDIR} --max-mem 30G /home/jvierstra/proj/footprinting.2019/processed/fHeart-DS24812/interval.all.bedgraph | starch - > /home/jvierstra/proj/footprinting.2019/processed/fHeart-DS24812/interval.all.bedgraph.starch

bgzip -c /home/jvierstra/proj/footprinting.2019/processed/fHeart-DS24812/interval.all.bedgraph > /home/jvierstra/proj/footprinting.2019/processed/fHeart-DS24812/interval.all.bedgraph.gz
tabix -0 -p bed /home/jvierstra/proj/footprinting.2019/processed/fHeart-DS24812/interval.all.bedgraph.gz

rm -rf ${TMPDIR}