delly CNV error

82 views
Skip to first unread message

PatS

unread,
Aug 1, 2024, 11:26:55 AM8/1/24
to delly-users
Hi group,

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
 
Reply all
Reply to author
Forward
0 new messages