logging.info("Skipping untargeted chromosomes %s",
" ".join(sorted(chroms_to_skip)))
skip_idx = accessible.chromosome.isin(chroms_to_skip)
accessible = accessible[~skip_idx]
else:
// Chromosome accessible sequence regions not known -- use heuristics
// (chromosome length is endpoint of last probe; skip initial
// <magic number> of bases that are probably telomeric)