Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Convertig blastplus to blast output

    Hi all,

    here is the next blast problem (and I'm actually not sure there is a solution for this). I need to perform a blast against nr with some species/gis excluded. Formatdb offers the possibility to create a database subset based on a positive gi list; however, if this list is huge, the blast takes very long.
    The blast+ application offers both a positive and a negative gi list argument and performs much more efficient.
    However, I use some programs which expect the old plain-text blastall output as input and fail to parse the blast+ output (which differs in some details from blastall output).

    So, this is my question: does somebody know of a possibility to convert the blast+ output to the old blast output format?

    Thanks in advance and best regards!

  • #2
    There is a reason that for years (*) the NCBI (and projects like BioPerl and Biopython) have been telling people not to try to parse the plain text BLAST output - the NCBI have a long history of tweaking the text presentation and thus breaking parsers.

    Can you not use the tabular or XML output formats?

    (*) I think it must be years by now...

    Comment


    • #3
      Thanks for your reply. I think you are perfectly right. It's just that some of the programs I use seem to expect plain text BLAST output. I might try to tweak the BLAST+ output so the parser can read it...

      However, I'm having trouble with BLAST+ right now.

      blastx -db nr -frame_shift_penalty 15 -query myquery.fa
      throws an exception:

      Error: NCBI C++ Exception:
      "/am/ncbiapdata/release/blast/src/2.2.25/Linux64-Suse-icc/c++/ICC1010-ReleaseMT64--Linux64-Suse-icc/../src/objects/seq/../seqalign/Seq_align.cpp", line 670: Error: CreateDensegFromStdseg(): Std-seg segment lengths not accurate!
      The query is:

      >myquery
      GTTCTACATAAGCAACAGGTATATGCCCAACGTTATTATTCTTCTATATAAAAGATAAGCAAATC
      There is no exception if there is no frame-shift-penalty argument or if the sequence is different (e.g. the last 5 nucleotides removed). Someone knows what I'm doing wrong?

      Comment


      • #4
        Originally posted by sammy07 View Post
        Thanks for your reply. I think you are perfectly right. It's just that some of the programs I use seem to expect plain text BLAST output. I might try to tweak the BLAST+ output so the parser can read it...
        I remember having to tweak the Biopython parser to take BLAST+ plain text output.
        Originally posted by sammy07 View Post
        However, I'm having trouble with BLAST+ right now.
        ...
        Someone knows what I'm doing wrong?
        I've not seen that before. I'm trying the example here (although it is likely our copies of the NR database won't be exactly the same).

        Update - Same Error: CreateDensegFromStdseg(): Std-seg segment lengths not accurate! here from BLAST 2.2.25+ under both Linux 32 bit, Linux 64 bit and Mac OS X.

        You should report this to the NCBI.
        Last edited by maubp; 08-11-2011, 06:17 AM.

        Comment


        • #5
          Originally posted by maubp View Post
          You should report this to the NCBI.
          I've reported this by email.

          Update: Wayne Matten at the NCBI replied this was a known bug, and the BLAST developers are still working on a solution.
          Last edited by maubp; 08-15-2011, 12:51 AM.

          Comment


          • #6
            Great, thanks a lot!!

            Comment


            • #7
              BLASTX 2.2.26+
              ...
              Error: CreateDensegFromStdseg(): Std-seg segment lengths not accurate!


              Doesn't seem that it has changed with blast 2.2.26+. It's a pity...

              Comment


              • #8
                The problem was elegantly solved by the blast release 2.2.27+. The parameter -frame_shift_penalty no longer exists.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Best Practices for Single-Cell Sequencing Analysis
                  by seqadmin



                  While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                  06-06-2024, 07:15 AM
                • seqadmin
                  Latest Developments in Precision Medicine
                  by seqadmin



                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                  Somatic Genomics
                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                  05-24-2024, 01:16 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Today, 07:23 AM
                0 responses
                8 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-17-2024, 06:54 AM
                0 responses
                11 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-14-2024, 07:24 AM
                0 responses
                24 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 06-13-2024, 08:58 AM
                0 responses
                17 views
                0 likes
                Last Post seqadmin  
                Working...
                X