Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • BLASTX Issues

    I am attempting to run blastx from the command line, and it seems to be very slow (~5-6 minutes for 10 100bp sequences). I compared this to blastx on the NCBI website, which took only 30 seconds. This indicates that there is an issue with the command line usage of blastx, but I don't know where to start examining this issue.

    I made sure that the db used in the command line and on the ncbi website was the same (nr). I am using all default parameters.

    Any help would be appreciated.

    Thanks

  • #2
    That is an apples to oranges comparison.

    Hardware NCBI uses for Blast searches is going to be tremendously more powerful than anything you are using locally.

    What kind of hardware are you running your searches on? How much RAM are you allocating to your searches locally?

    Comment


    • #3
      The machine I run the searches on is a 48 core 2.7ghz with 264 Gb RAM, but not sure how much of that is allocated for the search (nor do I know how to specify allocation). I run the search with default parameters so num_threads is set to 1 by default. Even though, I don't think it should be that slow. Or am I wrong?

      Comment


      • #4
        If you have 48 cores then you would certainly want to up that num_threads option. You may need to experiment with that number to find the sweet spot for your hardware configuration. Start at the half-way point (and go up and down).

        If you have admin access on the machine you may want to make a RAMDISK and copy the database files there (http://askubuntu.com/questions/15286...ake-a-ram-disk) If you are not using linux then look for an appropriate set of instructions for your OS.

        Comment


        • #5
          Increasing the threads definitely helps. I guess I'll fiddle around with that. Thanks GenoMax!

          Comment


          • #6
            Hi ..
            I am trying to get only the "no hits" result after perform the Blastx on the command line

            I've tried:

            <my Blastx command> | tee redirect.result33.txt | cut -d ' ' -f 1 | uniq | comm -31 - ids33_query.txt

            Any help would be appreciated.

            Thanks
            Last edited by debem; 12-09-2014, 05:42 PM.

            Comment


            • #7
              Hi debem,

              Input your results in to two files.csv and try:
              sdiff 1.csv 2.csv > difer.csv

              Open the file difer.csv and select the lines with "<".
              This is your differents IDs.

              Big hub,
              Rhewter.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Recent Developments in Metagenomics
                by seqadmin





                Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                09-23-2024, 06:35 AM
              • seqadmin
                Understanding Genetic Influence on Infectious Disease
                by seqadmin




                During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.

                Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...
                09-09-2024, 10:59 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 04:51 AM
              0 responses
              8 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 10-01-2024, 07:10 AM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-30-2024, 08:33 AM
              0 responses
              18 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 09-26-2024, 12:57 PM
              0 responses
              16 views
              0 likes
              Last Post seqadmin  
              Working...
              X