[Bioperl-l] Problem in OMIM parser

Mark A. Jensen maj at fortinbras.us
Sat Dec 6 06:20:11 UTC 2008


Hi Neeti - 

I looked at the source, and found that that exception is thrown 
when there are not exactly 18 fields ( separated by 17 "|" symbols ) 
on a genemap file line. The Bioperl tests ran fine for OMIMparser 
on the test data file. I had a hunch, so I removed the final linebreak 
from the test file; that is, from this...
.
.
.
1.132|12|2|97|1p35.1|P28|P|Inner dynein arm, clamydomonas, ...[\n]
1.133|6|10|98|1p35|BAI2|P|Brain-specific angiongenesis inhibi...[\n]
1.134|8|21|91|1p35|EBVS1|P|Epstein-Barr virus integration sit...[\n]
1.135|3|21|93|1p35|G1P3, IFI616|C|Interferon, alpha-inducibl...[\n]
[EOF]

to this...
.
.
.
1.132|12|2|97|1p35.1|P28|P|Inner dynein arm, clamydomonas, ...[\n]
1.133|6|10|98|1p35|BAI2|P|Brain-specific angiongenesis inhibi...[\n]
1.134|8|21|91|1p35|EBVS1|P|Epstein-Barr virus integration sit...[\n]
1.135|3|21|93|1p35|G1P3, IFI616|C|Interferon, alpha-inducibl...[EOF]

and ran the test again. Immediately, I got

**
ok 1 - use Bio::Phenotype::OMIM::OMIMparser;

------------- EXCEPTION: Bio::Root::Exception -------------
MSG: Gene map file "t\data\omim_genemap_test" is not in the expected format
STACK: Error::throw
STACK: Bio::Root::Root::throw c:/Perl/site/lib/Bio/Root/Root.pm:374
STACK: Bio::Phenotype::OMIM::OMIMparser::_read_genemap c:/Perl/site/lib/Bio/Phen
otype/OMIM/OMIMparser.pm:883
STACK: Bio::Phenotype::OMIM::OMIMparser::genemap_file_name c:/Perl/site/lib/Bio/
Phenotype/OMIM/OMIMparser.pm:408
STACK: Bio::Phenotype::OMIM::OMIMparser::new c:/Perl/site/lib/Bio/Phenotype/OMIM
/OMIMparser.pm:224
STACK: t/Phenotype/OMIMparser.t:16
-----------------------------------------------------------
# Looks like you planned 174 tests but only ran 1.
# Looks like your test died just after 1.
**

just your exception. So you might begin by making sure you're not chopping the final
newline character from your input files, or visiting those files to make sure there is 
a linebreak after the final line.

Hope this a start.
cheers MAJ



More information about the Bioperl-l mailing list