#!/bin/bash
declare -A MCR

MCR[ ITR31-wt3P90.bam ]=37447047
MCR[ ITR32-wtpool190514.bam ]=49251794

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


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