Illumina Paired End FASTQ

[kjsalimian]

Hi guys,

I'm very new to this. I have a paired end data set from HiSeq1000. I want to take the first 10,000 or 100,000 reads out of ~40mil reads to use for tests rather than putting the entire 40 mil reads through the tests.

What is the easiest way to generate files of only the first 100,000 reads??

Thanks.

[swbarnes2]

I'm not really familar with HiSeq data foramts, but there are probably some kind of coordinates, maybe by tile. If the data below from the checked answer is correct, you could use a grep to get only the entries with '[sequencer name]_[run]:[lane]:10:'. That's be a fairly random.

Or simpler, head -400000 to get the top 1000k reads?

The ones in the very beginning are probably at the edge of the flow-cell, and will have more bad quality reads. Pulling from the middle might get you more good reads.


http://biostar.stackexchange.com/que...d-naming-tiles

[dpryan]

Have a look at this thread. You'll often find that the first batch of reads in a file are crap (at least from some machines), so you're better off randomly selecting them.