LCA

[TZGyalsey]

I tried to run the following code:

for file in *.sam.gz; do
filename=$(basename "$file" .sam.gz )
$metaDMG lca --names "$nam" --nodes "$nod" --acc2tax "$acc" --sim_score_low 0.95 --sim_score_high 1.0 --min_mapq 10 --bam "$file" --lca_rank family --used_reads [0,1] --weight_type [0,1]
done

I got this error:

-> metaDMG version: v0.4-48-gddf8f79-dirty (htslib: 1.17-28-ge13611a) build(Jan  4 2024 17:06:48)
-> /work/home/acz498gefh/programmes/metaDMG-cpp/metaDMG-cpp lca --names /work/home/acz498gefh/ref/names.dmp --nodes /work/home/acz498gefh/ref/nodes.dmp --acc2tax /work/home/acz498gefh/ref/combined_taxid_accssionNO_20223103.txt --sim_score_low 0.95 --sim_score_high 1.0 --min_mapq 10 --bam TTH25.col.fq.merged.sorted.sam --lca_rank family --used_reads [0,1] --weight_type [0,1] 
-> Will output lca results in file:     'outnames.lca.gz'
-> [thread1] Will read header
-> Will output lca distribution in file:        'outnames.stat.gz'

[E::hts_hopen] Failed to open file TTH25.col.fq.merged.sorted.sam.gz [E::hts_open_format] Failed to open file "TTH25.col.fq.merged.sorted.sam.gz" : Exec format error /opt/gridview/slurm/spool_slurmd/job8381806/slurm_script: line 23: 55885 Segmentation fault $metaDMG lca --names "$nam" --nodes "$nod" --acc2tax "$acc" --sim_score_low 0.95 --sim_score_high 1.0 --min_mapq 10 --bam "$file" --lca_rank family --used_reads [0,1] --weight_type [0,1]

how should I solve this problem? [E::hts_hopen] Failed to open file TTH25.col.fq.merged.sorted.sam.gz [E::hts_open_format] Failed to open file "TTH25.col.fq.merged.sorted.sam.gz" : Exec format error

I checked I installed htslib.