SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
reciprocal blast renesh Bioinformatics 7 02-10-2017 08:33 AM
Reciprocal Blast Guidance... James88 Bioinformatics 4 11-05-2014 02:45 AM
Reciprocal blast help Shorash Bioinformatics 5 07-21-2014 10:39 PM
detect orthologs using reciprocal BLAST nicedad Bioinformatics 2 06-07-2012 07:37 AM
Need Help Regarding Reciprocal best hit in BLAST Bharat Bioinformatics 1 03-11-2010 03:04 AM

Reply
 
Thread Tools
Old 08-23-2015, 03:02 AM   #1
ac91
Member
 
Location: porto

Join Date: Aug 2015
Posts: 13
Default Reciprocal blast (RBH) - help

Hi!

Im 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 Im also interested in other options as e-values etc...

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

thanks!!
ac91 is offline   Reply With Quote
Old 08-23-2015, 04:19 AM   #2
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

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?
maubp is offline   Reply With Quote
Old 08-23-2015, 05:11 AM   #3
ac91
Member
 
Location: porto

Join Date: Aug 2015
Posts: 13
Default

First of all, thanks for your fast answer.

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

Yes, I have some basic knowledge in python programing.
ac91 is offline   Reply With Quote
Old 08-23-2015, 12:25 PM   #4
maubp
Peter (Biopython etc)
 
Location: Dundee, Scotland, UK

Join Date: Jul 2009
Posts: 1,543
Default

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).
maubp is offline   Reply With Quote
Old 08-23-2015, 12:34 PM   #5
ac91
Member
 
Location: porto

Join Date: Aug 2015
Posts: 13
Default

Thank you very much!!!
ac91 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:40 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO