discard_bs_noise_reads: Discard BS-Seq noisy reads

In andreaskapou/processHTS: Process HTS Files and Data

Description

discard_bs_noise_reads discards low coverage and (really) high reads from BS-Seq experiments. These reads can be thought as noise of the experiment.

Usage

discard_bs_noise_reads(bs_data, min_bs_cov = 2, max_bs_cov = 1000)

Arguments

bs_data	A GRanges object containing the BS-Seq data.
min_bs_cov	The minimum number of reads mapping to each CpG site.
max_bs_cov	The maximum number of reads mapping to each CpG site.

Value

The clean GRanges object without noisy observations

Author(s)

C.A.Kapourani [e-mail]

Examples

bs_data_noisy  <- rrbs_data
bs_data <- discard_bs_noise_reads(bs_data_noisy, 4, 14)

andreaskapou/processHTS documentation built on May 12, 2019, 3:33 a.m.