Hi all,
background: I need to be able to recognize all possible sequence identifiers present in preformatted NCBI nucleotide databases. I've implemented regular expression following https://www.ncbi.nlm.nih.gov/Sequin/acc.html, but it is not enough. Other accessions (e.g. PDB) are also present. So I would like to have examples of all possible formats I can encounter. But I was not able to find any list which would describe what actually can be inside those databases.
One possible solution, I thought would be to use ENTREZ to retrieve the accessions for me. There is blastdbinfo database which lists the avalible databases. But I not able to get elink to link anywhere.
Lets for example focus on refseq_genomes.
The database is available with following command:
esearch -query refseq_genomes[DB] -db blastdbinfo
So given that I want nucleotide sequence accessions present in that database what the elink statement should be?
esearch -query refseq_genomes[DB] -db blastdbinfo | ... SOME ELINK .... | efetch --format acc
For ENTREZ experts here - How do I tell which database links where?
I know I can download the databases and use blastdbcmd to obtain the accessions, but It should be possible to obtain the accessions in some better way.
Thank you
1 answer
To anybody who might want to retrieve accessions from BLAST database without downloading whole database.
I didn't find any route through ENTREZ, however, I've found that you don't need the whole BLAST database to retrieve the accessions for the sequences.
Only the .nhr, .nin and .nal (.ndb if .nal is not present) files are required to call the blastdbcmd -entry all -outfmt "%a".
According to this documentation (http://nebc.nerc.ac.uk/bioinformatics/documentation/blast/formatdb.html)
The .nhr are headers for the fasta files and .nin are the indices. The nal file describes the database.
Warning - to my knowledge this is not documented and can stop working at any time. Tested on db v4 and v5.
Yes I agree, I've also tested it on db v5.
Log in to answer this question.
-
Odd behaviour - NCBI Eutils
written by npb27 •
-
What is wrong with my NCBI esearch command for 16S rRNA?
written by Morgan S. •
-
Get corresponding SRA accessions for very large list of Biosample accesions?
written by goatsrunfaster •
-
How can I find genomes that have linked transcriptome data ? (bacteria)
written by hugo.avila
-
efetch results incomplete
written by j.j.e.vanhooff •
-
Using Entrez to find the taxonomy for an accession number
written by Frieda •
-
E-utilities for obtain gene sequences from the gene database
written by dllopezr •
-
NCBI esearch on multiple database
written by peteladrien •
-
Fetching multiple FASTA from NCBI
written by mikael.lenz.strube •
-
Running Entrez Direct-like commands on local or on custom mirror
written by p6patel •
For a given
db, you can find all available link names and a brief description as follows:The Entrex Link Descriptions webpage also lists this information but I am not sure how up-to-date that is. It looks like
blastdbcmdmay be the best solution for you.Thank you for the link. According to that it looks like there is no direct link between
blastdbinfoand e.g.nuccore.