[Bioperl-l] need help urgently
Neeti Somaiya
neetisomaiya at gmail.com
Mon Sep 7 10:04:06 UTC 2009
I tried using EntrezGene instead of GenBank, as is given in the link
that you sent :
http://www.bioperl.org/wiki/HOWTO:Beginners#Retrieving_a_sequence_from_a_database
http://doc.bioperl.org/releases/bioperl-current/bioperl-live/Bio/DB/EntrezGene.html
use Bio::DB::EntrezGene;
my $db = Bio::DB::EntrezGene->new;
my $seq = $db->get_Seq_by_id(2); # Gene id
# or ...
my $seqio = $db->get_Stream_by_id([2, 4693, 3064]); # Gene ids
while ( my $seq = $seqio->next_seq ) {
print "id is ", $seq->display_id, "\n";
}
This doesnt seem to work.
-Neeti
Even my blood says, B positive
On Fri, Sep 4, 2009 at 1:09 PM, Emanuele Osimo<e.osimo at gmail.com> wrote:
> Hello,
> have you tried this?
> http://bio.perl.org/wiki/HOWTO:Getting_Genomic_Sequences#Using_Bio::DB::GenBank_when_you_have_genomic_coordinates
>
> Emanuele
>
> On Fri, Sep 4, 2009 at 08:49, Neeti Somaiya <neetisomaiya at gmail.com> wrote:
>>
>> Hi,
>>
>> I have an input list of gene names (can get gene ids from a local db
>> if required).
>> I need to fetch sequences of these genes. Can someone please guide me
>> as to how this can be done using perl/bioperl?
>>
>> Any help will be deeply appreciated.
>>
>> Thanks.
>>
>> -Neeti
>> Even my blood says, B positive
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
>
More information about the Bioperl-l
mailing list