bbmix (version 1.0.0)

poolreads: Pool randomly selected reads from different files

Description

Pool randomly selected reads from different files

Usage

poolreads(files, N = 1000, d = 10, out)

Value

save files

Arguments

files

names for files to extract reads

N

number of reads to extract

d

depth for reads

out

file name to save reads

Examples

counts_f <- system.file("extdata/input", "NA12878.chr22.Q20.allelicCounts.txt",
package = "bbmix",
mustWork = TRUE)

## In this example we only use one file and we take a pool of 10 reads

out <- tempfile()

poolreads(files=counts_f,
N=10,
d=10,
out = out)

unlink(out)