Skip to main content
Join the official 2020 Python Developers SurveyStart the survey!

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

Project description

# Blast-Besties

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

# Table of contents

* [Example usage](#example-usage)
* [Input requirements](#input-requirements)
* [Standard options](#standard-options)
* [License](#license)

# Example usage

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

blastp -qcov_hsp_perc 95 -query A_prot.fa -subject B_prot.fa -out -evalue 0.001 -outfmt 6 -use_sw_tback
blastp -qcov_hsp_perc 95 -query B_prot.fa -subject A_prot.fa -out -evalue 0.001 -outfmt 6 -use_sw_tback

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

./ -e 0.001 -l 100 -s 100 -a ** -b *

## 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: BLAST-Besties [-h] [-v] -a BLASTAVB -b BLASTBVA [-l MINLEN] [-e EVAL]


# 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
-b, --blastBvA Blast tab result file for fastaB query against fastaA

# 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.

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for blastbesties, version 1.1.1
Filename, size File type Python version Upload date Hashes
Filename, size blastbesties-1.1.1.tar.gz (5.0 kB) File type Source Python version None Upload date Hashes View

Supported by

Pingdom Pingdom Monitoring Google Google Object Storage and Download Analytics Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page