[Bioperl-l] Fastq - is it flush, or not?
Nestor Zaburannyi
nestor at linuxmail.org
Mon Jul 11 08:47:03 UTC 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