gapfillerRun: A function to confirm fusion break point by de novo assembly

Description Usage Arguments Value Author(s) See Also Examples

View source: R/denovoFusionCheck.R

Description

A function that uses GapFiller to confirm, by de novo assembly, the presence of the fusion break point. The function needs as input the fusion transcript generated by chimeraSeqs function and two fastq files corresponding to the reads mapping over the fusion transcript.

Usage

1
2
3
gapfillerRun(fusion.fa, seed1,  seed2, gapfiller=NULL, seed.ins=200, seed.var=50, 
block.length=5, overlap=20, max.length=5000, slack=30, 
k=6, global.mismatch=5)

Arguments

fusion.fa

fasta file with the fusin transcript

seed1

The R1 fastq of a pair-end

seed2

The R2 fastq of a pair-end

gapfiller

path to GapFiller executable program

seed.ins

seed reads insert size

seed.var

seed reads insert variation

block.length

length of perfect match

overlap

minimum suffix-prefix overlap

max.length

print only contigs <= max-length long

slack

number of overlaps: suffix-prefix overlap range is given by overlap, overlap + slack

k

length of the word used to hash

global.mismatch

maximum number of mismatches between mate and contig

Value

The program will write in a temporary directory contigs.fasta and contig.stats, which are used to evaluate if the de novo assembly allows the identification of the fusion break point. The function returns a list of three objects. The list is returned only in case that some of de novo assemblies cover the breakpoint junction. The list is made of:

contigs

which is a PairwiseAlignments object

junction.contigs

which is a DNAStringSet encompassing the sequences present in the contigs object

fusion

which is a DNAStringSet object encompassing the fusion transcript

Author(s)

Raffaele A Calogero

See Also

chimeraSeqs, gapfillerInstallation

Examples

1
#	tmp <- gapfillerRun(fusion.fa=paste(path.package("chimera", quiet = FALSE),"/examples/uc002xvp.1-243_uc002iyu.4-1031.fa",sep=""),seed1=paste(path.package("chimera", quiet = FALSE),"/examples/uc002xvp.1-243_uc002iyu.4-1031_R1.fastq",sep=""),  seed2=paste(path.package("chimera", quiet = FALSE),"/examples/uc002xvp.1-243_uc002iyu.4-1031_R1.fastq",sep=""))

Bioconductor-mirror/chimera documentation built on June 1, 2017, 5:30 a.m.