In Chebuu/fastalign:

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

library(fastalign)

data(thcasaa)

print(thcasaa@metadata)

summary(thcasaa)
head(names(thcasaa))

inactive <- 'sp|Q33DQ2|THCAI_CANSA Inactive tetrahydrocannabinolic acid synthase OS=Cannabis sativa OX=3483 PE=3 SV=1'

multiPairwise(
  pattern = thcasaa,
  subject = thcasaa[[inactive]]
)

Chebuu/fastalign documentation built on July 19, 2020, 8:39 p.m.