Description Usage Arguments Value Author(s) Examples
Note that samtools idxstats may or may not include secondary alignments. This will depend on how the reads were aligned. See the documention for samtools idxstats and the aligner used to create the bam file. The bam file must have been aligned to a reference genome containing chromosone GL000220.1 total reads mapping to GL0000220.1 or know rRNA regions defined by Ensembl
1 | GetrRNAStatistics(BamFilePath)
|
BamFilePath |
location of aligned bam file. Corresponding bai index file must be in the same directory. |
An integer for the total number of rRNA read counts in a given bam file.
Nathaniel J. Madrid, Jason Byars
1 2 3 4 5 | BamFilePaths <- dir(path="/home/ubuntu/Projects/RNAseqPackage/bams", pattern=".*.bam$", full.names=T)
# samtools idxstats must be installed on the local system.
# Note that samtools idxstats may or may not include secondary alignments. This will depend on how the reads
# were aligned. See the documention for samtools idxstats and the aligner used to create the bam file.
rRNAStatsDF = GetrRNAStatistics(BamFilePaths[1])
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.