#!/bin/bash
#
#SBATCH --output=/home/jvierstra/proj/footprinting.2019/processed/fBrain-DS23849/logs/%J.out
#SBATCH --mem=8G
#SBATCH --cpus-per-task=1
#SBATCH --partition=queue1 

set -e -o pipefail

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

thresholds=(0.2 0.1 0.05 0.01 0.001 0.0001)
thresh=${thresholds[${SLURM_ARRAY_TASK_ID}-1]}

cat /home/jvierstra/proj/footprinting.2019/processed/fBrain-DS23849/interval.all.bedgraph 	| awk -v OFS="\t" -v thresh="${thresh}" '$8 <= thresh { print $1, $2-3, $3+3; }' 	| sort-bed --tmpdir ${TMPDIR} --max-mem 6G - 	| bedops -m - 	| awk -v OFS="\t" -v thresh="${thresh}" '{ $4="."; $5=thresh; print; }' > /home/jvierstra/proj/footprinting.2019/processed/fBrain-DS23849/interval.all.fps.${thresh}.bed

bgzip -c /home/jvierstra/proj/footprinting.2019/processed/fBrain-DS23849/interval.all.fps.${thresh}.bed > /home/jvierstra/proj/footprinting.2019/processed/fBrain-DS23849/interval.all.fps.${thresh}.bed.gz
tabix -0 -p bed /home/jvierstra/proj/footprinting.2019/processed/fBrain-DS23849/interval.all.fps.${thresh}.bed.gz

rm -rf ${TMPDIR}