#!/usr/bin/env zsh

################################################################
export PATH=~/local/bin:$PATH
export PATH=~/.local/bin:$PATH

export PYTHONPATH=$PYTHONPATH:~/.local/lib/python2.7/site-packages
export PYTHONPATH=$PYTHONPATH:~/local/lib/python2.7/site-packages

### Job name
#BSUB -J PeakCalling
 
### File / path where STDOUT & STDERR will be written
#BSUB -o PeakCalling
  
### Request the time you need for execution in minutes
### The format for the parameter is: [hour:]minute,
### that means for 80 minutes you could also use this: 1:20
#BSUB -W 5:00
   
### Request memory you need for your job in TOTAL in MB
#BSUB -M 10240

### Make sure the jobs can use the HPCWORK file system
#BSUB -R "select[hpcwork]"

### Assign the job to the specified project
#BSUB -P lect0023

mkdir -p ~/Peaks
target=/hpcwork/lect0023/bioinfolab_2018/Data/ATAC-seq/H1hesc/ATAC.bam
macs2 callpeak -t $target -n H1hesc --outdir ~/Peaks -f BAM -g hs --keep-dup all --call-summits --nomodel