This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Entrez Direct to retrieve Protein and Coding sequences from NCBI accesion

Hi everyone!

I have a bunch of fasta files of gene sequences that I download from ncbi trough entrez direct tool. I am wondering if it is posible to obtain the protein and coding sequences of these genes using the accession, that is in this format: NZ_CP006694.1:1104181-1105143

where the data following the : is the sequence section where the gene is located.

Can you help me with that?

Thank you so much

entrez ncbi coding sequences retrieve protein
• 2,551 views • link
updated by vkkodali_ncbi vkkodali_ncbi • written by dllopezr • dllopezr

An R-based solution would be to use the Bioconductor package BiomaRt; please see my post here. Since you have the exact chromosomal position already, you can easily covert this to sequences. You can find the appropriate filters (chromosome, start and end position) using listFilters(ensembl), and the attribute (protein / dan sequence) using listAttributes(ensembl).

• 1 views • link
written by thomaskuilman thomaskuilman

1 answer

I think you can use Edirect for this as follows:

efetch -db nuccore -id 'NZ_CP006694.1' -seq_start 1104181 -seq_stop 1105143 -format fasta_cds_aa
• 1 views • link
written by vkkodali_ncbi vkkodali_ncbi

Log in to answer this question.

More posts like this