I need to randomly subsample a fastq file. I'm trying to use the shortbread package.
This are the commands I tried:
xx=FastqSampler("~/Desktop/Silano2014_RawReads/1_S1_R1.fastq", 100) xx class: FastqSampler file: 1_S1_R1.fastq status: n=100 current=0 added=0 total=0 buffer=0 ordered: FALSE set.seed(123L); p1 = yield(xx) class: ShortReadQ length: 100 reads; width: 151 cycles
Then I tried to write this new file:
But the file it created doesn't look like a fastq
head subset.fastq ?]K??<????ث?aٲhAm I doing something wrong?