[Bioperl-l] Fastq - is it flush, or not?

Nestor Zaburannyi nestor at linuxmail.org
Mon Jul 11 04:47:03 EDT 2011


Dear All.

I need to reverse-complement a huge *.fastq file. I try to

while (my $seq = $in->next_seq)
    {
    $out->write_seq($seq->revcom);
    }

However, it stops with:

------------- EXCEPTION -------------
MSG: Can not get a reverse complement. The object is not flush.
STACK Bio::Seq::Meta::Array::revcom /usr/share/perl5/Bio/Seq/Meta/Array.pm:648


Adding

$seq->force_flush('1');

helps, but is it necessary and proper way? Sequences seem to be "flush".

Sincerely yours
Nestor




More information about the Bioperl-l mailing list