[Bioperl-l] Bio::SeqFeature::Generic tag trouble.
Govind Chandra
govind.chandra at bbsrc.ac.uk
Fri Oct 28 09:17:17 EDT 2005
Hi,
The script below is giving me the output shown further below. Could
someone please point out what is it that I am doing wrong.
Perl is 5.8.7. BioPerl version is 1.5.0
Linux is RedHat 8.0
Thanks
Govind
Microbiology
John Innes Centre
Norwich NR4 7UH
UK
### script begins ###
use Bio::SeqIO;
$seqout=Bio::SeqIO->new('-fh' => \*STDOUT, '-format' => 'embl');
$outobj=Bio::Seq->new(
'-seq' => 'cccgcggagcgggtaccacatcgctgcgcgatgtgcaagcgaacacccgcgctgc');
$nft=Bio::SeqFeature::Generic->new(
'-start' => 10,
'-end' => 25,
'-strand' => -1,
'-primary' => 'CDS',
'-tag' => {'locus_tag' => 'something',
'gene' => 'hoaX'}
);
$outobj->add_SeqFeature($nft);
$seqout->write_seq($outobj);
### script ends ###
### output begins ###
ID standard; DNA; UNK; 55 BP.
XX
AC unknown;
XX
XX
FH Key Location/Qualifiers
FH
FT CDS complement(10..25)
FT /locus_tag="Bio::Annotation::SimpleValue=HASH(0x85daeac)"
FT /gene="Bio::Annotation::SimpleValue=HASH(0x85e1adc)"
XX
SQ Sequence 55 BP; 10 A; 21 C; 18 G; 6 T; 0 other;
cccgcggagc gggtaccaca tcgctgcgcg atgtgcaagc gaacacccgc gctgc 55
//
### output ends ###
More information about the Bioperl-l
mailing list