lh3 / seqtk

Toolkit for processing sequences in FASTA/Q formats

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

seqtk produces different number of reads for paired end files

ruy-jauregui opened this issue · comments

Hello!
I run these commands:
seqtk sample -s100 cJ288_1.fq 90000 > ./subsample/cJ288_1.fq
seqtk sample -s100 cJ288_2.fq 90000 > ./subsample/cJ288_2.fq

When I count the reads in the output I get:
cJ288_1.fq:89995
cJ288_2.fq:89998

What's going on?

Probably the input fastq are not paired.