HTSeq—a Python framework to work with high-throughput sequencing data
bedtools : a flexible suite of utilities for comparing genomic features .