seqan14.cc

#include <seqan/seq_io.h>

int main(int argc, char **argv)
{
    std::string id, seq, qual;

    if (argc < 2) {
        return 1;
    }

    seqan::SequenceStream seqFileIn;
    seqan::open(seqFileIn, argv[1]);
    size_t seqs = 0;

    while (!seqan::atEnd(seqFileIn)) {
        seqan::readRecord(id, seq, qual, seqFileIn);
        std::cout << seqs++ << ' ' <<  id << '\t' << seq << '\t' << qual << '\n';
    }

    return 0;
}