prepare_blast_pairs: Convert BLAST+ alignment table to two outputs. First, a full...

In brandonsie/phipmake: PhIP-seq Data Analysis Pipeline

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.

Description

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.

Usage

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

brandonsie/phipmake documentation built on March 15, 2023, 3:24 p.m.