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

Repository seq_filter_by_mapping
Owner: peterjc
Synopsis: Filter sequencing reads using SAM/BAM mapping files
This tool is a short Python script (using Biopython library functions) which divides
a FASTA, FASTQ, or SFF file in two, those sequences which do or do not map according
to given SAM/BAM file(s).

Example uses include mapping of FASTQ reads against a known contaminant in order
to remove reads prior to a de novo assembly.

See https://github.com/peterjc/pico_galaxy/tree/master/tools/seq_filter_by_mapping
Type: unrestricted
Revision: 6:566e456371c6
This revision can be installed: True
Times cloned / installed: 592

Repository README files - may contain important installation or license information

Galaxy tool to filter FASTA, FASTQ or SFF sequences by SAM/BAM mapping

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

This tool is a short Python script (using Biopython library functions) which divides a FASTA, FASTQ, or SFF file in two, those sequences which do or do not map according to given SAM/BAM file(s).

Example uses include mapping of FASTQ reads against a known contaminant in order to remove reads prior to a de novo assembly.

This tool is available from the Galaxy Tool Shed at:

See also related tools:

Automated Installation

This should be straightforward using the Galaxy Tool Shed, which should be able to automatically install the dependency on Biopython and samtools 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:

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

The suggested location is a dedicated tools/seq_filter_by_mapping/ 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="seq_filter_by_mapping/seq_filter_by_mapping.xml" />

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 seq_filter_by_mapping

You will also need to install Biopython 1.54 or later. That's it.

History

Version Changes
v0.0.1
  • Initial version.
v0.0.2
  • Fixed some error messages.
v0.0.3
  • Report counts for FASTQ as done for FASTA and SFF files.
v0.0.4
  • Use the format_source=... tag.
  • Reorder XML elements (internal change only).
  • Planemo for Tool Shed upload (.shed.yml, internal change only).
v0.0.5
  • Python script cleanups (internal change only).
  • Depends on Biopython 1.67 via legacy Tool Shed package or bioconda.
  • Use <command detect_errors="aggressive"> (internal change only).
  • Single quote command line arguments (internal change only).
v0.0.6
  • Python 3 compatible print function.
v0.0.7
  • Script works on Python 2 and 3 (fixed input file mode)
v0.0.8
  • Bumped Biopython dependency version for Python 3 fixes.

Developers

Development is on this GitHub repository: https://github.com/peterjc/pico_galaxy/tree/master/tools/seq_filter_by_mapping

Much of the code was copied from my older tool: https://github.com/peterjc/pico_galaxy/tree/master/tools/seq_filter_by_id

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/seq_filter_by_mapping/
...

or:

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

To just build and check the tar ball, use:

$ planemo shed_upload --tar_only tools/seq_filter_by_mapping/
...
$ tar -tzf shed_upload.tar.gz
test-data/SRR639755_mito_pairs.fastq.gz
test-data/SRR639755_sample_by_coord.sam
test-data/SRR639755_sample_lax.fastq
test-data/SRR639755_sample_strict.fastq
tools/seq_filter_by_mapping/README.rst
tools/seq_filter_by_mapping/seq_filter_by_mapping.py
tools/seq_filter_by_mapping/seq_filter_by_mapping.xml
tools/seq_filter_by_mapping/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_samtools_0_1_19 revision c9bd782f5342 owned by iuc

Name Version Type
samtools 0.1.19 package

Contents of this repository

Name Description Version Minimum Galaxy Version
from SAM/BAM file 0.0.8 16.01

Categories
Assembly - Tools for working with assemblies
Fasta Manipulation - Tools for manipulating fasta data
Fastq Manipulation - Tools for manipulating fastq data
Next Gen Mappers - Tools for the analysis and handling of Next Gen sequencing data
SAM - Tools for manipulating alignments in the SAM format
Sequence Analysis - Tools for performing Protein and DNA/RNA analysis