Error in unite() when specifying destrand = TRUE #349
Description
Hi,
As written in the title, it showed:
[E::hts_idx_push] Chromosome blocks not continuous
Error in value[[3L]](cond) : index build failed
We have double checked that the methylRawListDB object does contain strand information (i.e + and -). We have also traced back the building blocks of unite() function through: destrandFun() -> applyTbxByChunk() -> getTabixByChunk(). We guess that getTabixByChunk() is presumably the origin of error.
Do you have any thoughts on fixing the issues?
Example of our methylRawListDB object
> myobj[[1]] |> head(n=30)
chr start end strand coverage numCs numTs
1 chr1 29268 29268 - 12 0 12
2 chr1 29273 29273 - 12 1 11
3 chr1 29286 29286 - 12 0 12
4 chr1 29299 29299 - 11 0 11
5 chr1 29301 29301 - 11 0 11
6 chr1 29307 29307 - 10 0 10
7 chr1 29311 29311 - 13 1 12
8 chr1 29348 29348 - 14 0 14
9 chr1 29351 29351 - 13 0 13
10 chr1 29354 29354 - 14 0 14
11 chr1 29360 29360 - 12 0 12
12 chr1 29368 29368 - 14 0 14
13 chr1 29374 29374 - 14 2 12
14 chr1 29390 29390 - 13 1 12
15 chr1 29393 29393 - 12 0 12
16 chr1 29404 29404 - 13 0 13
17 chr1 29408 29408 - 13 0 13
18 chr1 29410 29410 - 13 0 13
19 chr1 29416 29416 - 13 0 13
20 chr1 181056 181056 - 16 16 0
21 chr1 181073 181073 - 15 15 0
22 chr1 181094 181094 - 16 16 0
23 chr1 181102 181102 - 14 14 0
24 chr1 181119 181119 + 13 13 0
25 chr1 181139 181139 + 20 20 0
26 chr1 181140 181140 - 15 15 0
27 chr1 181147 181147 + 20 20 0
28 chr1 181148 181148 - 16 16 0
29 chr1 181150 181150 + 20 14 6
30 chr1 181151 181151 - 16 16 0Thank you!
Regards,
Ben