Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Basic question about local BLASTP

    Hello.

    I am using local BLASTP program to blast my protein sequences (around 600 sequences) against other organism's whole proteome database. What I want to do is to filter the output by some parameters. I have applied a reciprocal best_hit strategy, and the parameters were set as 'max identity > 30%', e<1e-10 and best_hit_overhang 0.25). However, in some cases, for example, the max identity between one best-hit set of query and subject sequence was as high as 40%, the query coverage was only 2%. Therefor, I want to exclude the best-hit sequences like that by setting a parameter like 'query coverage > 30%'.

    The web based blast on NCBI could show the query coverage, but I do not know what is the command line for local BLASTP to output "query coverage". Please give me some ideas to resolve this problem. Thanks!
    Last edited by Tsuyoshi; 04-22-2013, 05:18 PM.

  • #2
    You can try this:

    Update to the newest BLAST+ package.
    Use
    Code:
    -outfmt
    , which allows the output to be in custom format. You can find documentation on custom output format here (BLAST Command Line Applications User Manual, Appendix C, option 'outfmt'). From what I see the values to use are
    Code:
    -outfmt 6 qcovs
    , which should add query coverage per sequence value at the end of each output line and save it in tab-delimited format.

    Comment


    • #3
      Originally posted by mknut View Post
      You can try this:

      Update to the newest BLAST+ package.
      Use
      Code:
      -outfmt
      , which allows the output to be in custom format. You can find documentation on custom output format here (BLAST Command Line Applications User Manual, Appendix C, option 'outfmt'). From what I see the values to use are
      Code:
      -outfmt 6 qcovs
      , which should add query coverage per sequence value at the end of each output line and save it in tab-delimited format.
      Thank you mknut.

      I tried adding qcovs, however, the output files did not contain the qcovs information. For example, if I input the command as outfmt "qacc sacc evalue bitscore qcovs", the output files only contained qacc sacc evalue bitoscore. I am thinking whether qcovs is only available for blast but not for blastp.

      Comment


      • #4
        Originally posted by Tsuyoshi View Post
        I am thinking whether qcovs is only available for blast but not for blastp.
        The manual seems to indicate that those options are available for all blast+ programs. Have you tried 'std qcovs'?

        Or perhaps the command needs to be (note: single quotes):

        Code:
        -outfmt 6 [COLOR="Red"]'[/COLOR]qcovs[COLOR="Red"]'[/COLOR]
        Last edited by GenoMax; 04-23-2013, 04:47 AM.

        Comment


        • #5
          Table C1: Options common to all BLAST+ search applications.
          So yes, it seems that this command should be applicable to BLASTp. It also seems that the proper syntax to use is with double quotes, containing numerical argument as well as the rest (single quotes might work) so:
          Code:
          -outfmt "6 qcovs"
          If this does not work, you might want to think about calculating the coverage % with a script, taking query lenght from your BLASTp query and alignment length from 'length' of BLASTp output.

          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, Yesterday, 06:54 AM
          0 responses
          10 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-14-2024, 07:24 AM
          0 responses
          18 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-13-2024, 08:58 AM
          0 responses
          16 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 06-12-2024, 02:20 PM
          0 responses
          17 views
          0 likes
          Last Post seqadmin  
          Working...
          X