GUIDEseq (version 1.2.1)

GUIDEseq-package: Analysis of GUIDE-seq

Description

The package includes functions to retain one read per unique molecular identifier (UMI), filter reads lacking integration oligo sequence, identify peak locations (cleavage sites) and heights, perform target and off target search of the input gRNA. This package leverages CRISPRseek and ChIPpeakAnno packages.

Details

Package:

GUIDEseq

Type:

Package

Version:

1.0

Date:

2015-09-04

License:

GPL (>= 2)

Function GUIDEseqAnalysis integrates all steps of GUIDE-seq analysis into one function call

References

Shengdar Q Tsai and J Keith Joung et al. GUIDE-seq enables genome-wide profiling of off-target cleavage by CRISPR-Cas nucleases. Nature Biotechnology 33, 187 to 197 (2015)

See Also

GUIDEseqAnalysis

Examples

    if(interactive())
    {
        umiFile <- system.file("extdata", "UMI-HEK293_site4_R1.txt", 
            package = "GUIDEseq")
        alignFile <- system.file("extdata","bowtie2.HEK293_site4.sort.bed" , 
            package = "GUIDEseq")
        gRNA.file <- system.file("extdata","gRNA.fa", package = "GUIDEseq")
        guideRes <- GUIDEseqAnalysis(
            alignment.inputfile = alignFile ,
            umi.inputfile = umiFile, gRNA.file = gRNA.file)
        names(cleavages)
   }