hi,
when u click at the link which i mentioned, you will find:
'search [Nucleotide] for [ ]'
when u'll paste a gi id e.g. gi|5104821|dbj|AP000061.1| in space after 'for',
a page will open where u will see highlighted accession id: AP000061.
when u'll click at this, a default information page will open which contain gene and sometimes chromosome name for this sequence.in this case, many genes r displayed, n all are required.
there cases when even gene name is not available......
what i need is these genes n/or chromosome names if available if i supply a list of such gi ids' and get info in a list in the same order as that of gi ids.....
could u plz help me?
Before going any further, have you looked at <http://eutils.ncbi.nlm.nih.gov/entrez/query/static/eutils_help.html>? There are some perl scripts there that demonstrate how to query the database. Do any of them do what you want or something similar?
Regards, Randy.
PS: The program below parses the data into the fields you require.
--
#!/usr/bin/perl
use strict; use warnings;
use Data::Dumper;
my @data;
while (defined(my $line = <DATA>)) {
chomp( $line );
$line =~ s/^\s+//;
$line =~ s/\s+$//;
next unless length( $line );
my %fields;
@fields{ qw(query_id
subject_id
percent_identity
alignment_length
mismatches
gap_openings
q_start
q_end
s_start
s_end
e_value
bit_score
)
} = split( /\s+/, $line );
push( @data, \%fields );
}print Dumper( [EMAIL PROTECTED] );
__DATA__
gi|37182815|gb|AY358849.1| gi|28592069|gb|U63637.2|BTU63637 100.00 17 0 0 552 568 3218 3234 1.1 34.19
gi|37182815|gb|AY358849.1| gi|14318385|gb|AC089993.2| 95.24 21 1 0 435 455 56604 56624 1.1 34.19
gi|37182815|gb|AY358849.1| gi|14318385|gb|AC089993.2| 100.00 16 0 0 260 275 89982 89967 4.2 32.21
gi|37182815|gb|AY358849.1| gi|7385112|gb|AF222766.1|AF222766 100.00 17 0 0 345 361 242 226 1.1 34.19
-- To unsubscribe, e-mail: [EMAIL PROTECTED] For additional commands, e-mail: [EMAIL PROTECTED] <http://learn.perl.org/> <http://learn.perl.org/first-response>
