def count_elements_in_all_bounding_regions(track_name, genome='testgenome', allow_overlaps=False): bounding_regions = BoundingRegionShelve(genome, track_name, allow_overlaps).getAllBoundingRegions() track = PlainTrack(track_name) num_elements = 0 for tv in [track.getTrackView(region) for region in bounding_regions]: num_elements += count_elements(tv) return num_elements
def get_track_view(track_name, genome_region): track = PlainTrack(track_name) return track.getTrackView(genome_region)