[Bioperl-l] Build new Swiss-Prot Seq
hafiz hafiz
khufaz83 at yahoo.com
Thu Sep 16 07:59:48 EDT 2004
Hello , everybody can help me, i want to buid a new
sequence in swissprot and i have used 3 modul such as
Bio::SeqFeature, Bio::Seq, Bio::Annotation::Comments,
Bio::Annotation::Simplevalue, Bio::Annotation
Bio::Annotation::references and
Bio::Annotation::Dblink, so its true or not.
This is my souces code;
$feat = new Bio::SeqFeature::Generic ( -start => 10,
-end => 100,
-strand => -1,
-primary => 'repeat',
-source_tag =>
'repeatmasker',
-display_name => 'alu
family',
-score => 1000,
-tag => {
new => 1,
author =>
'someone',
sillytag => 'this
is silly!' } );
#my $feat = new Bio::SeqFeature::Generic(-start
=>10,-end => 22, -strand =>1,
# -tag =>{
author => 'john',note =>'TATA box'});
$feat->add_tag_value("match1","PF000123 e-7.2");
$feat->add_tag_value("match2","PF002534 e-7.2");
my @arr = $feat->get_all_tags;
for my $tag (@arr){
print $tag,":",$feat->get_tag_values($tag),"
";
}
my $seq_obj =Bio::Seq
->new(-seq=>"attcccccchjkhjkhkjh",
-display_id=>"BIO52" );
$seq_obj->add_SeqFeature($feat);
my $comment =Bio::Annotation::Comment->new;
$comment->text("this looks like a good TATA box");
my $coll =new Bio::Annotation::Collection;
$coll->add_Annotation('comment',$comment);
$seq_obj->annotation($coll);
my $col = new Bio::Annotation::Collection;
my $sv = new Bio::Annotation::SimpleValue(-value =>
'someval');
$col->add_Annotation('tagname', $sv);
$reg = Bio::Annotation::Reference->new( -title =>
'title line',
-location
=> 'location line',
-authors
=> 'author line',
-medline
=> 998122 );
$link1 = new Bio::Annotation::DBLink(-database =>
'TSC',
-primary_id =>
'TSC0000030'
);
#or
$link2 = new Bio::Annotation::DBLink();
$link2->database('dbSNP');
$link2->primary_id('2367');
# DBLink is-a Bio::AnnotationI object, can be added
to annotation
# collections, e.g. the one on features or seqs
$feat->annotation->add_Annotation('dblink', $link2);
$file3="/home/database/Bioperl/test.gb";
my $io= Bio::SeqIO ->new(-format =>"swiss",-file =>
">$file3");
$io ->write_seq($seq_obj);
________________________________________________________________________
Yahoo! Messenger - Communicate instantly..."Ping"
your friends today! Download Messenger Now
http://uk.messenger.yahoo.com/download/index.html
More information about the Bioperl-l
mailing list