Convert BLAST+ alignment table to two outputs. First, a full symmetrical two column table of aligning peptides. Second, a subset of the first, including only intra-protein alignments, for which protein identity is specified by an annotation file.

prepare_blast_pairs(blast.file, annot, name1, name2)

Arguments

blast.file

directory path to blast file e.g. 1/PhageomeLa_000_BLAST.txt

annot

Annotation file.

name1

file path & name for all-alignment output file tsv

name2

file path & name for intraprotein-alignment output file tsv