This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Take A Subset Of A Fastq Paired-End Sample

Hi,

I have two paired-end fastq compressed files coming from HiSeq RNA-SEq experiment, ie., pair.1.fastq.gz and pair.2.fastq.gz.

The files are very large so I wanted to just take a few million/thousand reads from each of them (by their respective pairs) and use that file for trying/debuuging purposes.

The results should be two paired-end files, i.e., pair.test.1.fastq.gz and pair.test.2.fastq.gz.

I'd be happy to hear some suggestions on how to do this or hear about tools available, thanks!

paired-end fastq rna-seq illumina
• 15,274 views • link
updated by Ram Ram • written by dfernan dfernan

thanks Pierre, I didn't realize someone else asked about it!

• 0 views • link
written by dfernan dfernan

2 answers

Hi,

Assuming that the reads are in same order in both of the files. I would do like this:

zcat pair.1.fastq.gz | sed -n 1,4000000p > pair_1_millions.fastq
zcat pair.2.fastq.gz | sed -n 1,4000000p > pair_2_millions.fastq

Thanks,
Rahul

• 280 views • link
updated by Ram Ram • written by Rahul Sharma Rahul Sharma

Hi, thanks a lot, however, I am not sure if the reads are in the same order, I'd like to add that I am pairing them correctly...

• 0 views • link
updated by Ram Ram • written by dfernan dfernan

I believe that a better option for paired-end data is to use fastq-sample from fastq-tools:

fastq-sample -n 5000000 pair_R1.fastq.gz pair_R2.fastq.gz -o pair_5M_R
• 0 views • link
written by sebastian.gregoricchio • sebastian.gregoricchio

Log in to answer this question.

More posts like this