view fasta-to-karyotype.py @ 2:890ef899a3d7 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:24 -0400
parents 8464186c8503
children b84d385679e2
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],
        )
    )