-
Notifications
You must be signed in to change notification settings - Fork 29
Description
I have the next bedMethyl file
chr1 10468 10469 m 4 . 10468 10469 255,0,0 4 75.00 3 1 0 1 2 1 1
chr1 10470 10471 m 6 . 10470 10471 255,0,0 6 83.33 5 1 0 1 0 1 1
chr1 10483 10484 m 6 . 10483 10484 255,0,0 6 66.67 4 2 0 2 1 0 0
chr1 10488 10489 m 7 . 10488 10489 255,0,0 7 71.43 5 2 0 1 1 0 0
chr1 10492 10493 m 7 . 10492 10493 255,0,0 7 100.00 7 0 0 1 1 0 0
chr1 10496 10497 m 6 . 10496 10497 255,0,0 6 100.00 6 0 0 1 2 0 0
chr1 10524 10525 m 6 . 10524 10525 255,0,0 6 83.33 5 1 0 0 3 0 0
chr1 10541 10542 m 4 . 10541 10542 255,0,0 4 100.00 4 0 0 0 5 0 0
chr1 10562 10563 m 8 . 10562 10563 255,0,0 8 87.50 7 1 0 0 1 0 0
chr1 10570 10571 m 8 . 10570 10571 255,0,0 8 100.00 8 0 0 0 1 0 0
When I try to parse it using bsseq::read.bedMethyl, I`ve got the next error message
Error in .local(x, ...) : invalid strand levels in 'x': .
12.
if (length(invalid_levels) != 0L)
stop("invalid strand levels in 'x': ",
paste(invalid_levels, collapse=", ")) at strand-utils.R#33
11.
.local(x, ...)
10.
function(x, ...) standardGeneric("strand") at strand.R#5
9.
strand(strand)
8.
eval(jsub, SDenv, parent.frame())
7.
jval = eval(jsub, SDenv, parent.frame()) at data.table.R#1383
6.
x[,:=(strand, strand(strand))] at read.bedMethyl.R#145
5.
x[,:=(strand, strand(strand))] at read.bedMethyl.R#145
4.
dt <- .readbedMethylAsDT(file = file,
col_spec = "BSseq",
check = TRUE,
verbose = verbose) at read.bedMethyl.R#10
3.
loci_from_first_file <- .readbedMethylAsFWGRanges(
file = files[[1L]],
rmZeroCov = rmZeroCov,
strandCollapse = strandCollapse, ... at read.bedMethyl.R#77
2.
loci <- .constructFWGRangesFrombedMethylFiles(
files = files,
rmZeroCov = rmZeroCov,
strandCollapse = strandCollapse, ... at read.bedMethyl.R#452
1.
bsseq::read.bedMethyl(samplesheet$modkit[1], strandCollapse = TRUE)
Why point is incorrect strand, if specification of bedMethyl are told us otherwise. Can I kindly ask you to fix it, please?
bsseq version 1.45.1