in reply to Re^2: .qual File Writing Complication
in thread .qual File Writing Complication
Just to note, as others have indicated putting the code here that generated that stack trace does help to trace the problem down. Looks like the problem arises when you are creating a new Bio::Seq::Quality, on line 165 of the code you submitted.
According to the Bio::Seq::Quality you are using the wrong named attribute for the quality information:
my $new_qual = Bio::Seq::Quality->new(-seq => $qual_string, -id => $seq->id(), -desc => $seq->desc);That should be:
my $new_qual = Bio::Seq::Quality->new(-qual => $qual_string, -id => $seq->id(), -desc => $seq->desc);This class also accepts sequence data to go with those scores; not sure if it requires it, though.
|
|---|