Question: Retrieve amino acid sequences from old locusID?
0
gravatar for taylrtam
2.1 years ago by
taylrtam0
taylrtam0 wrote:

Hi,

I have a list of discontinued locusID tags that I want to retrieve the amino acid sequences for.

The list is as follows:

Bphyt_0466

Bphyt_4414

...

Bphyt_4439

Is there any way to batch retrieve the amino acid sequences?

ADD COMMENTlink modified 2.1 years ago by Sparrow_kop220 • written 2.1 years ago by taylrtam0
1
gravatar for Sparrow_kop
2.1 years ago by
Sparrow_kop220
China
Sparrow_kop220 wrote:

Hi, you can use R 'KEGGREST' library to extract the amino seq from the KEGG database: The below is the demo script:

library(KEGGREST)
library(Biostrings)

aaseq_query <- keggGet("bpy:Bphyt_2906", "aaseq")
aaseq <- as.character(unlist(aaseq))
aaseq

[1] "MSEKEIDQVLVERVQKGDKAAFELLVSKYHRKILRLISRLVRDPAEVEDVAQDAFIKAYRALPQFRGESAFYTWLYRIAVNTAKNYLATQGRRAPTSTEADAEEAETFSDADQLRDINTPESMLMSKQIAETVNAAMAVLPEELRTAITLREIEGLSYEEIAEMMGCPIGTVRSRIFRAREAIAAKLRPLLDTPEGKRW"

The other way is to use web scrapping by python or others, use the below url :

http://www.genome.jp/dbget-bin/www_bget?bpy:Bphyt_2906

and you can replace the gene id at the last position in the url , and the write a script to parse the web to extract the seq.

ADD COMMENTlink written 2.1 years ago by Sparrow_kop220
0
gravatar for Sparrow_kop
2.1 years ago by
Sparrow_kop220
China
Sparrow_kop220 wrote:

Hi, you can use R 'KEGGREST' library to extract the amino seq from the KEGG database: The below is the demo script:

library(KEGGREST)
library(Biostrings)

aaseq_query <- keggGet("bpy:Bphyt_2906", "aaseq")
aaseq <- as.character(unlist(aaseq))
aaseq

[1] "MSEKEIDQVLVERVQKGDKAAFELLVSKYHRKILRLISRLVRDPAEVEDVAQDAFIKAYRALPQFRGESAFYTWLYRIAVNTAKNYLATQGRRAPTSTEADAEEAETFSDADQLRDINTPESMLMSKQIAETVNAAMAVLPEELRTAITLREIEGLSYEEIAEMMGCPIGTVRSRIFRAREAIAAKLRPLLDTPEGKRW"

The other way is to use web scrapping by python or others, use the below url :

http://www.genome.jp/dbget-bin/www_bget?bpy:Bphyt_2906

and you can replace the gene id at the last position in the url , and the write a script to parse the web to extract the seq.

ADD COMMENTlink written 2.1 years ago by Sparrow_kop220
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1366 users visited in the last hour