Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Could I get the species name list from a blast database?

    I was thinking is there anyway to get the species name list from a blast db, I saw they are binary files.

    I found this:

    blastdbcmd -db [db name] -info, but actually I don't have the blastdbcmd installed on our server.

    Any thoughts? Cheers,

    ======================================================
    Here is what I tried on blastdbcmd, no species name found anyway.

    Database: Kog
    4,825 sequences; 2,336,026 total residues

    Date: Jun 2, 2009 10:24 AM Longest sequence: 5,019 residues

    Volumes:
    /usr/local/agrf/data/ncbi-blast/linked_databases/Kog
    Last edited by arkilis; 05-06-2014, 06:32 PM. Reason: more details

  • #2
    BLAST databases can record the species (via an NCBI taxid) for each sequence, so there can be many species in a single database. You can use the out and outfmt options to blastdbcmd to get this information:

    Code:
     -outfmt <String>
       Output format, where the available format specifiers are:
       		%f means sequence in FASTA format
       		%s means sequence data (without defline)
       		%a means accession
       		%g means gi
       		%o means ordinal id (OID)
       		%i means sequence id
       		%t means sequence title
       		%l means sequence length
       		%h means sequence hash value
       		%T means taxid
       		%e means membership integer
       		%L means common taxonomic name
       		%S means scientific name
       		%K means taxonomic super kingdom
       		%P means PIG
       		%m means sequence masking data.
       		   Masking data will be displayed as a series of 'N-M' values
       		   separated by ';' or the word 'none' if none are available.
       	If '%f' is specified, all other format specifiers are ignored.
       	For every format except '%f', each line of output will correspond
       	to a sequence.
    e.g.

    Code:
    $ blastdbcmd -db pdbaa -outfmt "%a %T %K %S %L" -entry all | head
    1VS5_R 562 Bacteria Escherichia coli Escherichia coli
    1VS7_R 562 Bacteria Escherichia coli Escherichia coli
    3I1M_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3I1O_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3I1Q_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3I1S_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3I1Z_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3I21_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3KC4_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    3OR9_R 83333 Bacteria Escherichia coli K-12 Escherichia coli K-12
    For this to work, in addition to the BLAST database, you also need the taxonomy information. Download ftp://ftp.ncbi.nlm.nih.gov/blast/db/taxdb.tar.gz then unzip this next to your other BLAST databases.
    Last edited by maubp; 05-07-2014, 02:41 AM. Reason: Adding example

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin


      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
      Yesterday, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    39 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    41 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    35 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    55 views
    0 likes
    Last Post seqadmin  
    Working...
    X