Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • FASTA database of human_g1k_v37.fasta

    Hi all,
    I want to BLAST using the human_g1k_v37_fasta as database. When i run formatdb on this file and run blast subsequent to that, i get a message stating index file not found.

    Can someone please let me know if it's possible to convert this file to FASTA database and if so the command to be used.

    Thanks in advance.

    Ashwin

  • #2
    In principle there is no reason why human_g1k_v37_fasta cant be used as a BLAST database.

    I suggest reading this: Ten Simple Rules for Getting Help from Online Scientific Communities


    For a start, which version of BLAST are you using, which OS are you using, and what were the commands you ran to build the database, and to execute the BLAST search, and what was the EXACT error message?

    Comment


    • #3
      Follow up info

      Thanks for link to that article.

      I am using blast-2.2.25+ on Ubuntu and formatdb -n human_g1k_v37_fasta.fasta

      The exact error i get when i try to run blastn is index file not found for database human_g1k_v37_fasta

      Thanks.

      Ashwin

      Comment


      • #4
        What files did formatdb create? And why are you using formatdb - that's from legacy BLAST, you should be using makeblastdb for BLAST+

        Comment


        • #5
          Files and still getting error

          Right, so i used makeblastdb as follows

          makeblastdb -in human_g1k_v37.fasta -title human_g1k_v37

          it made three files each with
          .psq, .phr and .pin file extensions

          i am trying to execute blastn as follows

          blasntn -query {query file} -db human_g1k_v37 -out Results.txt -strand plus

          I am getting error 'No index or alias found for database human_g1k_v37'

          thanks

          Comment


          • #6
            Originally posted by ashkot View Post
            Right, so i used makeblastdb as follows

            makeblastdb -in human_g1k_v37.fasta -title human_g1k_v37

            it made three files each with
            .psq, .phr and .pin file extensions
            From the extensions that made a protein database (which is the default behaviour IIRC). There is an option to makeblastdb to make a nucleotide database which you'll need to use.

            Comment


            • #7
              Update

              Thanks, a lot, BLAST is working now. It had to do with the database being not available in the right format.

              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, 07:49 AM
              0 responses
              13 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-20-2024, 07:23 AM
              0 responses
              14 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-17-2024, 06:54 AM
              0 responses
              16 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 06-14-2024, 07:24 AM
              0 responses
              25 views
              0 likes
              Last Post seqadmin  
              Working...
              X