addCoverage: Import sequencing data

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/qsea.createSet.R

Description

This function imports the alignment files (in sam/bam format) and counts the reads per genomic window or directly imports coverage files (in wiggle/bigwiggle format)

Usage

1
2
addCoverage(qs, fragment_length, uniquePos=TRUE, minMapQual=1, paired=FALSE, 
parallel=FALSE)

Arguments

qs

qseaSet object, e.g. produced by the createQseaSet() function

fragment_length

For single end data, provide the expected fragment length

paired

If set to TRUE, data is considered to be paired end sequencing, and the actual fragments size is used.

uniquePos

If set to TRUE, fragments with same position and orientation are considered to be PCR duplicates and replaced by one representative.

minMapQual

The minimal mapping quality for reads to be considered. Note that the definition of mapping quality depends on the alignment tool.

parallel

Switch for parallel computing, using BiocParallel

Details

The coverage is imported from the files specified in the file_name column of the sample table, provided for the createQseaSet() function. In case of alignment files, the reads are counted for the window at the center of the sequencing fragment. For single end data, Filetypes is detected automatically from the file suffix.

Value

The function returns the qseaSet object, extended by the number of reads per window for all samples

Author(s)

Mathias Lienhard

See Also

crateQseaSet

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
library("BSgenome.Hsapiens.UCSC.hg19")

bam_hESCs_1 = system.file("extdata", 
    "hESCs.MeDIP.Rep1.chr22.bam", package="MEDIPSData")
bam_hESCs_2 = system.file("extdata", 
    "hESCs.MeDIP.Rep2.chr22.bam", package="MEDIPSData")
sample_table=data.frame(sample_name=paste0("hESCs_", 1:2), 
    file_name=c(bam_hESCs_1,bam_hESCs_2), 
    group=rep("hESC",2), stringsAsFactors=FALSE)
qseaSet=createQseaSet(sampleTable=sample_table, 
        BSgenome="BSgenome.Hsapiens.UCSC.hg19", 
        chr.select="chr22", 
        window_size=500)
qseaSet=addCoverage(qseaSet, fragment_length=300)

qsea documentation built on Nov. 8, 2020, 8:28 p.m.