process_qseq.py – Given a directory of per-swath qseq files, this script generates a single fastq per lane.

Warning! A bug has been found in process_qseq.py in QIIME 1.3.0 which results in many of the sequences not being written to the fastq. A patch will be made available shortly. In the meantime, e-mail gregcaporaso@gmail.com for a patch if you need to use this script now. The issue has been corrected in the svn version of QIIME. The bug is described here.

Description:

Usage: process_qseq.py [options]

Input Arguments:

Note

[REQUIRED]

-i, --input_dir

The input directory

-o, --output_dir

The output directory

-r, --read

The read number to consider

[OPTIONAL]

-l, --lanes

The lane numbers to consider, comma-separated [defaut: 1,2,3,4,5,6,7,8]

-b, --bases

The number of bases to include (useful for slicing a barcode) [defaut: all]

Output:

Generate fastq files from all lanes of read 1 data in the current directory.

process_qseq.py -i ./ -o ./fastq/ -r 1

Generate fastq files from all lanes of read 2 data in the current directory, truncating the sequences after the first 12 bases.

process_qseq.py -i ./ -o ./fastq/ -r 2 -b 12