This is a test version of Biostars. For the public version, visit https://www.biostars.org.
ERROR: "Fragment file is not indexed"

Running a pipeline with the following argument:

chrom_assay <- CreateChromatinAssay(
   counts=counts,
   sep = c(":", "-"),
   genome = 'mm10', 
   fragments = '.../fragments.tsv.gz.tbi', 
   min.cells = 10,
   min.features = 200)

I keep receiving the error `Fragment file is not indexed. Any suggestions?

atac-seq cellranger snatac-seq

1 answer

There is no such a directory as .../ - maybe you meant ../ for one directory up? Either that is wrong, or you don't have the original fragments file (fragments.tsv.gz) in the same directory.

https://github.com/timoast/signac/issues/7

Right; I was demonstrating that I was providing the path of the fragment file. I have tried .tsv, .tsv.gz, .tbi files, etc. Same error.

So both files ../fragments.tsv.gz.tbi and ../fragments.tsv.gz do exist and you are allowed to read them?

Log in to answer this question.