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")