createCompositeFile: Function to create composite Strand-seq file

Description Usage Arguments Author(s)

Description

This script will move through .bam files in a folder read in each individual file and go thru each chr determine the chr is WW or CC based on WCcutoff reverse compliment all reads in the WW file append to a new composite file for that chr order the composite file of each chr based on pos

Usage

1
2
createCompositeFile(file.list, chromosomes = NULL, pairedEndReads = TRUE,
  min.mapq = 10, keep.duplicate.reads = FALSE, WC.cutoff = 0.9)

Arguments

chromosomes

If only a subset of the chromosomes should be binned, specify them here.

pairedEndReads

Set to TRUE if you have paired-end reads in your file.

min.mapq

Minimum mapping quality when importing from BAM files.

keep.duplicate.reads

A logical indicating whether or not duplicate reads should be kept.

WC.cutoff

Percentage of WW or CC reads to consider chromosome being WW or CC

datapath

Location where bam files to be processed are stored

Author(s)

Ashley Sanders, David Porubsky


drashley/InversionAnalysis_HGSVC documentation built on May 6, 2019, 8:49 a.m.