blastbesties 1.1.2

Rapid discovery of reciprocal best blast pairs from BLAST output files.

Blast-Besties

Rapid discovery of reciprocal best blast pairs from BLAST output files.

Table of contents

• Blast-Besties
  • Table of contents
  • Install Blast-Besties
  • Example usage
    • Input requirements
    • Standard options
  • License

Install Blast-Besties

There are 3 options available for installing Blast-Besties locally:

1. Clone from this repository and install as a local Python package.

git clone https://github.com/Adamtaranto/blast-besties.git && cd blast-besties && pip install -e .

1. Pip install directly from this git repository.

pip install git+https://github.com/Adamtaranto/blast-besties.git

3. Install from PyPi.

pip install blastbesties

Example usage

Run BLASTp for each protein set AvsB and BvsA. Require vaild alignments to cover 90% of the query sequence and with an e-value of < 0.001.

blastp -qcov_hsp_perc 90 -query A_prot.fa -subject B_prot.fa -out AvB.tab -evalue 0.001 -outfmt 6 -use_sw_tback
blastp -qcov_hsp_perc 90 -query B_prot.fa -subject A_prot.fa -out BvA.tab -evalue 0.001 -outfmt 6 -use_sw_tback

Report reciprocal best match pairs where each hit meets criteria e-value <= 0.001, hit-length >= 40, bitscore >= 100.

blastbesties -e 0.001 -l 40 -s 100 -a AvB.tab -b BvA.tab -o pairs.tab

Input requirements

• Takes two BLAST output files as input. SetA vs SetB, and SetB vs SetA.
• BLAST output format 6 (tabular)
• Hits must be sorted by query name then descending match quality - this is default BLAST behaviour.

Standard options

Usage: blastbesties [-h] [-v] -a BLASTAVB -b BLASTBVA [-l MINLEN] [-e EVAL]
                     [-s BITSCORE] [-o OUTFILE] [-d OUTDIR]

Options:

  # Info
  -h, --help        Show this help message and exit.
  -v, --version     Show program's version number and exit.
  
  # Input 
  -a, --blastAvB    Blast tab result file for fastaA query against fastaB
                    subject.
  -b, --blastBvA    Blast tab result file for fastaB query against fastaA
                    subject.
  
  # Settings
  -l, --minLen      Minimum length of hit to consider valid.
  -e, --eVal        Maximum e-value to consider valid pair.
  -s, --bitScore    Minimum bitscore to consider valid pair.

  # Output
  -o, --outFile     Write reciprocal blast pairs to this file.
  -d, --outDir      Directory for new sequence files to be written to.

License

Software provided under MIT license.