Repository revision
repository tip
Select a revision to inspect and download versions of Galaxy utilities from this repository.

Repository sample_seqs
Name: sample_seqs
Owner: peterjc
Synopsis: Sub-sample sequences files (e.g. to reduce coverage)
Takes an input file of sequences (typically FASTA or FASTQ, but also
Standard Flowgram Format (SFF) is supported), and returns a new sequence
file sub-sampling from this (in the same format).
Type: unrestricted
Revision: 22:75b50b0a6250
This revision can be installed: True
Times cloned / installed: 496

Repository README files - may contain important installation or license information

Galaxy tool to sub-sample sequence files

This tool is copyright 2014-2017 by Peter Cock, The James Hutton Institute (formerly SCRI, Scottish Crop Research Institute), UK. All rights reserved. See the licence text below (MIT licence).

This tool is a short Python script (using Biopython library functions) to sub-sample sequence files (in a range of formats including FASTA, FASTQ, and SFF). This can be useful for preparing a small sample of data to test or time a new pipeline, or for reducing the read coverage in a de novo assembly.

This tool is available from the Galaxy Tool Shed at:

Automated Installation

This should be straightforward using the Galaxy Tool Shed, which should be able to automatically install the dependency on Biopython, and then install this tool and run its unit tests.

Manual Installation

There are just two files to install to use this tool from within Galaxy:

  • sample_seqs.py (the Python script)
  • sample_seqs.xml (the Galaxy tool definition)

The suggested location is in a dedicated tools/sample_seqs folder.

You will also need to modify the tools_conf.xml file to tell Galaxy to offer the tool. One suggested location is in the filters section. Simply add the line:

<tool file="sample_seqs/sample_seqs.xml" />

You will also need to install Biopython 1.62 or later.

If you wish to run the unit tests, also move/copy the test-data/ files under Galaxy's test-data/ folder. Then:

./run_tests.sh -id sample_seqs

That's it.

History

Version Changes
v0.0.1
  • Initial version.
v0.1.1
  • Using optparse to provide a proper Python command line API.
v0.1.2
  • Interleaved mode for working with paired records.
  • Tool definition now embeds citation information.
v0.2.0
  • Option to give number of sequences (or pairs) desired. This works by first counting all your sequences, then calculates the percentage required in order to sample them uniformly (evenly). This makes two passes through the input and is therefore slower.
v0.2.1
  • Was missing a file for the functional tests.
  • Included testing of stdout messages.
  • Includes testing of failure modes.
v0.2.2
  • Reorder XML elements (internal change only).
  • Use format_source=... tag.
  • Planemo for Tool Shed upload (.shed.yml, internal change only).
v0.2.3
  • Do the Biopython imports at the script start (internal change only).
  • Clarify paired read example in help text.
v0.2.4
  • Depends on Biopython 1.67 via legacy Tool Shed package or bioconda.
  • Style changes to Python code (internal change only).
v0.2.5
  • Use <command detect_errors="aggressive"> (internal change only).
  • Single quote command line arguments (internal change only).

Developers

This script and related tools are being developed on this GitHub repository: https://github.com/peterjc/pico_galaxy/tree/master/tools/sample_seqs

For pushing a release to the test or main "Galaxy Tool Shed", use the following 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 tools/sample_seqs/
...

or:

$ planemo shed_update -t toolshed --check_diff tools/sample_seqs/
...

To just build and check the tar ball, use:

$ planemo shed_upload --tar_only tools/sample_seqs/
...
$ tar -tzf shed_upload.tar.gz
test-data/MID4_GLZRM4E04_rnd30_frclip.pair_sample_N5.sff
test-data/MID4_GLZRM4E04_rnd30_frclip.sff
test-data/MID4_GLZRM4E04_rnd30_frclip.sample_C1.sff
test-data/MID4_GLZRM4E04_rnd30_frclip.sample_N5.sff
test-data/MID4_GLZRM4E04_rnd30_frclip.pair_sample_N5.sff
test-data/ecoli.fastq
test-data/ecoli.pair_sample_N100.fastq
test-data/ecoli.sample_C10.fastq
test-data/ecoli.sample_N100.fastq
test-data/get_orf_input.Suis_ORF.prot.fasta
test-data/get_orf_input.Suis_ORF.prot.pair_sample_C10.fasta
test-data/get_orf_input.Suis_ORF.prot.pair_sample_N100.fasta
test-data/get_orf_input.Suis_ORF.prot.sample_C10.fasta
test-data/get_orf_input.Suis_ORF.prot.sample_N100.fasta
tools/sample_seqs/README.rst
tools/sample_seqs/sample_seqs.py
tools/sample_seqs/sample_seqs.xml
tools/sample_seqs/tool_dependencies.xml

Licence (MIT)

Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

Dependencies of this repository

Repository package_numpy_1_9 revision 1e2b7e63dec2 owned by iuc (prior install required)
Repository package_atlas_3_10 revision 35d6b30ef137 owned by iuc (prior install required)
Repository package_python_2_7 revision a03693de8e0c owned by iuc (prior install required)
Repository package_zlib_1_2_8 revision 3348f1acee2a owned by iuc (prior install required)
Repository package_readline_6_3 revision de22495ddf2b owned by iuc (prior install required)
Repository package_bzlib_1_0 revision 29f79feefc77 owned by iuc (prior install required)
Repository package_xz_5_0_5 revision c9d30f1de10a owned by iuc (prior install required)
Repository package_sqlite_3_8_3 revision 9ba893c2353f owned by iuc (prior install required)
Repository package_gdbm_1_11 revision 77ea468f199e owned by iuc (prior install required)
Repository package_ncurses_6_0 revision 15e2190c87ca owned by iuc (prior install required)
Repository package_ncurses_6_0 revision 15e2190c87ca owned by iuc (prior install required)
Repository package_libtool_2_4 revision 5ccfb33b5622 owned by iuc (prior install required)

Name Version Type
biopython 1.67 package

Contents of this repository

Name Description Version Minimum Galaxy Version
e.g. to reduce coverage 0.2.5 16.01

Categories
Assembly - Tools for working with assemblies
Fasta Manipulation - Tools for manipulating fasta data
Fastq Manipulation - Tools for manipulating fastq data
Sequence Analysis - Tools for performing Protein and DNA/RNA analysis