I am trying to run delly cnv on my short-read data, but it fails right at the beginning with this error:
[W::hts_idx_load3] The index file is older than the data file: sample1.bai
[E::fai_build_core] Format error, unexpected "1" at line 1
.command.sh: line 2: 45 Segmentation fault (core dumped)
delly cnv -a -g hg38.fa -m mappability.bed -c sample1.cov.gz -o sample.bcf sample1.bam I checked all input files if chromosome names match, everything is ok.
Could the bed file be the reason? I was also trying to create mappability map from the files on
https://gear-genomics.embl.de/data/delly/,
Homo_sapiens.GRCh37.dna.primary_assembly.fa.r101.s501.blacklist , but I am not sure how exactly to extract the map from this files?
And as for the indexing error, my bai and fai files should be also alright, since I ran the delly somatic SV call with them and everything worked just fine.
Since I am using nextflow for delly pipeline (docker), I wanted to run delly cnv in the terminal on conda env, but I was not able to update delly to higher version then 0.7.6, even when using mamba and bioconda channel...
Thank you very much for your help.
Best,
Patricie