How to split fastq into small fastq based on barcode?

[peterrjp]

Dear All,
I have a big fastq file which include 34 samples. I want to split it into 34 small fastq files based on barcode sequences.
I tried to do it with the script "split_libraries_fastq" in Qiime. However, a barcode read fastq file should be used in this script. I don't know how to get this barcode read fastq file, so I can't use this script to solve the problem.
Is there any other method to solve the problem? Thank you!

Peter

[wlchew]

FastX-toolkit?

http://hannonlab.cshl.edu/fastx_toolkit/

[Vinz]

Is this MiSeq or HiSeq data?
Are you sure, the index was sequenced? If you do not specify it in the sample sheet, no index is sequenced. Then there is no way to determine the samples and split the fastq afterwards.
Open the fastq with an appropriate text viewer and see if barcode information is present in the header of each read.

[GenoMax]

Two things you need to clarify. Is the barcode "inline" (part of the sequence read) or was it sequenced separately (as in Illumina multiplexing)?

As Vinz said above if you are not sure then you can "cat" or "zcat" the sequence file (pipe through "more") and then post a few sequences here so someone can help.

[GW_OK]

As Vinz said above if you are not sure then you can "cat" or "zcat" the sequence file (pipe through "more") and then post a few sequences here so someone can help.

Oh, man, use 'less' or 'zless' (press 'q' to return to command line). Way better than 'cat' piped through 'more'.

[acoada]

try perl.
you can store barcode list in a hash.

[peterrjp]

Dear All,
Thank you for your replies! I have solve the problem with the FASTX-ToolKit software. Besides, I also found the way to get a barcode.fastq with the script split_libraries_fastq.py (http://qiime.org/tutorials/extractin...astq_data.html).

Peter