readsToFragments {Basic4Cseq} | R Documentation |
This function maps aligned reads to fragment ends of the virtual fragment library to calculate the coverage of the fragments. The number of reads at the start and end of a fragment is provided, as well as the average of both fragment ends.
readsToFragments(expData, fragmentLib)
expData |
Experiment data of class |
fragmentLib |
Fragment library for the given genome and cutting enzyme combination |
Data frame containing fragment-based data, i.e. a fragment's position and read coverage
Carolin Walter
data(liverData) file <- system.file("extdata", "vfl_aagctt_catg_mm9_54_vp.csv", package="Basic4Cseq") rawFragments(liverData) = readsToFragments(liverData, file) head(rawFragments(liverData))