#!/bin/sh
[ $# != 1 ] && echo "usage: mkRandomGoldLift db " >&2 && exit 1;
db=$1
grep random /cluster/data/$db/chrom.sizes | while read R
do
        chr=`echo $R | awk '{print $1}'`
        size=`echo $R | awk '{print $2}'`
        hgsql -N -e "select frag, fragEnd, chrom, chromStart, chromEnd from gold where chrom=\"$chr\";" $db | \
awk '
BEGIN {size="'$size'"}
{
        printf "%s\t%s\t%s\t%s\t%s\n", $4, $1, $2, $3, size
}
'
done
