# bash script for scaling theta bed file using watterson's formula

numChroms=$1

origBed=preprocessedDrosophilaData/theta-1K-unscaled.bed
scaledBed=preprocessedDrosophilaData/theta-1K.bed

cat $origBed | awk -v numChroms=$numChroms -v OFS="\t" \
   'BEGIN{wattersonsA=0; for(n=1;n<numChroms;n++) {wattersonsA+=(1.0/n);}} \
    { \
       print $1,$2,$3,$4/wattersonsA,$6,$5; \
    }' > $scaledBed