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

Re: [Seqan-dev] Reading fastq records - missing records?

<-- thread -->
<-- date -->
  • From: Alexander Vowinkel <vowinkel.alexander@gmail.com>
  • To: SeqAn Development <seqan-dev@lists.fu-berlin.de>
  • Date: Wed, 3 Jun 2015 09:52:46 -0500
  • Reply-to: SeqAn Development <seqan-dev@lists.fu-berlin.de>
  • Subject: Re: [Seqan-dev] Reading fastq records - missing records?

Hi Fadhlan,

why should memory be an issue?
I can load 400 of 10000 sized batches.
I always clear the Sets that save the sequences
after they were used.

Best,
Alexander

2015-06-02 23:36 GMT-05:00 Fadhlan Rizal <fadhlan.fari@gmail.com>:
i dont know why, but maybe because your memory doesn't enough?

2015-06-01 23:14 GMT+07:00 Alexander Vowinkel <vowinkel.alexander@gmail.com>:
Hi Team,

here is another question:
Following code reads a fastq-file.
The file I'm testing on has 263333413 records:
$ zcat C563WACXX_1.fastq.gz | wc -l
1053333652

The problem is, that the code seems to only read
the first 4000000 records. I don't understand why.

Can you advise me here?

Thank you,
Alexander

PS: I removed try/catch for better readability

seqan::StringSet<seqan::CharString> ids;
seqan::StringSet<seqan::Dna5QString> seqs;
seqan::StringSet<seqan::CharString> quals;
unsigned batchSize = 10000;
do {
seqan::clear(ids); seqan::clear(seqs); seqan::clear(quals);
readRecords(ids, seqs, quals, inputFile, batchSize);
processRecords(ids, seqs, quals, barcodeGraph, vertexFilePtrMap);
totalProcessed += seqan::length(ids);
std::cout << "\r";
std::cout << "reads processed in total: " << totalProcessed;
std::cout.flush();
} while(seqan::length(ids) > 0);

_______________________________________________
seqan-dev mailing list
seqan-dev@lists.fu-berlin.de
https://lists.fu-berlin.de/listinfo/seqan-dev



_______________________________________________
seqan-dev mailing list
seqan-dev@lists.fu-berlin.de
https://lists.fu-berlin.de/listinfo/seqan-dev


<-- thread -->
<-- date -->
  • References:
    • [Seqan-dev] Reading fastq records - missing records?
      • From: Alexander Vowinkel <vowinkel.alexander@gmail.com>
    • Re: [Seqan-dev] Reading fastq records - missing records?
      • From: Fadhlan Rizal <fadhlan.fari@gmail.com>
  • seqan-dev - June 2015 - 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