seqtk produces different number of reads for paired end files

lh3 / seqtk

Toolkit for processing sequences in FASTA/Q formats

MIT License

1.38k stars 308 forks source link

seqtk produces different number of reads for paired end files #198

Closed ruy-jauregui closed 1 year ago

ruy-jauregui commented 1 year ago

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?

lh3 commented 1 year ago

Probably the input fastq are not paired.