Entering edit mode
Jeremy Leipzig
▴
70
@jeremy-leipzig-4924
Last seen 6.1 years ago
BiocInstaller::biocLite(("RNAseqData.HNRNPC.bam.chr14")) BiocInstaller::biocLite("Rsamtools") library(Rsamtools) library("RNAseqData.HNRNPC.bam.chr14") fl <- RNAseqData.HNRNPC.bam.chr14_BAMFILES[1] res <- pileup(fl) head(res) > res <- pileup(fl) > > head(res) [1] seqnames pos strand nucleotide count <0 rows> (or 0-length row.names) > table(res$strand, res$nucleotide) A C G T N = - + + 0 0 0 0 0 0 0 0 - 0 0 0 0 0 0 0 0 * 0 0 0 0 0 0 0 0 > sessionInfo() R version 3.4.1 (2017-06-30) Platform: x86_64-redhat-linux-gnu (64-bit) Running under: Amazon Linux AMI 2017.09 Matrix products: default BLAS/LAPACK: /efs/R/lib/libRblas.so locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 LC_PAPER=en_US.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats4 parallel stats graphics grDevices utils datasets methods base other attached packages: [1] RNAseqData.HNRNPC.bam.chr14_0.12.0 Rsamtools_1.26.2 Biostrings_2.42.1 [4] XVector_0.14.1 GenomicRanges_1.26.4 GenomeInfoDb_1.10.3 [7] IRanges_2.8.2 S4Vectors_0.12.2 BiocGenerics_0.20.0 loaded via a namespace (and not attached): [1] pcaPP_1.9-72 Rcpp_0.12.11 DEoptimR_1.0-8 compiler_3.4.1 BiocInstaller_1.24.0 RColorBrewer_1.1-2 [7] zlibbioc_1.20.0 bitops_1.0-6 tools_3.4.1 dotCall64_0.9-5 lattice_0.20-35 graph_1.52.0 [13] mvtnorm_1.0-6 spam_2.1-2 hexbin_1.27.1 cluster_2.0.6 IDPmisc_1.1.17 fields_9.0 [19] maps_3.2.0 grid_3.4.1 robustbase_0.92-8 Biobase_2.34.0 rrcov_1.4-3 R6_2.2.2 [25] BiocParallel_1.8.2 latticeExtra_0.6-28 corpcor_1.6.9 matrixStats_0.52.2 flowFramePlus_0.1.0 MASS_7.3-48 [31] assertthat_0.2.0 flowCore_1.40.6 KernSmooth_2.23-15 flowViz_1.38.0 RCurl_1.95-4.9
same after update to Rsamtools_1.30.0 and RNAseqData.HNRNPC.bam.chr14_0.16.0