#!/bin/bash


cat <<__SCRIPT__ > /home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/enrichment/run.sh
#!/bin/bash
#SBATCH --output=/home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/enrichment/logs/%A.%a.out
#SBATCH --mem=16G
#SBATCH --cpus-per-task=1
#SBATCH --partition=queue0
set -e -o pipefail

source /home/jvierstra/.bashrc

conda activate py3.9_default

MOTIF=\$(cat /home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/pfms.txt | head -n \${SLURM_ARRAY_TASK_ID} | tail -n 1)

python motif.py  \
	/home/jvierstra/public_html/projects/encode4-allelic-imbalance/allelic_mapping/01132022/all.hets.combined.pvals.txt.gz \
	/home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/by_pfm/\${MOTIF}.txt \
	/home/jvierstra/data/motifs/moods/GRCh38_no_alts.alldbs.1e-4/pfm/\${MOTIF}.pfm \
	/home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/enrichment/\${MOTIF}
__SCRIPT__


njobs=$(wc -l < /home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/pfms.txt)

JOB0=$(sbatch --export=ALL \
	--job-name=motif_enrich \
	--array=1-${njobs} \
	/home/jvierstra/public_html/projects/encode4-allelic-imbalance/motifs/02032022/enrichment/run.sh)
echo $JOB0