view fasta-to-karyotype.py @ 2:014a21767ac4 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
author iuc
date Tue, 17 Sep 2019 16:54:57 -0400
parents ef5f8bbf7730
children 740057a5126d
line wrap: on
line source

#!/usr/bin/env python
import sys

from Bio import SeqIO

# Process fasta data, extracting only headers
for idx, seq in enumerate(SeqIO.parse(sys.argv[1], "fasta")):
    sys.stdout.write(
        "chr - {seq_id} {seq_id} 0 {length} {theme}-{color}\n".format(
            seq_id=seq.id,
            idx=idx,
            length=len(seq),
            color=((idx + 1) % 12),
            theme=sys.argv[2],
        )
    )