# HG changeset patch # User peterjc # Date 1486657104 18000 # Node ID eef06667d973eb951d85c72789073410e1110252 # Parent 839cb62120a47e565644dc5ad0edcdb7e5257bc0 planemo upload for repository https://github.com/peterjc/pico_galaxy/tree/master/tools/align_back_trans commit 5fb6ee31f33e8e24c1ae5d56c93ac3e0157dad72-dirty diff -r 839cb62120a4 -r eef06667d973 tools/align_back_trans/README.rst --- a/tools/align_back_trans/README.rst Wed Feb 01 06:52:08 2017 -0500 +++ b/tools/align_back_trans/README.rst Thu Feb 09 11:18:24 2017 -0500 @@ -70,6 +70,7 @@ - Use ``format_source=...`` tag. - Planemo for Tool Shed upload (``.shed.yml``, internal change only). v0.0.7 - Minor Python code style improvements (internal change only). +v0.0.8 - Python coding style change to avoid lamba (internal change only). ======= ====================================================================== diff -r 839cb62120a4 -r eef06667d973 tools/align_back_trans/align_back_trans.py --- a/tools/align_back_trans/align_back_trans.py Wed Feb 01 06:52:08 2017 -0500 +++ b/tools/align_back_trans/align_back_trans.py Thu Feb 09 11:18:24 2017 -0500 @@ -18,12 +18,14 @@ """ import sys -from Bio.Seq import Seq -from Bio.Alphabet import generic_protein + +from Bio import AlignIO +from Bio import SeqIO + from Bio.Align import MultipleSeqAlignment -from Bio import SeqIO -from Bio import AlignIO +from Bio.Alphabet import generic_protein from Bio.Data.CodonTable import ambiguous_generic_by_id +from Bio.Seq import Seq if "-v" in sys.argv or "--version" in sys.argv: print "v0.0.7" @@ -130,19 +132,23 @@ def alignment_back_translate(protein_alignment, nucleotide_records, key_function=None, gap=None, table=0): """Thread nucleotide sequences onto a protein alignment.""" # TODO - Separate arguments for protein and nucleotide gap characters? - if key_function is None: - key_function = lambda x: x if gap is None: gap = "-" aligned = [] - for protein in protein_alignment: - try: - nucleotide = nucleotide_records[key_function(protein.id)] - except KeyError: - raise ValueError("Could not find nucleotide sequence for protein %r" - % protein.id) - aligned.append(sequence_back_translate(protein, nucleotide, gap, table)) + try: + if key_function is None: + for protein in protein_alignment: + nucleotide = nucleotide_records[protein.id] + aligned.append(sequence_back_translate(protein, nucleotide, gap, table)) + else: + for protein in protein_alignment: + nucleotide = nucleotide_records[key_function(protein.id)] + aligned.append(sequence_back_translate(protein, nucleotide, gap, table)) + except KeyError: + raise ValueError("Could not find nucleotide sequence for protein %r" + % protein.id) + return MultipleSeqAlignment(aligned) diff -r 839cb62120a4 -r eef06667d973 tools/align_back_trans/align_back_trans.xml --- a/tools/align_back_trans/align_back_trans.xml Wed Feb 01 06:52:08 2017 -0500 +++ b/tools/align_back_trans/align_back_trans.xml Thu Feb 09 11:18:24 2017 -0500 @@ -1,4 +1,4 @@ - + Gives a codon aware alignment biopython