Displaying 1 result from an estimated 1 matches for "mapq".
Did you mean:
map
2012 Dec 17
1
Code works standalone, yet same code fails when part of package
...an error.
Please help!
Mick
> library(Rsamtools)
> source("viRome/R/read.bam.R")
> read.bam
function(bamfile=NULL, chr=NULL, start=1, end=1e07, what=c("qname", "flag",
"rname", "strand", "pos", "qwidth", "mapq", "cigar", "mrnm", "mpos",
"isize", "seq"), tag=c("NM"), removeN=TRUE) {
which <- RangesList(chr = IRanges(start,end))
names(which) <- chr
param <- ScanBamParam(which = which, what = what, tag...