# HG changeset patch # User peterjc # Date 1494436192 14400 # Node ID 530c8d6fedd890bdfd9d599f232397b6c54055dd # Parent 9b074c1db68e321ab34158a60104b11c0f795e28 v0.0.15 - internal changes diff -r 9b074c1db68e -r 530c8d6fedd8 tools/seq_primer_clip/README.rst --- a/tools/seq_primer_clip/README.rst Thu Feb 02 11:52:37 2017 -0500 +++ b/tools/seq_primer_clip/README.rst Wed May 10 13:09:52 2017 -0400 @@ -69,8 +69,10 @@ - Reorder XML elements (internal change only). - Planemo for Tool Shed upload (``.shed.yml``, internal change only). - Fixed input file help text. -v0.0.14 - Updated to point at Biopython 1.67 (latest version in Tool Shed). +v0.0.14 - Depends on Biopython 1.67 via legacy Tool Shed package or bioconda. - Explicit dependency on ``galaxy_sequence_utils``. +v0.0.15 - Use ```` (internal change only). + - Single quote command line arguments (internal change only) ======= ====================================================================== @@ -89,17 +91,17 @@ Planemo commands (which requires you have set your Tool Shed access details in ``~/.planemo.yml`` and that you have access rights on the Tool Shed):: - $ planemo shed_update -t testtoolshed --check_diff ~/repositories/pico_galaxy/tools/seq_primer_clip/ + $ planemo shed_update -t testtoolshed --check_diff tools/seq_primer_clip/ ... or:: - $ planemo shed_update -t toolshed --check_diff ~/repositories/pico_galaxy/tools/seq_primer_clip/ + $ planemo shed_update -t toolshed --check_diff tools/seq_primer_clip/ ... To just build and check the tar ball, use:: - $ planemo shed_upload --tar_only ~/repositories/pico_galaxy/tools/seq_primer_clip/ + $ planemo shed_upload --tar_only tools/seq_primer_clip/ ... $ tar -tzf shed_upload.tar.gz test-data/MID4_GLZRM4E04_rnd30.fasta diff -r 9b074c1db68e -r 530c8d6fedd8 tools/seq_primer_clip/seq_primer_clip.py --- a/tools/seq_primer_clip/seq_primer_clip.py Thu Feb 02 11:52:37 2017 -0500 +++ b/tools/seq_primer_clip/seq_primer_clip.py Wed May 10 13:09:52 2017 -0400 @@ -29,8 +29,10 @@ NOTE: Currently it uses Python's regular expression engine for finding the primers, which for my needs is fast enough. """ + +import re import sys -import re + from galaxy_utils.sequence.fasta import fastaReader, fastaWriter from galaxy_utils.sequence.fastq import fastqReader, fastqWriter @@ -147,8 +149,8 @@ # We'll use a set to remove any duplicate patterns # if letter not in "NX": pattern = seq[:i] + "N" + seq[i + 1:] - assert len(pattern) == len(seq), "Len %s is %i, len %s is %i" \ - % (pattern, len(pattern), seq, len(seq)) + assert len(pattern) == len(seq), ("Len %s is %i, len %s is %i" + % (pattern, len(pattern), seq, len(seq))) yield make_reg_ex(pattern) if mm >= 2: for i, letter in enumerate(seq): @@ -158,8 +160,8 @@ # We'll use a set to remove any duplicate patterns # if letter not in "NX": pattern = seq[:i] + "N" + seq[i + 1:i + 1 + k] + "N" + seq[i + k + 2:] - assert len(pattern) == len(seq), "Len %s is %i, len %s is %i" \ - % (pattern, len(pattern), seq, len(seq)) + assert len(pattern) == len(seq), ("Len %s is %i, len %s is %i" + % (pattern, len(pattern), seq, len(seq))) yield make_reg_ex(pattern) diff -r 9b074c1db68e -r 530c8d6fedd8 tools/seq_primer_clip/seq_primer_clip.xml --- a/tools/seq_primer_clip/seq_primer_clip.xml Thu Feb 02 11:52:37 2017 -0500 +++ b/tools/seq_primer_clip/seq_primer_clip.xml Wed May 10 13:09:52 2017 -0400 @@ -1,18 +1,14 @@ - + Trim off 5' or 3' primers galaxy_sequence_utils biopython - Bio - - - - - - seq_primer_clip.py --version - -seq_primer_clip.py $input_file $input_file.ext $primer_fasta $primer_type $mm $min_len $keep_negatives $output_file + +python $__tool_directory__/seq_primer_clip.py --version + + +python $__tool_directory__/seq_primer_clip.py $input_file $input_file.ext $primer_fasta $primer_type $mm $min_len $keep_negatives $output_file