view fasta-to-karyotype.py @ 3:fd59d783248a draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit f04fe62fb7dbad75e597da8cd674cfe36ed3d80a"

#!/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],
        )
    )