Skip to content

Instantly share code, notes, and snippets.

View icaoberg's full-sized avatar

icaoberg icaoberg

View GitHub Profile
@icaoberg
icaoberg / example02.sh
Last active February 3, 2020 23:20 — forked from arq5x/go.sh
[bedtools] Compute average scores for share intervals
# icaoberg - this example is fork that uses a bedtools in a Singularity container
CONTAINER=../../singularity-bedtools.simg
echo "chr1 10 50 10" > a.bed
echo "chr1 20 40 20" > b.bed
echo "chr1 30 33 30" > c.bed
# Find the sub-intervals shared and unique to each file.
if [ -f $CONTAINER ]; then
singularity run --app bedtools $CONTAINER multiinter -i a.bed b.bed c.bed | column -t
@icaoberg
icaoberg / .gitignore
Last active February 3, 2020 23:21 — forked from arq5x/make-master-hmm.sh
[bedtools] For Gemini: Create a master ChromHMM track from the 9 distinct cell types.
*.txt
*.bedg
#!/usr/bin/env python
import bz2
import datetime
import os
import sys
import time
import urllib2
import warnings