Methylation profile ERROR for large genome

[Yang990-sys]

Hello, Help please,

ERROR: can't load csi index，try to rebuild bai index and report bug

MethylDackel version: 0.6.1 (using HTSlib version 1.21) Installed by conda

script:

bsmap -a sample.R1.clean.fq.gz  -b  sample.R2.clean.fq.gz -d ref_bis/all.fa -p 20 -v 0.05 -o align.bam 2>bsmap_report.txt

sambamba sort -m 8GB --tmpdir tmp -t 10 -o ${prefix}_sorted.bam ${prefix}.bam
sambamba markdup --overflow-list-size 1000000 --tmpdir tmp -t 10 ${prefix}_sorted.bam dedup.bam 2> markdup_report.txt

MethylDackel extract ref_bis/all.fa dedup.bam --opref prefix

Thanks a lot