Error message with medaka pyabpoa and libgsl.so.25

124 Views Asked by At

I am using medaka in hybracter pipeline. Hybracter creates different conda environments for each step of its pipeline, and notably uses medaka for polishing. When I run the hybracter pipeline I get this error message :

[Mon Jan 15 13:18:40 2024]
Finished job 11.
3 of 18 steps (17%) done
Checking program versions
This is medaka 1.11.3
Program    Version    Required   Pass     
bcftools   Not found  1.11       False    
bgzip      1.17       1.11       True     
minimap2   2.26       2.11       True     
samtools   1.18       1.11       True     
tabix      1.17       1.11       True     
[Mon Jan 15 13:18:41 2024]
Error in rule medaka_round_1:
    jobid: 25
    input: /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/pre_polish/hybracter_RB04_test_chromosome_plus_plasmids.fasta, /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/qc/hybracter_RB04_test_filt_trim.fastq.gz
    output: /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/complete/medaka_rd_1/hybracter_RB04_test/consensus.fasta, /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/versions/hybracter_RB04_test/medaka_complete.version
    log: /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/stderr/medaka_round_1/hybracter_RB04_test.log (check log file(s) for error details)
    conda-env: /home/horigene/anaconda3/envs/hybracterENV/lib/python3.12/site-packages/hybracter/workflow/conda/348d2d1aa0134c8a68065158ce832239_
    shell:
        
        medaka_consensus -i /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/qc/hybracter_RB04_test_filt_trim.fastq.gz -d /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/pre_polish/hybracter_RB04_test_chromosome_plus_plasmids.fasta -o /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/complete/medaka_rd_1/hybracter_RB04_test -m r1041_e82_400bps_sup_v4.2.0  -t 16 2> /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/stderr/medaka_round_1/hybracter_RB04_test.log
        medaka --version > /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/versions/hybracter_RB04_test/medaka_complete.version
        touch /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/complete/medaka_rd_1/hybracter_RB04_test/calls_to_draft.bam
        rm /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/complete/medaka_rd_1/hybracter_RB04_test/calls_to_draft.bam
        touch /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/complete/medaka_rd_1/hybracter_RB04_test/consensus_probs.hdf
        rm /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/processing/complete/medaka_rd_1/hybracter_RB04_test/consensus_probs.hdf
        
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)
Logfile /mnt/data/horigene/Corentin/20231214_24_genomes_GW/hybracter/stderr/medaka_round_1/hybracter_RB04_test.log:
================================================================================
Cannot import pyabpoa, some features may not be available.
Cannot import pyabpoa, some features may not be available.
Cannot import pyabpoa, some features may not be available.
Cannot import pyabpoa, some features may not be available.
bcftools: error while loading shared libraries: libgsl.so.25: cannot open shared object file: No such file or directory
================================================================================

Shutting down, this might take some time.
Exiting because a job execution failed. Look above for error message
Complete log: .snakemake/log/2024-01-15T131838.221728.snakemake.log
WorkflowError:
At least one job did not complete successfully.

Have you ever had this problem or do you know how to fix it?

Thanks in advance

Corentin

I tried to install pyabpoa using the command conda install pyabpoa but the error message remains. I tried to update bcftools with the conda install bcftools command but I still get the error message

0

There are 0 best solutions below