Mercurial > repos > iuc > circos

annotate fasta-to-karyotype.py @ 4:6d45ebe16927 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit f5db2969e8602eff32ec799d237875b911578fe5"
author iuc
date Thu, 09 Jan 2020 15:02:34 +0000
parents 890ef899a3d7
children b84d385679e2
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
8464186c8503 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit a41493893bdcbe330434db9c5851719012b62fa8
iuc
parents:
diff changeset
1 #!/usr/bin/env python
8464186c8503 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit a41493893bdcbe330434db9c5851719012b62fa8
iuc
parents:
diff changeset
2 import sys
8464186c8503 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit a41493893bdcbe330434db9c5851719012b62fa8
iuc
parents:
diff changeset
3
8464186c8503 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit a41493893bdcbe330434db9c5851719012b62fa8
iuc
parents:
diff changeset
4 from Bio import SeqIO
8464186c8503 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit a41493893bdcbe330434db9c5851719012b62fa8
iuc
parents:
diff changeset
5
2
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
6 # Process fasta data, extracting only headers
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
7 for idx, seq in enumerate(SeqIO.parse(sys.argv[1], "fasta")):
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
8 sys.stdout.write(
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
9 "chr - {seq_id} {seq_id} 0 {length} {theme}-{color}\n".format(
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
10 seq_id=seq.id,
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
11 idx=idx,
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
12 length=len(seq),
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
13 color=((idx + 1) % 12),
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
14 theme=sys.argv[2],
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
15 )
890ef899a3d7 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/circos commit 076837a2e9c2b6ececcea4aa286ea7a412387a96"
iuc
parents: 0
diff changeset
16 )