[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