#!/bin/bash
#SBATCH --job-name=zSMM-33
#SBATCH --output=zSMM.33.out
#SBATCH --error=zSMM.33.err
#SBATCH --time=24:00:00
#SBATCH -p akundaje,sfgf
#SBATCH --cpus-per-task=1
#SBATCH --mem=50000
python ~/code/bisulfite/dSMF-footprints.py /oak/stanford/groups/akundaje/marinovg/ENCODE4/DNAzoo/2021-03-17-Novogene-NovaSeq/King_eccDNA_NOME.deduped.pooled.sorted.bam ~/genomes/hg20/sequence/hg38.fa GC MYC-ecDNA-regions.300bp-tiling-50bp-step.annotated.bed-chunks32 0 1 2 3 MYC-ecDNA-regions-300bp-tiling-50bp-step -label 3 -heatmap ~/code/plotting/heatmap.py 10 10 binary 10,100 -minCov 0.8 -unstranded