[Bioperl-l] extracting 'comments' from gff-ver 2

Charles Hauser chauser at duke.edu
Fri Apr 4 10:20:50 EST 2003


I am trying to generate a plot of the distribution of ESTs on genomic
scaffolds using the GFF output from exonerate.  I need to grab the
clone_id  in the comments field '833005F02.y1'. 

the method $feature->comments does not exist, how can I access these
data?  

-regards

Charles




use Bio::Tools::GFF;
<snip>
my $gffio = Bio::Tools::GFF->new(-fh => \*FH,
				 -gff_version => 2);
my $feature;
while($feature = $gffio->next_feature()) {
    if($feature->primary_tag eq 'gene') {
	print 
	    $feature->seqname,  "\n",
	    $feature->comments,  "\n";



##gff-version 2
##source-version exonerate:est2genome 0.6.7
##date 2003-03-24
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
scaffold_1	exonerate:est2genome	gene	367897	368609	1884	+	.	gene_id	1	;	sequence	833005F02.y1	;	gene_orientation	
+	;






More information about the Bioperl-l mailing list