Move callable region calculation to use bedtools with groupby thanks to tips from @arq5x. Provides improved speed and memory usage over pysam based implementation. Fixes for pybedtools temp directory specification