Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Reciprocal blast (RBH) - help

    Hi!

    I´m new in the bioinformatic world so be patience if my question is dummy.

    I have two fasta files and I run blast reciprocal best hits (RBH) for blastp using the script from peterjc that you can be found here: https://github.com/peterjc/galaxy_bl...ools/blast_rbh

    So, the output gives me:A_id, B_id, A_length, B_length, A_qcovhsp, B_qcovhsp, length, pident and bitscore.

    But I´m also interested in other options as e-values etc...

    Anyone knows if there is an away of change the script that I´m using to obtain other results or if there is another script that can give these options ?

    thanks!!

  • #2
    I'm glad it's useful (our preprint covering this and the Galaxy BLAST+ wrappers is out, and the formal publication in GigaScience is coming soon - so you'll be able to formally cite the script).

    The e-values depend on the size of the FASTA files (or the temporary databases made from them), so are often not directly useful for comparison between datasets - that's why I deliberately only reported the bit score.

    But yes, it ought to be fairly simple for you to add the evalue field to the BLAST commands and collect that in the final output. Do you know any Python programming?

    Comment


    • #3
      First of all, thanks for your fast answer.

      I had the idea that e-values were the main criteria to "eliminate" some hits. I´m wrong?

      Yes, I have some basic knowledge in python programing.

      Comment


      • #4
        Maybe read http://dx.doi.org/10.1371/journal.pcbi.1000160 (linked to from the tool help in the Galaxy wrapper XML file).

        My RBH tool does not use the e-value at all, the "best" hit has the highest bit score (which in turn means it has the lowest e-value because BLAST calculates that from the bit score).

        Comment


        • #5
          Thank you very much!!!

          Comment

          Latest Articles

          Collapse

          • seqadmin
            The Impact of AI in Genomic Medicine
            by seqadmin



            Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
            02-26-2024, 02:07 PM
          • seqadmin
            Multiomics Techniques Advancing Disease Research
            by seqadmin


            New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

            A major leap in the field has
            ...
            02-08-2024, 06:33 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 02-23-2024, 04:11 PM
          0 responses
          57 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-21-2024, 08:52 AM
          0 responses
          67 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-20-2024, 08:57 AM
          0 responses
          56 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 02-14-2024, 09:19 AM
          0 responses
          65 views
          0 likes
          Last Post seqadmin  
          Working...
          X