View source: R/synchronizeReadDirRegion.R
synchronizeReadDirRegion | R Documentation |
This function aims to synchronize read directionality of a specified genomic region.
synchronizeReadDirRegion(
bamfiles = "",
region = NULL,
pairedEndReads = TRUE,
min.mapq = 10,
filterAltAlign = TRUE,
alpha = 0.1,
genot.frac.reads = 0.25,
genot.region.ends = "right",
desired.genot = c("cc", "ww")
)
bamfiles |
A list of BAM files to be processed. |
region |
A |
pairedEndReads |
Set to |
min.mapq |
Minimum mapping quality when importing from BAM files. |
filterAltAlign |
Set to |
alpha |
Estimated level of background in Strand-seq reads. |
genot.frac.reads |
A fraction of the total number of reads to be genotype from both or a single end of a given 'region' based on 'genot.region.ends'. |
genot.region.ends |
Define if 'left' or 'right' of 'both ends of a given region shall be genotyped. |
desired.genot |
Select desired genotypes: 'ww', 'cc', 'wc' to be allowed on selected region ends. If 'genot.region.ends' = 'both', both ends have to have the same genotype. |
A GRanges-class
object that reads synchronized by directionality.
David Porubsky
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.