#!/bin/bash
set -o errexit
set -o nounset
set -o pipefail

WIG=$1
THRESHOLD=$2
CHROM_SIZES=$3

wiggletools write_bg - gt $THRESHOLD $WIG | awk '{ print $1, $2, $3 }' | tr ' ' '\t' | sort -k1,1 -k2,2n | bedtools merge -i - > $WIG.conserved.bed
wiggletools write_bg - lt -$THRESHOLD $WIG | awk '{ print $1, $2, $3 }' | tr ' ' '\t' | sort -k1,1 -k2,2n | bedtools merge -i - > $WIG.accelerated.bed
bedToBigBed $WIG.conserved.bed $CHROM_SIZES $WIG.conserved.bb
bedToBigBed $WIG.acclerated.bed $CHROM_SIZES $WIG.accelerated.bb