Description Usage Arguments Author(s)
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
1 2 | createCompositeFile(file.list, chromosomes = NULL, pairedEndReads = TRUE,
min.mapq = 10, keep.duplicate.reads = FALSE, WC.cutoff = 0.9)
|
chromosomes |
If only a subset of the chromosomes should be binned, specify them here. |
pairedEndReads |
Set to |
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 |
Ashley Sanders, David Porubsky
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.