Learn R Programming

Basic4Cseq (version 1.8.0)

readsToFragments: Determine fragment coverage of a 4C-seq fragment library

Description

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.

Usage

readsToFragments(expData, fragmentLib)

Arguments

expData
Experiment data of class Data4Cseq with information on the 4C-seq experiment, including raw 4C-seq read data
fragmentLib
Fragment library for the given genome and cutting enzyme combination

Value

Data frame containing fragment-based data, i.e. a fragment's position and read coverage

Examples

Run this code
  data(liverData)
  file <- system.file("extdata", "vfl_aagctt_catg_mm9_54_vp.csv", package="Basic4Cseq")
  rawFragments(liverData) = readsToFragments(liverData, file)
  head(rawFragments(liverData))

Run the code above in your browser using DataLab