[Bioperl-l] patch: fix parsing of FASTA results

Aaron Mackey amackey at virginia.edu
Wed Jun 9 14:50:48 UTC 2010


Hopefully someone a bit more git-savvy than me can drop this little ditty in
to Bio/SearchIO/fasta.pm?

[ajm6q at lc4 SearchIO]$ diff -u fasta.pm~ fasta.pm
--- fasta.pm~    2010-04-09 12:59:00.000000000 -0400
+++ fasta.pm    2010-06-09 10:46:11.350773222 -0400
@@ -623,7 +623,7 @@
             $_ = $self->_readline();
             my ( $score, $bits, $e ) = /Z-score: \s* (\S+) \s*
                                (?: bits: \s* (\S+) \s+ )?
-                               (?: E|expect ) \s* \(\) :? \s*(\S+)/ox;
+                               (?: E|expect ) \s* \((?:\d+)?\) :?
\s*(\S+)/ox;
             $bits = $score unless defined $bits;

             my $v = shift @hit_signifs;



More information about the Bioperl-l mailing list