-
Notifications
You must be signed in to change notification settings - Fork 7
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Preprocessing warning message #10
Comments
That's odd for a BAM - HTSlib seems to be treating this file as though it's a CRAM, then throwing an error when it tries to download reference sequences from ENA. Is there any chance I could get you to please provide one of the commands that's generating the error? |
Hi, thanks for the answer, actually some of the files in my cohort are cram and others are bam, however, in the manual I saw that I should use "--mode=bam" for both cram and bam; furthermore the error occurs only for some cram files and not other cram. Here's the command I use:
Thanks |
Yep! You've got the correct command. The issue is in htslib and how it reads CRAMs (which superSTR relies on, as well as samtools and many other pieces of software). htslib is failing to download the reference sequence that you've aligned to for some reason from ENA; this happens sometimes transiently I've found, but can also happen in server/workflow contexts. If you're aligned to GRCh38, running the following should work (from here):
Full documentation is available here. If you're running in an HPC, you'll need to make sure that REF_CACHE and REF_PATH are correctly set in the job environment. Unfortunately, you'll need to re-run those CRAMs for which the error was generated :( |
Oh, as a note to myself - I'll add some comments on this in the documentation before closing this out. If there's any other issues you hit with superSTR let me know! |
Hi, I ran superSTR preprocessing on a cohort of WGS samples (in bam format). The output seems ok for all samples, however, for some of them I get this message in the log:
Is it safe to proceed with the analysis anyway?
Thanks.
The text was updated successfully, but these errors were encountered: