Instructions on loading BAM files and indices in MMAPPR2data package.
MMAPPR2data 1.16.0
We first load the MMAPPR2data
package:
library(MMAPPR2data, quietly = TRUE)
This package contains the following two functions, which
provide easy access to the BAM files and their indices,
returning BamFile
objects:
exampleWTbam()
## class: BamFile
## path: /tmp/RtmpUvKRNa/Rinst10c68d26c72ae/MMAPPR2data/extdata/wt.bam
## index: /tmp/RtmpUvKRNa/Rinst10c68d26c72ae/MMAPPR2data/extdata/wt.bam.bai
## isOpen: FALSE
## yieldSize: NA
## obeyQname: FALSE
## asMates: FALSE
## qnamePrefixEnd: NA
## qnameSuffixStart: NA
exampleMutBam()
## class: BamFile
## path: /tmp/RtmpUvKRNa/Rinst10c68d26c72ae/MMAPPR2data/extdata/mut.bam
## index: /tmp/RtmpUvKRNa/Rinst10c68d26c72ae/MMAPPR2data/extdata/mut.bam.bai
## isOpen: FALSE
## yieldSize: NA
## obeyQname: FALSE
## asMates: FALSE
## qnamePrefixEnd: NA
## qnameSuffixStart: NA
Annotation data for the region is also included with the package and can be accessed with these two functions:
goldenFasta()
## [1] "/tmp/RtmpUvKRNa/Rinst10c68d26c72ae/MMAPPR2data/extdata/slc24a5.fa.gz"
goldenGFF()
## [1] "/tmp/RtmpUvKRNa/Rinst10c68d26c72ae/MMAPPR2data/extdata/slc24a5.gff.gz"
For details on the source of these files, and on their construction
see ?MMAPPR2data
and the inst/scripts/
folder.
sessionInfo()
## R version 4.3.1 (2023-06-16)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.18-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] MMAPPR2data_1.16.0 BiocStyle_2.30.0
##
## loaded via a namespace (and not attached):
## [1] crayon_1.5.2 cli_3.6.1 knitr_1.44
## [4] rlang_1.1.1 xfun_0.40 jsonlite_1.8.7
## [7] S4Vectors_0.40.0 RCurl_1.98-1.12 Biostrings_2.70.1
## [10] htmltools_0.5.6.1 stats4_4.3.1 sass_0.4.7
## [13] rmarkdown_2.25 evaluate_0.22 jquerylib_0.1.4
## [16] bitops_1.0-7 fastmap_1.1.1 GenomeInfoDb_1.38.0
## [19] IRanges_2.36.0 yaml_2.3.7 bookdown_0.36
## [22] BiocManager_1.30.22 compiler_4.3.1 codetools_0.2-19
## [25] XVector_0.42.0 BiocParallel_1.36.0 digest_0.6.33
## [28] R6_2.5.1 GenomeInfoDbData_1.2.11 parallel_4.3.1
## [31] GenomicRanges_1.54.0 bslib_0.5.1 tools_4.3.1
## [34] Rsamtools_2.18.0 zlibbioc_1.48.0 BiocGenerics_0.48.0
## [37] cachem_1.0.8
Thanks to Mike Love’s alpineData package for vignette structure inspiration.