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

unexpected command line complaint with -r #34

Closed
AMChalkie opened this issue May 31, 2019 · 6 comments
Closed

unexpected command line complaint with -r #34

AMChalkie opened this issue May 31, 2019 · 6 comments
Assignees

Comments

@AMChalkie
Copy link

java -Xmx6g -jar JACUSA_v2.0.0-RC5.jar call-1 -b Sppl2a.bed -f V -r Sppl2a.jacusa2.vcf myfile.bam
INFO 00:00:00 Computing overlap between sequence records.

JACUSA2 Version: 2.0.0-RC5 call-1 -b Sppl2a.bed -f V -r Sppl2a.jacusa2.vcf myfile.bam

java.lang.IllegalArgumentException: Cannot set a file type if the output is not to a file.
at htsjdk.variant.variantcontext.writer.VariantContextWriterBuilder.setOutputFileType(VariantContextWriterBuilder.java:185)
at jacusa.io.format.call.VCFcallWriter.(VCFcallWriter.java:59)
at jacusa.io.format.call.VCFcallFormat.createWriter(VCFcallFormat.java:41)
at lib.worker.WorkerDispatcher.(WorkerDispatcher.java:53)
at lib.util.AbstractMethod.getWorkerDispatcherInstance(AbstractMethod.java:60)
at lib.util.AbstractTool.run(AbstractTool.java:64)
at jacusa.JACUSA.main(JACUSA.java:99)
at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:62)
at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
at java.lang.reflect.Method.invoke(Method.java:498)
at org.eclipse.jdt.internal.jarinjarloader.JarRsrcLoader.main(JarRsrcLoader.java:61)

@AMChalkie
Copy link
Author

Note that without the -f V, an output file is produced.

@ellalalalalalala
Copy link

ellalalalalalala commented Jul 14, 2021

Hi AMChalkie, how did you solve the error? If I remove the -f V option, will my output then still be VCF, given the file ending of the output file name? In the manual it says, that the default output is BED. I did not have any comparable problems with the combination of the "-f V" and "-r file.out" option, when I used Jacusa1...

@ellalalalalalala
Copy link

I tried by removing -f V and labeling the output file as FILE.vcf... no error anymore but results in an empty output file :(

@piechottam piechottam self-assigned this Jul 15, 2021
@piechottam
Copy link
Collaborator

I have found the problem and I will publish a new release with the fix today.

Please post issues in:
https://github.com/dieterich-lab/JACUSA2

This repository is for JACUSA1

@ellalalalalalala
Copy link

Hi piechottam, great, thanks a lot for your fast reply and help 🥇 Sorry for posting in the wrong repository, ich just found this thread and it was exactly my problem. :)

@piechottam
Copy link
Collaborator

Follow dieterich-lab/JACUSA2#39

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