[BioRuby-cvs] bioruby/lib/bio/appl/blast format0.rb,1.25,1.26
Naohisa Goto
ngoto at dev.open-bio.org
Tue Feb 12 02:13:34 UTC 2008
Update of /home/repository/bioruby/bioruby/lib/bio/appl/blast
In directory dev.open-bio.org:/tmp/cvs-serv9861/lib/bio/appl/blast
Modified Files:
format0.rb
Log Message:
* Bug fixes reported by Shuji Shigenobu.
* Failed to parse query length for long query (>= 10000 letters)
because comma is inserted for digit separator by blastall.
* Failed to parse e-value for some BLASTX results
Index: format0.rb
===================================================================
RCS file: /home/repository/bioruby/bioruby/lib/bio/appl/blast/format0.rb,v
retrieving revision 1.25
retrieving revision 1.26
diff -C2 -d -r1.25 -r1.26
*** format0.rb 27 Dec 2007 17:28:57 -0000 1.25
--- format0.rb 12 Feb 2008 02:13:31 -0000 1.26
***************
*** 265,270 ****
q << sc.scan(/.*/)
sc.skip(/\s*^ ?/)
! end until !sc.rest or r = sc.skip(/ *\( *(\d+) *letters *\)\s*\z/)
! @query_len = sc[1].to_i if r
@query_def = q.join(' ')
end
--- 265,270 ----
q << sc.scan(/.*/)
sc.skip(/\s*^ ?/)
! end until !sc.rest or r = sc.skip(/ *\( *([\,\d]+) *letters *\)\s*\z/)
! @query_len = sc[1].delete(',').to_i if r
@query_def = q.join(' ')
end
***************
*** 969,973 ****
while sc.rest?
sc.skip(/\s*/)
! if sc.skip(/Expect(?:\(\d\))? *\= *([e\-\.\d]+)/) then
ev = sc[1].to_s
ev = '1' + ev if ev[0] == ?e
--- 969,973 ----
while sc.rest?
sc.skip(/\s*/)
! if sc.skip(/Expect(?:\(\d+\))? *\= *([e\-\.\d]+)/) then
ev = sc[1].to_s
ev = '1' + ev if ev[0] == ?e
More information about the bioruby-cvs
mailing list