Skip to content
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

Closed
shandu-m opened this issue Apr 28, 2022 · 9 comments
Closed

Comments

@shandu-m
Copy link

(1) I downloaded the workflow files locally with

wget https://github.com/broadinstitute/viral-pipelines/archive/v2.1.0.2.tar.gz

(2) I ran fastq_to_ubam with the provided test data:

miniwdl run ../pipeline/viral-pipelines-2.1.0.2/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

(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

@tomkinsc
Copy link
Member

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?

@dpark01
Copy link
Member

dpark01 commented Apr 29, 2022

Well I think what's weird is that she's using the same input files from the repo's test/input/ directory that our CI test already does and passes regularly on GHActions on both miniwdl and Cromwell. I mean she's using an older release than usual, but I'm pretty sure all these tests passed back then too…

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.

@shandu-m
Copy link
Author

shandu-m commented May 3, 2022

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

I get the same error using the fastq file attached here
fq_18.fastq.zip

@dpark01
Copy link
Member

dpark01 commented May 3, 2022 via email

@shandu-m
Copy link
Author

shandu-m commented May 3, 2022

I am running on a Mac and appear to be using gRPC FUSE for file sharing. Should I be using legacy osxfs?

@dpark01
Copy link
Member

dpark01 commented May 3, 2022

Yeah, try that (you'll need to restart your Docker environment)

chanzuckerberg/miniwdl#145 (comment)

@shandu-m
Copy link
Author

shandu-m commented May 3, 2022

Hmm, tried that but now I get a DockerException

docker.errors.DockerException: Error while fetching server API version: 500 Server Error for http+docker://localhost/version: Internal Server Error (\"b'dial unix /Users/shandu/Library/Containers/com.docker.docker/Data/docker.raw.sock: connect: connection refused'\")

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?

@tomkinsc
Copy link
Member

tomkinsc commented May 3, 2022

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.

@shandu-m
Copy link
Author

shandu-m commented May 3, 2022

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!

@shandu-m shandu-m closed this as completed May 3, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants