#!/bin/bash
#
#SBATCH --output=/home/jvierstra/proj/footprinting.2019/processed/CD19+-DS17186/motif.per-nt/logs/%A.%a.out
#SBATCH --mem=12G
#SBATCH --cpus-per-task=4

set -o pipefail

TMPDIR=`mktemp -d`

ARGS=(`cat /home/jvierstra/proj/footprinting.2019/processed/CD19+-DS17186/motif.per-nt/slurm.motif.params | head -n ${SLURM_ARRAY_TASK_ID} | tail -n 1`)

mkfifo ${TMPDIR}/motif.bed
mkfifo ${TMPDIR}/fps.bed

tabix /home/jvierstra/data/motifs/fimo/GRCh38_no_alts.alldbs.1e-4/fimo.combined.1e-4.parsed.bed.gz ${ARGS[0]} | awk '$5<1e-5' > ${TMPDIR}/motif.bed &
tabix /home/jvierstra/proj/footprinting.2019/processed/CD19+-DS17186/interval.all.bedgraph.gz ${ARGS[0]} | cut -f1-3,7 > ${TMPDIR}/fps.bed &

bedmap --delim "\t" --echo --indicator ${TMPDIR}/fps.bed ${TMPDIR}/motif.bed 	| awk '$4>0' | grep -v "NAN" | cut -f4- | sort -k1,1g > /home/jvierstra/proj/footprinting.2019/processed/CD19+-DS17186/motif.per-nt/motif.${ARGS[0]}

rm -rf ${TMPDIR}