How to use efetch to download specific gene from whole genome NCBI by linux terminal?
1
0
Entering edit mode
2.8 years ago
Ahmad ▴ 10

I need download specific gene from whole genome (like gene length from 11200 bp to 11350 bp) on linux terminal.

efetch NCBI terminal • 1.2k views
ADD COMMENT
0
Entering edit mode

You already have competent answers to your question, and this is not one of them. If I had to do the same for 500+ genes, I'd probably be looking for a similar solution as you have, because manual solutions would be so time-consuming. But for a single gene, I so would not trust the automated process because there are always things that can wrong. Unless this is for homework, how difficult is it to find a single gene on NCBI and manually excise the range of bases you need?

ADD REPLY
0
Entering edit mode

I have large number of genes, manually is difficult

ADD REPLY
0
Entering edit mode
2.8 years ago
vkkodali_ncbi ★ 3.8k

You can use EntrezDirect for this as follows:

$ esearch -db gene -query '11200:11350[gene length] AND "homo sapiens"[organism]'

This returns 192 genes that satisfy this criteria. You can then pipe the output of the above command to other EntrezDirect tools depending on your downstream analysis.

ADD COMMENT

Login before adding your answer.

Traffic: 1835 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6