This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Extract sequence from Genebank

Hello, Can someone tell the way to extract the fasta sequences of gene 16S? Can I extract them in the genebank? Thanks

genebank
• 1,833 views • link
written by anasofiamoreira94 anasofiamoreira94

A particular species or just 16S in general?

You can get the 16S microbial blast pre-formatted database here. You can extract the sequences from these file using blastdbcmd utility from BLAST+ package. Note: These sequences are from RefSeq targeted loci project.

• 0 views • link
written by GenoMax GenoMax

I have a number of particular species and I have to look for its fasta sequence according to the 16s. So my idea was to select only the gene 16s from mitochondrion in genebank.Do you think it is wrong? Thanks for the help.

• 0 views • link
written by anasofiamoreira94 anasofiamoreira94

Can you provide an example? There may be an Entrezdirect way to get this information.

• 0 views • link
written by GenoMax GenoMax

For instance, I want a fasta sequence of 16s for Acetes japonicus

• 0 views • link
written by anasofiamoreira94 anasofiamoreira94

Using EntrezDirect: For that organism there are only 14 sequences in Genbank and none of them are 16S. If you can use the following sequences then remove the last | grep part and replace with > file.fa.

$ esearch -db nuccore -query "Acetes japonicus [ORGN]" | efetch -format fasta | grep "^>"
>HQ630575.1 Acetes japonicus isolate AJ18_f_KK19 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630574.1 Acetes japonicus isolate AJ13_f_KK19 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630573.1 Acetes japonicus isolate AJ12_f_KG26 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630572.1 Acetes japonicus isolate AJ11_m_KG26 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630571.1 Acetes japonicus isolate AJ10_f_KG26 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630570.1 Acetes japonicus isolate AJ8_f_KG26 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630569.1 Acetes japonicus isolate AJ7_m_KG26 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630568.1 Acetes japonicus isolate AJ6_f_Tbhg18 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630567.1 Acetes japonicus isolate AJ5_m_Tbhg cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630566.1 Acetes japonicus isolate AJ4_f_Tbhg cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630565.1 Acetes japonicus isolate AJ3_m_tbhg cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630564.1 Acetes japonicus isolate AJ2_f_Tbhg cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>HQ630563.1 Acetes japonicus isolate AJ1_m_Tbhg cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
>KF977240.1 Acetes japonicus isolate DNZ050 cytochrome oxidase subunit I (COI) gene, partial cds; mitochondrial
• 0 views • link
written by GenoMax GenoMax

Ok,I will try that. But I'm having trouble using esearch, when I try to execute it, I get this error: WebEnv value not found in search output - WebEnv1 Db value not found in fetch input

• 0 views • link
written by anasofiamoreira94 anasofiamoreira94

Exact command above or something else?

• 0 views • link
written by GenoMax GenoMax

I used this command: esearch -db pubmed -query "lycopene cyclase" |efetch -format abstract

• 0 views • link
written by anasofiamoreira94 anasofiamoreira94

While this is unrelated to the original question, that command looks correct. I can get abstracts of papers by running the same command. There is going to be a lot of them so be careful.

• 0 views • link
written by GenoMax GenoMax

0 answers

No answers yet.

Log in to answer this question.

More posts like this