# HG changeset patch # User devteam # Date 1506797827 14400 # Node ID e0a8fba7ed2f045470f67f4a9df1da938c0c08e5 # Parent cfc3ad769dba052b3b834b16cc685e4fc93b6b74 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tool_collections/galaxy_sequence_utils/fastq_paired_end_interlacer commit f2582539542b33240234e8ea6093e25d0aee9b6a diff -r cfc3ad769dba -r e0a8fba7ed2f fastq_paired_end_interlacer.py --- a/fastq_paired_end_interlacer.py Thu Jan 23 12:31:16 2014 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,59 +0,0 @@ -#Florent Angly -import sys -from galaxy_utils.sequence.fastq import fastqReader, fastqWriter, fastqNamedReader, fastqJoiner - -def main(): - mate1_filename = sys.argv[1] - mate1_type = sys.argv[2] or 'sanger' - mate2_filename = sys.argv[3] - mate2_type = sys.argv[4] or 'sanger' - outfile_pairs = sys.argv[5] - outfile_singles = sys.argv[6] - - if mate1_type != mate2_type: - print "WARNING: You are trying to interlace files of two different types: %s and %s." % ( mate1_type, mate2_type ) - return - - type = mate1_type - joiner = fastqJoiner( type ) - out_pairs = fastqWriter( open( outfile_pairs, 'wb' ), format = type ) - out_singles = fastqWriter( open( outfile_singles, 'wb' ), format = type ) - - # Pairs + singles present in mate1 - nof_singles = 0 - nof_pairs = 0 - mate2_input = fastqNamedReader( open( mate2_filename, 'rb' ), format = type ) - i = None - for i, mate1 in enumerate( fastqReader( open( mate1_filename, 'rb' ), format = type ) ): - mate2 = mate2_input.get( joiner.get_paired_identifier( mate1 ) ) - if mate2: - out_pairs.write( mate1 ) - out_pairs.write( mate2 ) - nof_pairs += 1 - else: - out_singles.write( mate1 ) - nof_singles += 1 - - # Singles present in mate2 - mate1_input = fastqNamedReader( open( mate1_filename, 'rb' ), format = type ) - j = None - for j, mate2 in enumerate( fastqReader( open( mate2_filename, 'rb' ), format = type ) ): - mate1 = mate1_input.get( joiner.get_paired_identifier( mate2 ) ) - if not mate1: - out_singles.write( mate2 ) - nof_singles += 1 - - if (i is None) and (j is None): - print "Your input files contained no valid FASTQ sequences." - else: - print 'There were %s single reads.' % ( nof_singles ) - print 'Interlaced %s pairs of sequences.' % ( nof_pairs ) - - mate1_input.close() - mate2_input.close() - out_pairs.close() - out_singles.close() - - -if __name__ == "__main__": - main() diff -r cfc3ad769dba -r e0a8fba7ed2f fastq_paired_end_interlacer.xml --- a/fastq_paired_end_interlacer.xml Thu Jan 23 12:31:16 2014 -0500 +++ b/fastq_paired_end_interlacer.xml Sat Sep 30 14:57:07 2017 -0400 @@ -1,35 +1,78 @@ - - on paired end reads - - galaxy_sequence_utils - - fastq_paired_end_interlacer.py '$input1_file' '${input1_file.extension[len( 'fastq' ):]}' '$input2_file' '${input2_file.extension[len( 'fastq' ):]}' '$outfile_pairs' '$outfile_singles' - - - - - - - - - - - - - - - - - - - - - - - - - - + + on paired end reads + + galaxy_sequence_utils + + + + + + + + + + + + + + + + + + + + + + + reads['reads_selector'] == 'paired' + + + reads['reads_selector'] == 'paired' + + + reads['reads_selector'] == 'paired_collection' + + + reads['reads_selector'] == 'paired_collection' + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ]]> + + 10.1093/bioinformatics/btq281 + diff -r cfc3ad769dba -r e0a8fba7ed2f tool_dependencies.xml --- a/tool_dependencies.xml Thu Jan 23 12:31:16 2014 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,6 +0,0 @@ - - - - - -