[Bioperl-l] Bio::Seq::Quality description line problem
David Messina
dmessina at wustl.edu
Fri May 5 15:24:47 UTC 2006
Apologies if this is a repost -- mail troubles this morning.
Hilmar is correct.
From a cursory walk through the code in a debugger, it looks like
Bio::SeqIO::qual's write_seq method doesn't read the 'desc' out of
the Bio::Seq::Quality object.
I think there should be something like this:
if ($source->can('desc') and my $desc = $source->desc()) {
$desc =~ s/\n//g;
}
$header .= " $desc";
before line 218 in Bio::SeqIO::qual (where the header is printed):
$self->_print (">$header \n");
Dave
More information about the Bioperl-l
mailing list