Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • mknut
    replied
    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.

    Leave a comment:


  • GenoMax
    replied
    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.

    Leave a comment:


  • Tsuyoshi
    replied
    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.

    Leave a comment:


  • mknut
    replied
    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.

    Leave a comment:


  • Tsuyoshi
    started a topic Basic question about local BLASTP

    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.

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