/BackBLAST_Reciprocal_BLAST

This repository contains a reciprocal BLAST program for filtering down BLAST results to best bidirectional hits. It also contains a toolkit for finding and visualizing BLAST hits for gene clusters within multiple bacterial genomes.

Primary LanguagePythonMIT LicenseMIT

BackBLAST_Reciprocal_BLAST

This repository contains a reciprocal BLAST program for filtering down BLAST results to best bidirectional hits. It also contains a toolkit for finding and visualizing BLAST hits for gene clusters within multiple bacterial genomes.

  • BackBLAST.py - A script that uses NCBI BLAST to search for gene clusters within a within a bacterial genome genome. Non-orthalagous genes are filtred out by identifying and extracting only Bidirectional BLAST Hits using a graph-based algorithm. The algorithm is illustrated below:

BackBLAST Algorithm

  • Visualization - This repository also includes tools for visualizing the results from BackBLAST.py in the form of a R heatmap. Here is an example:

Example Results

For more thorough descriptions and information on usage please check the [wiki!] (https://github.com/LeeBergstrand/BackBLAST-Gene-Cluster-Finder/wiki)