This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieve amino acid sequences from old locusID?

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?

blast genome sequencing sequence gene
• 2,116 views • link
updated by Sparrow_kop Sparrow_kop • written by taylrtam • taylrtam

2 answers

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.

• 0 views • link
written by Sparrow_kop Sparrow_kop

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.

• 0 views • link
written by Sparrow_kop Sparrow_kop

Log in to answer this question.

More posts like this