Using the table returned from createIPRScanDomainTable, construct a domain fasta for a single accession number in the original fasta (i.e., the original fasta argument to createIPRScanDomainTable())
Source:R/fa2domain.R
convertIPRScanDomainTable2FA.Rd
Using the table returned from createIPRScanDomainTable, construct a domain fasta for a single accession number in the original fasta (i.e., the original fasta argument to createIPRScanDomainTable())
Value
AAStringSet A domain fasta containing all the domains for a single protein in the original fasta passed as an argument to createIPRScanDomainTable()
Examples
if (FALSE) { # \dontrun{
path_molevol_scripts <- file.path(Sys.getenv("DEV", unset = "/data/molevolvr_transfer/molevolvr_dev"), "molevol_scripts")
setwd(path_molevol_scripts)
source("R/fa2domain.R")
fasta <- Biostrings::readAAStringSet("./tests/example_protein.fa")
df_iprscan <- readIPRScanTSV("./tests/example_iprscan_valid.tsv")
accnum <- df_iprscan$AccNum[1]
df_iprscan_domains <- createIPRScanDomainTable(accnum, fasta, df_iprscan)
fasta_domains <- df_iprscan_domains |> convertIPRScanDomainTable2FA()
} # }