git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/gmapR@68552 bc3139a8-67e5-0310-9ffc-ced21a209358
... | ... |
@@ -1,3 +1,4 @@ |
1 |
+ |
|
1 | 2 |
% \VignetteIndexEntry{gmapR} |
2 | 3 |
% \VignetteKeywords{gmap,gsnap} |
3 | 4 |
% \VignettePackage{gmapR} |
... | ... |
@@ -256,9 +257,9 @@ bam_file <- file.path(system.file(package = "gmapR", mustWork=TRUE), |
256 | 257 |
breaks <- c(0L, 15L, 60L, 75L) |
257 | 258 |
bqual <- 56L |
258 | 259 |
mapq <- 13L |
259 |
-param <- BamTallyParam(cycle_breaks = as.integer(breaks), |
|
260 |
- high_quality_cutoff = as.integer(bqual), |
|
261 |
- minimum_mapq = as.integer(mapq), |
|
260 |
+param <- BamTallyParam(genome, cycle_breaks = breaks, |
|
261 |
+ high_quality_cutoff = bqual, |
|
262 |
+ minimum_mapq = mapq, |
|
262 | 263 |
concordant_only = FALSE, unique_only = FALSE, |
263 | 264 |
primary_only = FALSE, |
264 | 265 |
min_depth = 0L, variant_strand = 1L, |
... | ... |
@@ -266,7 +267,6 @@ param <- BamTallyParam(cycle_breaks = as.integer(breaks), |
266 | 267 |
indels = FALSE) |
267 | 268 |
|
268 | 269 |
tallies <-bam_tally(bam_file, |
269 |
- genome, |
|
270 | 270 |
param) |
271 | 271 |
@ |
272 | 272 |
|