Error while generating index files using samtools

Question

[bisht20diksha]

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?

Answer 1

[Devon Ryan]

samtools index -c test_sorted.bam