chimera (version 1.14.0)

chimeraSeqs: A function to generate the nucleotide sequences of a fusion event

Description

A function generating the nucleotide sequences of a chimera.

Usage

chimeraSeqs(fset, extend=1000, type="transcripts")

Arguments

fset

A fSet object.

extend

number of nucleotides used to extend a genomic region that is not an annotated gene. Default is 1000 nts

type

Chimera can be build at transcript level, i.e. transcript level will generate multiple fusion transcripts depending of the number of transcripts associated to each of the two genes involved in the fusion.

Value

A DNAStringSet encompassing the fusions generated using all the isoforms for each gene involved in the fusion.The name of each element of the DNAStringSet has the following format: gene1-lengthOfGeneFragment:gene2-lengthOfGeneFragment. In case the fusion junction is located in an intronic sequence, a warning is provided. The presence of a partial intron in the fusion is an indication that the fusion does not generate an active chimeric peptide.

See Also

fusionName, chimeraSeqSet

Examples

	tmp <- importFusionData("fusionmap", paste(find.package(package="chimera"),"/examples/mcf7.FMFusionReport", sep=""), org="hg19")
	fusion.names <- fusionName(tmp)
	fusion.names
	myset <- tmp[[13]]
	tmp.seq <- chimeraSeqs(myset, type="transcripts")
	#writeXStringSet(tmp.seq, paste(sub(":","_",fusion.names[[13]]),".fa",sep=""), format="fasta")