I'm just trying to index a sorted BAM with
samtools index but get the error:
[E::hts_idx_push] Region 536892210..536892211 cannot be stored in a csi index with min_shift = 14, n_lvls = 5. Try using min_shift = 14, n_lvls >= 6 samtools index: failed to create index for "868om_41_070_010.sorted.bam": Numerical result out of range
no .bai nor .cis is made by the command, as is suggested by enter link description here. I think i found where n_lvls is defined in hts.c (lines 1817, 1829, 1847). However I'm not sure I how I should proceed, i.e. what compiler and options should be used to compile hts.c after I change n_lvls.