search

zhengxwen / SeqArray

Data management of large-scale whole-genome sequence variant calls (Development version only)
http://www.bioconductor.org/packages/SeqArray
43 stars 12 forks source link

Duplicate indices in `seqSetFilter(f, variant.sel=)` #57

Closed zhengxwen closed 4 years ago

zhengxwen commented 4 years ago 
library(SeqArray)

# open the GDS file
f <- seqOpen(seqExampleFileName("gds"))

seqSetFilter(f, variant.sel=c(1:50, 1:50))

Fails to select the top 50 variants:

# of selected variants: 100

Incorrect read of selected variants

seqGetData(f, "chromosome")
  [1] "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1"
 [19] "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1"
 [37] "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" "1" ""  ""  ""  ""
 [55] ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""
 [73] ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""  ""
 [91] ""  ""  ""  ""  ""  ""  ""  ""  ""  ""