#SHELL=qrsh hostname; /bin/zsh
#SHELL=dish -s herbie 2 mulligan 2
shell=bash

MISC=../../scripts

AUTOSOMES=1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
CHROMOSOMES=$(AUTOSOMES) X 
ALIGN=/space/james/mm8_17way
REF=mm8
SPECIES=mm8,hg18,panTro2,rheMac2,rn4,canFam2,bosTau2
REORDER=3,0,1,2,4,5,6

REGIONS=$(CHROMOSOMES)

MODEL=test_1.063
MODSPEC=tree_pruned:N=10,D=0.01
MODORDER=2

chromlen.txt : /depot/data1/cache/chrominfo.txt 
	grep $(REF) $< | tr " " "," | prefix_lines "mouse," | ./nib_location.py /depot/data1/cache/mouse/mm8/seq/%s.nib > $@

preCRM_regions.bed :
	( for f in ~/work/ross_mm7_predictions/sequences/*_mm7.txt; do head -1 $$f | egrep "chr[0-9XY]+\:[0-9]+\-[0-9]+" -o | tr ":-" "\t\t" | paste - <( basename $${f%%locus_hg17.txt} ); done ) > $@

%.scores.bz2 : %.bed
	prefix_lines "mm7." < $< \
	| ~/projects/bx-python-trunk/scripts/maf_tile_2.py $(SPECIES) chromlen.txt $(ALIGN)/chr*.maf \
	| $(MISC)/maf_translate_chars.py \
	| rp_score_maf.py /dev/stdin $(MODEL).sm /dev/stdout -m $(MODEL).mapping -w 100 -s 1 -M $(MODSPEC) -r $(REORDER) \
	| bzip2 \
	> $@

include ../common.mk