[BioRuby] EMBL parsing
Anthony Underwood
email2ants at gmail.com
Wed May 2 16:38:31 UTC 2007
Dear all,
I am having a problem when parsing EMBL genome files.
embl = open(filename)
entry = embl.gets(Bio::EMBL::DELIMITER)
seq_obj = Bio::EMBL.new(entry)
puts seq_obj.sequence_length
This takes a long time (minutes) and reports a sequence_length of 0
When taking an equivalent genbank file and changing the code
appropriately it parses the file in seconds and reports the correct
length. I am new to bioruby having used bioperl until now. Please can
anybody let me know if they have had similar problems and any
possible solutions.
Many thanks
Anthony
More information about the BioRuby
mailing list