I was trying to use seqtk to randomly pickup 32,000,000 sequences from my sequence pool.
I first combined all of sequences from two samples with :
zcat sub592_1.headcropped.fastq.gz sub592_2.headcropped.fastq.gz > sub592.fastq.gz
then verified that the total number of reads in my file by :
`cat sub592.fastq.gz | echo $((`wc -l`/4))` #indicate I have 39,965,078
but when I am using seqtk to draw 32million reads from the pool, it did not do the job for me, instead, it just said “killed” at the end, could you help me point to a direction where could it go wrong?
seqtk sample ~/Desktop/rnaseq/trimmed/sub592.fastq.gz 3.2e+7 > s592_32.fastq
Thank you very much for your help and have a great day !