view split.py @ 3:e0734e88a104 draft default tip

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/vardict commit 80e3742f499e8efabad1fe7627201466d0bdd190"
author iuc
date Thu, 27 Jan 2022 15:19:23 +0000
parents
children
line wrap: on
line source

import sys


fai = sys.argv[1]
chunk_size = int(sys.argv[2])
overlap = int(sys.argv[3])  # Base pairs
with open(fai, 'r') as infile:
    for line in infile:
        name = line.split('\t')[0]
        stop = int(line.split('\t')[1])
        start = 1
        while start < stop:
            start = max(1, start - overlap)
            print('\t'.join([name, str(start),
                             str(min(start + chunk_size, stop))]))
            start += chunk_size