This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Retrieve gene length and cds locus informartion

I have been trying to retrieve cds start and end position and gene length for about 1000 protein accessions (eg, ALD89117.1, ALD89128.1, ALD89126.1, ANR02692.1,AVA17449.1) I have as input. Would someone be kind enough to share code or expertise telling how to get this done? Thanks

ncbi batch_entrez
• 1,619 views • link
updated by GenoMax GenoMax • written by MAPK MAPK

2 answers 

$ esearch -db protein -query "AVA17449.1" | elink -target nuccore | efetch -format ft
>Feature gb|MG256173.1|
473 2857    CDS
            product RNA-dependent RNA polymerase
            transl_table    4
            protein_id  gb|AVA17449.1|

$ esearch -db protein -query "ALD89117.1" | elink -target nuccore | efetch -format ft
>Feature gb|KP900907.1|
2812    572 CDS
            product RNA-dependent RNA polymerase
            transl_table    4
            protein_id  gb|ALD89117.1|
• 1 views • link
written by GenoMax GenoMax

You don't even tell us what those identifiers are, which would make the task of providing advice easier.

There are plenty of solutions around, for example using GenomicRanges and the "org.eg.db" packages, or if you have a GTF you can follow Tutorial: Extract Total Non-Overlapping Exon Length Per Gene With Bioconductor, you can use biomaRt as well - see this comment for hints. Finally, although Obtaining Exon Lengths: is really old, I think the answers still work.

• 1 views • link
written by h.mon h.mon

Thanks, I have edited the question with accessions.

• 1 views • link
written by MAPK MAPK

Log in to answer this question.

More posts like this