This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RSeQC : infer_experiment.py Error: Could not retrieve index file

I am trying to use RSeQC for identifying strandness but I am getting error as Could not retrieve index file

Code

infer_experiment.py --i test-sorted.bam -r output.bed12

Output:

[E::idx_find_and_load] Could not retrieve index file for 'test-sorted.bam'
Reading reference gene model output.bed12 ... Done
Loading SAM/BAM file ...  Total 200000 usable reads were sampled

Anyone knows how to address this issue?

rseqc rna-seq

1 answer

You need to have indexed your bam/sam file using samtools index

thanks, it's working now.

A small educational note: if an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they all work. If an answer was not really helpful or did not work, provide detailed feedback so others know not to use that answer.

upvote_bookmark_accept

Log in to answer this question.