#!/bin/bash
declare -A MCR

MCR[ bcatenincon.bam ]=18646048
MCR[ bcatenindox.bam ]=31288136
MCR[ HnrnpKcon.bam ]=24402936
MCR[ HnrnpKdox.bam ]=9101043

let min=9999999999;
for i in b*bam
do
echo $i
if [ ${MCR[ $i ]} -le $min ]; then let min=${MCR[ $i ]}; fi
done
echo "Sampling to smallest library containing" $min " reads"

for i in b*bam
do
sc=$(echo "scale=15; $min/${MCR[ $i ]}" | bc) 
echo $i ${MCR[ $i ]} $sc
/data/results/tools/proton/tophat2unstranded-sampled-bigwig1.sh $i /data/results/reference/hg19/Homo_sapiens/UCSC/hg19/Sequence/WholeGenomeFasta/hsa.hg19.genome $sc
#/data/results/tools/proton/tophat2unstranded-sampled-bigwig1.sh $i /data/results/reference/hg19/Homo_sapiens/UCSC/hg19/Sequence/WholeGenomeFasta/hsa.hg19.genome $sc
done


let min=9999999999;
for i in H*bam
do
echo $i
if [ ${MCR[ $i ]} -le $min ]; then let min=${MCR[ $i ]}; fi
done
echo "Sampling to smallest library containing" $min " reads"

for i in H*bam
do
sc=$(echo "scale=15; $min/${MCR[ $i ]}" | bc) 
echo $i ${MCR[ $i ]} $sc
/data/results/tools/proton/tophat2unstranded-sampled-bigwig1.sh $i /data/results/reference/hg19/Homo_sapiens/UCSC/hg19/Sequence/WholeGenomeFasta/hsa.hg19.genome $sc
#/data/results/tools/proton/tophat2unstranded-sampled-bigwig1.sh $i /data/results/reference/hg19/Homo_sapiens/UCSC/hg19/Sequence/WholeGenomeFasta/hsa.hg19.genome $sc
done