[Bioperl-l] Question about Fgenesh.pm
Andrew Nunberg
anunberg at oriongenomics.com
Fri Jul 23 12:58:09 EDT 2004
Hi,
I am trying to parse an fgenesh analysis file using Bio::Tools::Fgenesh from
bioperl-live
I had to make a change in the module at line 320
from:
$predobj->primary_tag($ExonTags{$flds[3]} . 'Exon');
to:
$predobj->primary_tag($ExonTags{$flds[4]} . 'Exon');
For the parser to work.
I then wish to retreive the name of the dna sequence that was fed to the
parser using the analysis_subject/query method
$fgenesh = Bio::Tools::Fgenesh->new(-file=>"results file");
$gene = $fgenesh->next_prediction;
$fgenesh->analysis_subject;
This method returns undef, although it appears in the code that the info is
stored there.
Lines 366-373 from Fgenesh.pm
if(/^(FGENESH)\s+([\d\.]+)/) {
$self->analysis_method($1);
$self->analysis_method_version($2);
if (/\s(\S+)\sgenomic DNA/) {
$self->analysis_subject($1);
}
next;
}
Any input would be appreciated
--
Andrew Nunberg
Bioinformagician
Orion Genomics
(314)-615-6989
www.oriongenomics.com
More information about the Bioperl-l
mailing list