quantify reads in consecutive bins
1 2 3 4 | countReads(samplenames, gtf, fragmentLength = 150, bamFolder,
outputDir = NA, modification = "m6A", binSize = 50,
strandToKeep = "opposite", paired = FALSE, threads = 1,
saveOutput = FALSE)
|
samplenames |
The names of each sample (prefix for bam files) |
gtf |
The gtf format gene annotation file |
fragmentLength |
The RNA fragment length (insert size of the library). |
bamFolder |
Path to the folder where bam file locates |
outputDir |
The directory to save output files |
modification |
The middle file name used to differentiate IP sample from Input sample. This should correspond to the file name of BAM files. For example, modification="m6A" for BAM files named as "samplenames.m6A.bam". |
binSize |
The size of consecutive bins to slice the transcripts |
strandToKeep |
According to library preparation protocol, choose which strand to count. Stranded RNA library usually seq the "ooposite" strand. Small RNA library seq the "same" strand. |
paired |
Logical indicating whether the input bam files are from paired end sequencing. Default is FALSE. If using paired end data, the read length will be estimated from the data and only good mate are counted. |
threads |
The number of threads to use for hyperthreading |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.