[Bioperl-l] Reciprocal best hits using Bioperl?

Jay Hannah jay at jays.net
Mon Jan 18 18:36:20 EST 2010


On Jan 17, 2010, at 10:02 AM, Bhakti Dwivedi wrote:
> Is there a Bio-perl module to parse the reciprocal best hits (query1-> hit1
> && hit1 -> query1)  from a blast table report?

If all the advice and resources in this thread have not dissuaded you from writing your own, you could glance at cross_blast() here as reference:

   https://clabsvn.ist.unomaha.edu/anonsvn/user/jhannah/UNO/seqlab/seqlab/tutorial.pod

About the (abandoned) project:

   http://clab.ist.unomaha.edu/CLAB/index.php/SeqLab_%28Perl%29

I wrote that in 2006 for clustering a few hundred proteins based on custom criteria.

Cheers,

Jay Hannah
http://clab.ist.unomaha.edu/CLAB/index.php/User:Jhannah




More information about the Bioperl-l mailing list