File: seq_reader.cpp

package info (click to toggle)
btllib 1.7.5%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 4,336 kB
  • sloc: cpp: 79,742; python: 941; sh: 302; makefile: 5
file content (22 lines) | stat: -rw-r--r-- 458 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
#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;
}