Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Processing Blast output for Blast2GO

    Hello all,

    I've been trying to figure out how the taxonomy ids are set up for NCBI Blast results, but have been having some trouble. I am trying to figure out how to include/exclude results from Blast searches based on the associated taxa that Blast results associated with each sequence blasted. In other words, I want to be able to pull out all the sequences that blast to E. coli or Stronglyocentros purpuratus or, more generally, Eukaryotes, whatever, and then be able to import them into a program like Blast2GO to run some stats on them. I thought perhaps I could pull out a series of sequence names to use for filtering using some unix or python commands, but I haven't been able to figure out how to use the taxonomy ids to do that. Has anyone done this before with a decent sized data set (a few thousand sequences) and have some good suggestions?

    Thanks in advance any advice.

  • #2
    You could use MEGAN, http://www-ab.informatik.uni-tuebing...oftware/megan/. It allows for selecting a part of the taxonomic tree, and fetching all the reads that have a hit to that part (select a node, 'select subtree' from the menu, export the reads through the menu). Drawback is that it needs the blast output file as input, but you can restrict the output to just the first x alignments and no description lines to keep the file smaller.

    Comment


    • #3
      I'll check it out. Thanks for the suggestion. I was able to come up with long way around this issue, but would much rather find a way to figure out the taxonomy codes so I could either code a solution or find one already out there. Does anyone know of where I could find a good description of how the hierarchy of the NCBI taxon ids is set-up? I looked around the website and couldn't find much that was useful.

      Comment


      • #4
        Originally posted by flxlex View Post
        You could use MEGAN, http://www-ab.informatik.uni-tuebing...oftware/megan/. It allows for selecting a part of the taxonomic tree, and fetching all the reads that have a hit to that part (select a node, 'select subtree' from the menu, export the reads through the menu). Drawback is that it needs the blast output file as input, but you can restrict the output to just the first x alignments and no description lines to keep the file smaller.
        Mine gives me an error

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • 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

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 06:37 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        67 views
        0 likes
        Last Post seqadmin  
        Working...
        X