help='Show this help message and exit.')
(options, args) = parser.parse_args(argv)
if len(args) < 2:
parser.error("genome_dir is required")
return options, args
################################################################################
if __name__ == '__main__':
options, args = process_command_line(sys.argv)
genome_dir = args[1]
assert os.path.exists(genome_dir)
outhandle = sys.stdout
if options.outfile:
outhandle = open(options.outfile, "w")
scanner = AstretchScan(genome_dir, options.f_ival)
i = 0
for chrom, start, end, strand in scanner.scan():
i += 1
outhandle.write("\t".join(map(str, [
chrom, start, end, i, 0, strand]))+"\n")
if options.outfile:
outhandle.close()
