Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BLAST+ creating custom blast database and using blast+ filtering features

    Hi,

    I would like to create a personal blast database of arbitrary sequences and be able to use all the features of BLAST+ to create subsets of databases based on identifiers or filter based on taxonomy.

    It looks like the formatting of the definition line in the input FASTA files is crucial to assign proper sequence identifiers.

    Using the General database identifier gnl|database|identifier or local identifier format lcl|identifier I wasn't able to use the blastdb_aliastool to create db subsets as it expects a GI list as input. I also didn't have any luck assigning taxonomy identifiers with the -taxid_map option of makeblastdb.

    What is the recommended way to format FASTA definition lines in order to be able to use all the filtering features of the BLAST+ tools.

    I was thinking of creating pseudo GenBank definitions for all my sequences: gi|<gi-number>|gb|<AccessionVersion>|<Accession>, where <gi-number> is a generated numeric value, and <Accession/Version> is my identifier. This works for the GI based filtering, however it seems like an ugly hack and I would prefer something more straight forward.

    How is the taxid_map file formatted? I've tried <gi-number>, <gb|<AccessionVersion>>, or <gb|Accession> as the sequence identifier, however they don't seem to be assigned properly and blastdbcmd with -outfmt %T gives me zero for all entries.

    Thanks for any help,

    Deniz

    PS: I've posted this already in http://I.SEQanswers.com but as it's still beta traffic is a bit low and I'm probably more likely to get an answer here in the forum.

  • #2
    Hi,
    I was wondering if you found a solution for this? I have been having all kinds of problems creating subset blast databases.
    Thanks

    Comment


    • #3
      GI lists are essentially simple to create, if you have experience with bash, perl, or any programming language that you can use to format text, you may be able to automatically pull the GIs from these files and put them into a GI list. A GI list is simply a text file with one number per line, and each number is a GI. There may be some utilities that do this automatically (i.e. FASTA->GI list), though I don't know of them.

      Alternatively, a nice way I've found for creating GI lists based on queries is to query NCBI for the data you want. You can do a search in NCBI, click send to on the top-right corner, and export it as a GI list. This is a nice, easy way of getting a GI list for blast subsets, though it's difficult to automate.

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Techniques and Challenges in Conservation Genomics
        by seqadmin



        The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

        Avian Conservation
        Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
        03-08-2024, 10:41 AM
      • seqadmin
        The Impact of AI in Genomic Medicine
        by seqadmin



        Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
        02-26-2024, 02:07 PM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 03-14-2024, 06:13 AM
      0 responses
      34 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-08-2024, 08:03 AM
      0 responses
      72 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-07-2024, 08:13 AM
      0 responses
      82 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-06-2024, 09:51 AM
      0 responses
      68 views
      0 likes
      Last Post seqadmin  
      Working...
      X