
I'm calling variants from exome sequencing data and I need to evaluate the efficiency of the capture and the coverage along the target regions. This sounds like a great use case for bedtools, your swiss-army knife for genomic arithmetic and interval manipulation.