extractSequencesR {VDJgermlines}

Extracts the germline genes of interest and stores them as a dataframe

Description

Extracts the germline genes of interest and stores them as a dataframe

Usage

extractSequencesR(species, chain, source, write.fasta, fasta.name,
  directory.string)

Arguments

species

Specifies desired species. Current species options can be found by running unique(VDJgermlines$species) after loading in data(VDJgermlines). "all_species" will select all available species. Species include ("alpaca","bovine","crab_eating_macaque","dog", "human","mouse","pig","platypus","rabbit",raintrout","rat","rhesus monkey","salmon", "sheep","zebrafish","atlanticcod","catfish","teleostei","goat","dolphin","mamonkey", "camel","runny")

chain

Specifies the desired chain. "all_chains" will select all chains. The various chains can be listed by running unique(VDJgermlines$chain).Antibody loci start with IG, TCR starts with TR.NExt comes either heavy, light, alpha, beta etc. Examples include IGHD IGHJ IGHV IGIJ IGIV IGKJ IGKV IGLJ IGLV TRAJ TRAV TRBD TRBJ TRBV TRDD TRDJ TRDV TRGJ TRGV

source

Specifies from which source the desired chains should be taken from. Current sources can be found by running unique(VDJgermlines$source) after loading in data(VDJgermlines). IMGT and vgenerepertoire are currently present.

write.fasta

A boolean (TRUE / FALSE) value that species if the output sequences should be also written as a fasta file.

fasta.name

A string specifying the name of the fasta file. If this is left blank this will be defaulted to "VDJgermlines.fasta". Manually including extension is not needed.

directory.string

A string stating the directory that the fasta file should be saved to. If left blank the fasta file will be stored in the current working directory.

Value

Returns a dataframe containing the species and chains of interest from the desired source. The columns of the dataframe include sequence, names, species, chain, accession, official species name, source, and strain. Furthermore, if write.fasta is set to TRUE then the output will additionally include a fasta file.

Examples

extractSequencesR(species = "dog",chain = "IGHD",source = "IMGT",write.fasta = FALSE)