Description Usage Arguments Value Author(s)
View source: R/createCompositeFile.R
This function will move through BAM files in a folder, read in each individual file and go through each chromosome, determine if the chromosome is WW or CC based on WCcutoff, reverse complement all reads in the WW file, append to a new composite file for that chromosome, order the composite file of each chromosome based on position.
1 2 3 4 5 6 7 8 9 10 11 | createCompositeFile(
file.list,
chromosomes = NULL,
pairedEndReads = TRUE,
pair2frgm = FALSE,
min.mapq = 10,
filtAlt = FALSE,
WC.cutoff = 0.9,
genoT = "fisher",
background = 0.05
)
|
file.list |
A list of BAM files to process. |
chromosomes |
If only a subset of the chromosomes should be binned, specify them here. |
pairedEndReads |
Set to |
pair2frgm |
Set to |
min.mapq |
Minimum mapping quality when importing from BAM files. |
filtAlt |
Set to |
WC.cutoff |
Percentage of WW or CC reads to consider chromosome being WW or CC |
genoT |
A method ('fisher' or 'binom') to genotype regions defined by a set of breakpoints. |
background |
The percent (e.g. 0.05 = 5%) of background reads allowed for WW or CC genotype calls. |
A GRanges-class
object.
Ashley Sanders, David Porubsky
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.