workflows/pangenome-generate/relabel-seqs.py


1
2
3
4
5
6
7
8
9
10
11
12
13

import sys

reads = $(inputs.readsFA)
subjects = $(inputs.subjects)

for i, r in enumerate(reads):
    with open(r["path"], "rt") as fa:
        fa.readline()
        print(">"+subjects[i])
        data = fa.read(8096)
        while data:
            sys.stdout.write(data)
            data = fa.read(8096)