#include #include #include using namespace seqan; int main(int argc, char const ** argv) { if (argc != 2) return 1; StringSet ids; StringSet seqs; SeqFileIn reader(argv[1]); while (!atEnd(reader)) { readRecords(ids, seqs, reader, 5000); } return 0; }