Re: [Seqan-dev] Reading FastQ files from GZ archive

["Holtgrewe, Manuel" <manuel.holtgrewe@fu-berlin.de>]

Hi JD,

good catch, you found a bug I just fixed. The problem was that assignQualityValues() did not resize the target string, it does now. Update your SVN checkout and it should work.

Also, there now is a demo for reading FASTQ from gzipped files in core/demos/stream_read_fastq_gz.cpp

HTH

Manuel

---

From: ngs geek [ngsgeek@gmail.com]
Sent: Friday, March 30, 2012 7:44 PM
To: seqan-dev@lists.fu-berlin.de
Subject: [Seqan-dev] Reading FastQ files from GZ archive

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