Main body (bottom of file, currently line 143), need to add --sci to the bedmap command. Without it, FDR values below 1e-6 are erroneously being passed through as 0.
annotate() (near bottom of file, currently line 131), need to change 10int(col9 + 0.5) to int(10col9 + 0.5)
Simple fixes needed in hsmerge.sh:
Main body (bottom of file, currently line 143), need to add --sci to the bedmap command. Without it, FDR values below 1e-6 are erroneously being passed through as 0.
annotate() (near bottom of file, currently line 131), need to change 10int(col9 + 0.5) to int(10col9 + 0.5)