This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error while generating index files using samtools

Hello. I am trying to identify differential alternative splicing between samples using splAdder. It requires sorted bam files and index files. I have sorted bam files individually using samtools script

samtools sort test.bam test_sorted

Now, I am trying to generate index file using script:

samtools index test_sorted.bam test_sorted.bai

But everytime following error appears:

[E::hts_idx_push] Region 536884677..536885767 cannot be stored in a bai index. Try using a csi index with min_shift = 14, n_lvls >= 6
samtools index: failed to create index for "cd_cnt_rt_rep1Aligned.sorted": Numerical result out of range

How to fix it?

spladder samtools

2 answers

samtools index -c test_sorted.bam
samtools index -b input.bam

That's the default...

Log in to answer this question.