bedtools genomecov -bg -strand + -ibam <file> -g <genomeindex> -scale <scalingfactor>   | awk '{print $1"\t"$2"\t"$3"\t-"$4}'     >    <name>.minus.bedGraph
bedtools genomecov -bg -strand  - -ibam <file> -g <genomeindex> -scale <scalingfactor> >  <name>.plus.bedGraph

bedGraphToBigWig <name>.plus.bedGraph    <genomeindex> ${file1%%.bam}.plus.bigWig
bedGraphToBigWig <name>.minus.bedGraph <genomeindex> ${file1%%.bam}.minus.bigWig 


thoc5.het == 1000000/1064914 =0.939043
scaling.brennecke.mutant == 1000000/1103542 = 0.906173