Identify parallelizable islands of mapped reads surrounded by unmapped regions. Add pybedtools as dependency for identifying blocks