Merge multiple fq read files

[hosseinv]

I have multiple Illumina Hi-seq 2000 fastq.gz files for each individual as follows;

sample1 lane1 read 1_001
sample1 lane1 read 1_002
sample1 lane1 read 2_001
sample1 lane1 read 2_002

sample1 lane2 read 1_001
sample1 lane2 read 1_002
sample1 lane2 read 2_001
sample1 lane2 read 2_002

sample1 lane3 read 1_001
sample1 lane3 read 1_002
sample1 lane3 read 2_001
sample1 lane3 read 2_002

they are all for a single individual. what script in Linux console would be best to merge all files in a final merged file? the idea is to do analysis in Galaxy.

THank you all.

[hosseinv]

They are paired-end reads.

[maubp]

Galaxy can concatenate files together, and might actually be easier to upload each file one by one to Galaxy (gzip them first), since uploading a single large file is harder.

If you do want to concatenate files at the command line, you can use the command 'cat', as in:

cat fileA fileB fileC > combinedfileABC

[hosseinv]

Thanks Peter for your answer
the thing that I want to know is the order of putting files in my cat command. if you have a look on my example files you'll see I've got sequences of one individual (sample 1) in 3 lanes and different number of reads both for read 1s (001 and 002) and read 2s (001 nad 002) in each lane.
plus I know how to merge two files in Galaxy, but don't know how to merge multiple files.

Thanks again.

[DFJ111]

Not sure why you have duplicates of sample1 lane1 read 1_002, sample1 lane2 read 1_002, and sample1 lane3 read 1_002. Typo from the sequencing lab? If all you want to do is QC (via FastQC for example) I don't think it matters what order they are in. Also, if you do want to do FastQC you should assess each lane separately. If you want to do something else to them, the order or merging (or whether to merge them at all) depends on what that "something else" is.

Also, you can merge multiple files in GALAXY, use "concatenate head-to-tail". I just use cat, it's quicker.

[hosseinv]

Thanks DFJ111

It's not a typo. The thing that I want to do is to map the sequences against the reference genes and find the polymorphism.

[maubp]

Which tool are you planning to use for the mapping, and does it require the paired reads in any specific order (e.g. interleaved in one file) or as separate files (forward and reverse reads)?

[hosseinv]

Hi maubp

I will be using BWA for mapping. and I think I'll hav to treat each of the reads individually and at some stage I can pool the BAM or SAM files together.

Thanks for your attention .