seq_reader.cpp

An example of reading a gzipped fastq file.

#include "btllib/seq_reader.hpp"

#include <iostream>

int
main()
{
  for (const auto& record : btllib::SeqReader(
         "my_reads.fq.gz", btllib::SeqReader::Flag::SHORT_MODE)) {
    std::cout << record.seq << '\n' << record.qual << '\n';
  }

  /*
   * In Python, this would be:
   *
   * with btllib.SeqReader("my_reads.fq.gz", btllib.SeqReaderFlag.SHORT_MODE) as
   reader: for record in reader: print(record.seq, record.qual)
   *
   */

  return 0;
}