[Bioperl-l] bp_extract_feature_seq.pl : --feature CDS error

Charles Hauser chauser at duke.edu
Thu Aug 28 15:18:26 EDT 2003


I updated my bioperl to cvs today and now am unable to extract CDS
features from a genbank report using the script Using
bp_extract_feature_seq.pl.


using:

 -i foo.gbk --format genbank --feature gene

script returns concatenated fasta data, but ends with:
Can't call method "isa" on an undefined value at
/usr/local/src/bioperl/core/Bio/SeqFeatureI.pm line 457, <GEN0> line
1781.

line 1781 is the '//'line at the end of the gbk report.


Attempting to extract CDS feature seqs  --feature CDS returns no data
just the error listed above.

Charles




More information about the Bioperl-l mailing list