Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • arkilis
    Senior Member
    • Jul 2013
    • 119

    #1

    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
  • maubp
    Peter (Biopython etc)
    • Jul 2009
    • 1544

    #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

    • kinpict.com
      Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by kinpict.com
      The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
      Today, 06:37 PM
    • SEQadmin2
      16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
      by SEQadmin2





      NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
      10-30-2026, 06:46 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 11-12-2026, 07:59 AM
    0 responses
    28 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 11-05-2026, 09:25 AM
    0 responses
    35 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 11-03-2026, 11:54 AM
    2 responses
    56 views
    0 reactions
    Last Post kinpict.com  
    Started by SEQadmin2, 10-28-2026, 09:57 AM
    0 responses
    175 views
    0 reactions
    Last Post SEQadmin2  
    Working...