Mercurial > repos > peterjc > fastq_filter_by_id
diff tools/fastq/fastq_filter_by_id.py @ 2:d570cc324779
Migrated tool version 0.0.4 from old tool shed archive to new tool shed repository
author | peterjc |
---|---|
date | Tue, 07 Jun 2011 17:24:08 -0400 |
parents | 10e963c79a45 |
children |
line wrap: on
line diff
--- a/tools/fastq/fastq_filter_by_id.py Tue Jun 07 17:23:49 2011 -0400 +++ b/tools/fastq/fastq_filter_by_id.py Tue Jun 07 17:24:08 2011 -0400 @@ -1,6 +1,9 @@ #!/usr/bin/env python """Filter a FASTQ file with IDs from a tabular file, e.g. from BLAST. +NOTE - This script is now OBSOLETE, having been replaced by a new verion +which handles FASTA, FASTQ and SFF all in one. + Takes five command line options, tabular filename, ID column numbers (comma separated list using one based counting), input FASTA filename, and two output FASTA filenames (for records with and without the given IDs). @@ -13,10 +16,10 @@ in column one, and the ID of the match from the database is in column two. Here sensible values for the column numbers would therefore be "1" or "2". -This script is copyright 2010 by Peter Cock, SCRI, UK. All rights reserved. +This script is copyright 2010-2011 by Peter Cock, SCRI, UK. All rights reserved. See accompanying text file for licence details (MIT/BSD style). -This is version 0.0.2 of the script. +This is version 0.0.4 of the script. """ import sys from galaxy_utils.sequence.fastq import fastqReader, fastqWriter @@ -86,7 +89,6 @@ #The [1:] is because the fastaReader leaves the @ on the identifer. if not record.identifier or record.identifier.split()[0][1:] not in ids: negative_writer.write(record) - positive_writer.close() negative_writer.close() else: stop_err("Neither output file requested")