view maf2bed.py @ 0:53c2be00bb6f draft default tip

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/jbrowse2 commit 0a86c88a95b0d1cc49d84544136de6556b95320f
author bgruening
date Wed, 05 Jun 2024 08:15:49 +0000
parents
children
line wrap: on
line source

#!/usr/bin/env python
# painfully converted from b0rken perl from
# https://unpkg.com/browse/jbrowse-plugin-mafviewer@1.0.6/dist/
# license is Apache2_license.txt included here

import sys

id = 0
buffer = ''
start = 0
end = 0
score = 0
chrom = ''

db = "%s." % sys.argv[1]
# Read input from stdin
for line in sys.stdin:
    line = line.strip()
    if not line:
        continue

    line = line.split()
    if line[0] == 's' and line[1].startswith(db):
        chrom = line[1]
        chrom = chrom.replace(db, '')
        start = int(line[2])
        end = int(line[2]) + int(line[3])
        line = line[1:]
        line = ':'.join(line)
        temp = line
        buffer = temp if buffer == '' else f"{buffer},{temp}"
    elif line[0] == 'a':
        score = int(line[1].split('=')[1])
        if id > 0:
            sys.stdout.write('\t'.join([chrom, '%d' % start, '%d' % end, f"{sys.argv[1]}_{id}", '%d' % score, buffer]) + '\n')
        id += 1
        buffer = ''
    elif line[0] == 's':
        line = line[1:]
        line = ':'.join(line)
        temp = line
        buffer = temp if buffer == '' else f"{buffer},{temp}"

sys.stdout.write('\t'.join([chrom, '%d' % start, '%d' % end, f"{sys.argv[1]}_{id}", '%d' % score, buffer]) + '\n')