[Seqan-dev] Reading FastQ files from GZ archive

[ngs geek <ngsgeek@gmail.com>]

Hi All:

I'm trying to filter reads contained in a gz file based on quality score; I was able to stream the id and the 

sequence with GZFile stream and using a singlepass reader (I don't want to keep the reads in memory).

Unfortunately, I can't access the quality string.

 RecordReader<Stream<GZFile>, SinglePass<> > reader(fastq_file);

  CharString id;

  String<Dna5Q> seq;

CharString qual;

while (!atEnd(reader))

  {

      if (readRecord(id, seq, reader, Fastq()) != 0)

      {

          std::cerr << "Problem with your FASTQ file." << std::endl;

          return 1;

      }

      // reading qual

  }

Any idea on this? I tried assignQualities etc. with no luck.

Thanks

JD