[Seqan-dev] Reading fastq records - missing records?

[Alexander Vowinkel <vowinkel.alexander@gmail.com>]

Hi Team,

here is another question:

Following code reads a fastq-file.

The file I'm testing on has 263333413 records:

$ zcat C563WACXX_1.fastq.gz | wc -l
1053333652

The problem is, that the code seems to only read

the first 4000000 records. I don't understand why.

Can you advise me here?

Thank you,

Alexander

PS: I removed try/catch for better readability

seqan::StringSet<seqan::CharString> ids;
seqan::StringSet<seqan::Dna5QString> seqs;
seqan::StringSet<seqan::CharString> quals;
unsigned batchSize = 10000;
do {
seqan::clear(ids); seqan::clear(seqs); seqan::clear(quals);
readRecords(ids, seqs, quals, inputFile, batchSize);
processRecords(ids, seqs, quals, barcodeGraph, vertexFilePtrMap);
totalProcessed += seqan::length(ids);
std::cout << "\r";
std::cout << "reads processed in total: " << totalProcessed;
std::cout.flush();
} while(seqan::length(ids) > 0);