This document is for an in-development version of Galaxy. You can alternatively view this page in the latest release if it exists or view the top of the latest release's documentation.
Source code for galaxy.datatypes.converters.fastqsolexa_to_fasta_converter
#!/usr/bin/env python """ convert fastqsolexa file to separated sequence and quality files. assume each sequence and quality score are contained in one line the order should be: 1st line: @title_of_seq 2nd line: nucleotides 3rd line: +title_of_qualityscore (might be skipped) 4th line: quality scores (in three forms: a. digits, b. ASCII codes, the first char as the coding base, c. ASCII codes without the first char.) Usage: %python fastqsolexa_to_fasta_converter.py <your_fastqsolexa_filename> <output_seq_filename> <output_score_filename> """ import sys assert sys.version_info[:2] >= (2, 4)[docs]def stop_err(msg): sys.stderr.write("%s" % msg) sys.exit()def __main__(): infile_name = sys.argv fastq_block_lines = 0 seq_title_startswith = '' with open(infile_name) as fh, open(sys.argv, 'w') as outfile: for i, line in enumerate(fh): line = line.rstrip() # eliminate trailing space and new line characters if not line or line.startswith('#'): continue fastq_block_lines = (fastq_block_lines + 1) % 4 line_startswith = line[0:1] if fastq_block_lines == 1: # line 1 is sequence title if not seq_title_startswith: seq_title_startswith = line_startswith if seq_title_startswith != line_startswith: stop_err('Invalid fastqsolexa format at line %d: %s.' % (i + 1, line)) outfile.write('>%s\n' % line[1:]) elif fastq_block_lines == 2: # line 2 is nucleotides outfile.write('%s\n' % line) else: pass if __name__ == "__main__": __main__()