FU Logo
  • Startseite
  • Kontakt
  • Impressum
  • Home
  • Listenauswahl
  • Anleitungen

[Seqan-dev] Trouble using SequenceStream

<-- thread -->
<-- date -->
  • From: Theodore Omtzigt <theo@stillwater-sc.com>
  • To: seqan-dev@lists.fu-berlin.de
  • Date: Thu, 24 Jan 2013 12:42:37 -0500
  • Organization: Stillwater Supercomputing, Inc.
  • Reply-to: SeqAn Development <seqan-dev@lists.fu-berlin.de>
  • Subject: [Seqan-dev] Trouble using SequenceStream

I am getting stuck on the simplest of tasks, reading an uncompressed
fastq file. This statement:
          seqan::SequenceStream seqIO(argv[1]);

returns a FILE_FORMAT_ERROR on a valid fastq text file. In
sequence_stream.h in the function _checkFormat, checkStreamFormat
returns false,

    // Guess file format, given the record reader.
    template <typename TStream, typename TSpec>
    SeqIOFileFormat_::Type _checkFormat(RecordReader<TStream, TSpec> &
recordReader)
    {
        AutoSeqStreamFormat formatTag;
        if (!checkStreamFormat(recordReader, formatTag))
            return SeqIOFileFormat_::FILE_FORMAT_ERROR;

The recursion in checkStreamFormat in guess_stream_format.h is confusing
me so I can't figure out what specific failure it is tripping over. I am
trying to debug this but need some help deciphering the recursion
architecture of the guessFormat code.

How do I know it is a valid fastq file? The file provided has been a
test file for non seqan code, and I have been able to use it in seqan
with the FragmentStore without trouble.




<-- thread -->
<-- date -->
  • Follow-Ups:
    • Re: [Seqan-dev] Trouble using SequenceStream
      • From: "Holtgrewe, Manuel" <manuel.holtgrewe@fu-berlin.de>
    • Re: [Seqan-dev] Trouble using SequenceStream
      • From: "zhoutianliangwen@genomics.org.cn" <zhoutianliangwen@genomics.org.cn>
  • seqan-dev - January 2013 - Archives indexes sorted by:
    [ thread ] [ subject ] [ author ] [ date ]
  • Complete archive of the seqan-dev mailing list
  • More info on this list...

Hilfe

  • FAQ
  • Dienstbeschreibung
  • ZEDAT Beratung
  • postmaster@lists.fu-berlin.de

Service-Navigation

  • Startseite
  • Listenauswahl

Einrichtung Mailingliste

  • ZEDAT-Portal
  • Mailinglisten Portal