total_reads=1882597
target_reads=10000000

scaling_factor=$(echo "scale=6; $target_reads / $total_reads" | bc)

zcat right_ventricle.Unique.stranded.bg.gz | \
awk -v scale_factor="$scaling_factor" 'BEGIN {OFS="\t"} { print $1, $2, $3, $4 * scale_factor }' | \
gzip > right_ventricle.Unique.stranded_10M.bg.gz
