blast.vs.self.sh

Description

Determine BLAST+ hits between pairs of different sequences of the same FASTA file, and retain only one out of two reciprocal hits.

Usage

blast.vs.self.sh <file>

Dependencies

blast.fasta.seqs.sh
filter.blast.vs.self.R

Arguments

# Required
1)  Path to FASTA file.

# Optional
--

Details

A BLAST+ database is created from the FASTA file, and the same sequences are used as query and target in a BLAST+ search with expect value (evalue) 1E-04.

The BLAST+ results are filtered to retain only hits between pairs of DIFFERENT sequences, and only ONE out of two reciprocal hits.

Value

A *.vs.self.blast file containing the BLAST+ results. A *.vs.self.blast.filtered file containing the filtered BLAST+ results.

Examples

blast.vs.self.sh consensus.fasta