[Bioperl-l] ace.pm
Wes Barris
wes.barris at csiro.au
Mon Feb 16 18:19:02 EST 2004
Hi,
ACE files generated by an application called tgicl have "CO"
lines of the form:
CO CL15Contig2 794 4 0 U
This line is not parsed properly by the ace.pm bioperl module.
Notice this line from Bio/Assembly/IO/ace.pm .
(/^CO Contig(\d+) (\d+) (\d+) (\d+) (\w+)/) && do { # New
contig found!
Bioperl expects the second "word" in the line to be "Contig\d+" where
the number is used as the "contigID". Is there a reason why
"contigID" must be a number? Why can't it be the whole second
"word" of the "CO" line?
--
Wes Barris
E-Mail: Wes.Barris at csiro.au
More information about the Bioperl-l
mailing list