Skip to contents

Convert an alignment position to residue position

Source: R/poi_to_res.R
poi_to_res.Rd

This function converts an alignment position to a position in the frame of the reference protein sequence, i.e., to the positions of the amino acids in the first sequence of the alignment.

Usage

poi_to_res(alignment, poi)

Arguments

alignment

An alignment.

poi

An alignment position.

Value

An integer vector of positions of the amino acid residues in the reference sequence.

Examples

align_ATM <- read_alignment('ATM')
align_ATM[, 1:5]
#> Hsap_ATM_AAB65827.1  1 M-SLV
#> Mmus_ATM_NP_031525.2 1 M-SLA
#> Sscr_ATM_AAT01608.1  1 M-SLA
#> Mdom_ATM_IARC        1 M-SLA
#> Ggal_ATM_edited      1 M-SLV
#> Xlae_ATM_AAT72929.1  1 M-SLA
#> Drer_ATM_IARC_v2     1 M-SLA
#> Bflo_ATM_IARC        1 MTDLL
#> Spur_ATM_ABY60856.1  1 MAEVL

# Convert the positions of the first five alignment positions to residue positions
poi_to_res(align_ATM, 1:5)
#> [1]  1 NA  2  3  4