Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • kulandaisamy
    Member
    • Mar 2015
    • 10

    standlone blast

    how to get query coverage in blastp?
  • skbrimer
    Member
    • Mar 2014
    • 55

    #2
    What you want in the outfmt option is qcovs, well at least its what I think you want.

    Comment

    • kulandaisamy
      Member
      • Mar 2015
      • 10

      #3
      i am running blastp in commandline.i didn't get query coverage results.how to got it?help me?

      Comment

      • GenoMax
        Senior Member
        • Feb 2008
        • 7142

        #4
        This manual documents the BLAST (Basic Local Alignment Search Tool) command line applications developed at the National Center for Biotechnology Information (NCBI).


        Check out the section on custom output formats.

        Detailed explanation of command line output format options for blast: https://www.biostars.org/p/88944/

        Comment

        • skbrimer
          Member
          • Mar 2014
          • 55

          #5
          How about this. Can you provide some information on how you are trying to run the blast command line tool? Are you getting an error? If you can give any of us an idea of what you are doing and perhaps why it is not working for you then maybe someone can give you the answer you are looking for. If you do not feel that providing anything other than "please help me" I suggest you type "blast -help" in the command line and get to reading.

          Comment

          • kulandaisamy
            Member
            • Mar 2015
            • 10

            #6
            blastall -p blastp -i plas_3d7.fasta -d approved_target.fasta -m 8 > blas.tab

            my output file contaiing only Query id, Subject id, % identity, alignment length, mismatches, gap openings, q. start, q. end, s. start, s. end, e-value, bit score.
            i want query coverage results

            Comment

            • kulandaisamy
              Member
              • Mar 2015
              • 10

              #7
              how to use -qcovs flag in blastp?

              Comment

              • skbrimer
                Member
                • Mar 2014
                • 55

                #8
                Originally posted by kulandaisamy View Post
                blastall -p blastp -i plas_3d7.fasta -d approved_target.fasta -m 8 > blas.tab
                I'm not sure how you would run the search in the older version of BLAST, NCBI recommends you download the new BLAST+ here is a link to a discussion https://www.biostars.org/p/9480/

                I have only used the newer version so the line of code would look something like....

                blastp -db nr -query plas_3d7.fasta -query_loc approved_target.fasta -out blas.tab -outfmt "6 qgi sgi qcovs" -remote

                Or something to this effect, I'm not sure what the -d flag is on the blastall version of blast but hopefully this can get you pointed in the right direction.

                Comment

                • kulandaisamy
                  Member
                  • Mar 2015
                  • 10

                  #9
                  ok.i tried to install blast+.but i didn't get it

                  Comment

                  • GenoMax
                    Senior Member
                    • Feb 2008
                    • 7142

                    #10
                    What version of blast are you running?
                    Code:
                    $ blastp -version

                    Comment

                    • kulandaisamy
                      Member
                      • Mar 2015
                      • 10

                      #11
                      blastp: 2.2.28+
                      Package: blast 2.2.28, build Jun 3 2013 11:32:37

                      Comment

                      • GenoMax
                        Senior Member
                        • Feb 2008
                        • 7142

                        #12
                        You are running blast+. Not the latest (v. 2.2.30) but close enough. Use @skbrimer's suggestion for the command above.

                        Comment

                        • kulandaisamy
                          Member
                          • Mar 2015
                          • 10

                          #13
                          ok..thnk you

                          Comment

                          • GenoMax
                            Senior Member
                            • Feb 2008
                            • 7142

                            #14
                            Don't use the -remote option if you are using blast locally.

                            Comment

                            • kulandaisamy
                              Member
                              • Mar 2015
                              • 10

                              #15
                              i tried but i didn't get it.give me some example.

                              Comment

                              Latest Articles

                              Collapse

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by SEQadmin2, Today, 11:58 AM
                              0 responses
                              6 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 06-05-2026, 10:09 AM
                              0 responses
                              21 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 06-04-2026, 08:59 AM
                              0 responses
                              34 views
                              0 reactions
                              Last Post SEQadmin2  
                              Started by SEQadmin2, 06-02-2026, 12:03 PM
                              0 responses
                              55 views
                              0 reactions
                              Last Post SEQadmin2  
                              Working...