-
Notifications
You must be signed in to change notification settings - Fork 27
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
[fastq_to_ubam] Cannot determine candidate qualities: no qualities found. #415
Comments
Ultimately this is an error from Picard rather than viral-pipelines, but I'm curious about what may be different about your input files. Would it be possible for you to post them or an excerpt? |
Well I think what's weird is that she's using the same input files from the repo's I mean this exact miniwdl invocation passed with the exact same input a few days ago in GHActions. It is almost as if the fastqs got corrupted after tarring or something. |
The input files are the provided ones which look like this:
and
I get the same error whether I run the script as in my post above, as well as when I run the script (from a newer version of viral-pipelines) like this:
I get the same error using the fastq file attached here |
Wait, Shandu are you running on a Mac? Is your Docker using "gRPC FUSE" or
"legacy osxfs" for file mounts? I wonder if this is the miniwdl/docker/mac
bug where all the input files look empty to everything....
…On Tue, May 3, 2022 at 10:15 AM shandukani ***@***.***> wrote:
The input files are the provided ones which look like this:
@myseq/1
TCAATAAAAAAAAAAAAGAAAGAAAAAAAAATTCTCCTCATTTTTGTTGT
+
""""""""""""""""""""""""""""""""""""""""""""""""""
and
@myseq/2
AATTATATTATTTCTTTGATAATTTCCTCTCCTCTTGTTTCTTTGTTTCT
+
"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#"#
I get the same error whether I run the script as in my post above, as well
as when I run the script (from a newer version of viral-pipelines) like
this:
miniwdl run https://raw.githubusercontent.com/broadinstitute/viral-pipelines/v2.1.8.0/pipes/WDL/workflows/fastq_to_ubam.wdl FastqToUBAM.fastq_1=../pipeline/viral-pipelines-2.1.0.2/test/input/in1.fastq \
FastqToUBAM.fastq_2=../pipeline/viral-pipelines-2.1.0.2/test/input/in2.fastq \
FastqToUBAM.sample_name=in \
FastqToUBAM.library_name=mylib
—
Reply to this email directly, view it on GitHub
<#415 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/ACA6RUWBOEFOFMALEXMNBQDVIEYJJANCNFSM5US3HYHQ>
.
You are receiving this because you commented.Message ID:
***@***.***>
--
Daniel J. Park, PhD *(he/him/his)*
Sr Group Leader, Viral Computational Genomics
Broad Institute of MIT and Harvard
Cambridge, MA, USA
Tel: +1-617-714-8526
***@***.***
https://www.broadinstitute.org/bios/daniel-park-phd
|
I am running on a Mac and appear to be using gRPC FUSE for file sharing. Should I be using legacy osxfs? |
Yeah, try that (you'll need to restart your Docker environment) |
Hmm, tried that but now I get a DockerException
Am I supposed to be explicitly running a docker container before trying to run the WDL file, or is that handled for me when I run the WDL? |
miniwdl will start (and stop) the docker containers as appropriate when working through workflow tasks. The socket refusing the connection makes me think the full docker engine process is not running. You may need to update/restart the docker engine and try again? As long as gRPC FUSE is disabled the file system mounting seems fairly reliable on macOS. |
Restarted docker again, as well as updated from 4.7.0 to 4.7.1, and it works (so issue was file system mounting and/or version issue)! Thank you both for your help! |
(1) I downloaded the workflow files locally with
(2) I ran fastq_to_ubam with the provided test data:
(3) I get the following stderr file, which seems to suggest that the qualities cannot be read in properly. The same error occurs when I use a different input fastq files
stderr.txt
The text was updated successfully, but these errors were encountered: