-
Notifications
You must be signed in to change notification settings - Fork 23
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
SIGSEGV: Illegal storage access. (Attempt to read from nil?) #27
Comments
can you show the command you are running? make sure you don't have repeated labels like:
|
Hi Najeeb, I'd like to get this resolved, can you give more information? |
Hi Brent , Extremely Sorry for late reply !! sh slivar.sh 54Samplesvcf2db.vcf.gz batch5.ped missinginped |
would you try runnng the command using this binary (just gunzip and chmod +x) so we can see more debug info? |
54 samples matched in VCF and PED to be evaluated
Traceback (most recent call last)
slivar.nim(200) slivar
slivar.nim(197) main
slivar.nim(105) expr_main
gnotate.nim(71) open
zipfiles.nim(194) extractFile
zipfiles.nim(182) extractFile
streams.nim(88) readData
SIGSEGV: Illegal storage access. (Attempt to read from nil?)
slivar version: 0.1.2
54 samples matched in VCF and PED to be evaluated
Traceback (most recent call last)
slivar.nim(200) slivar
slivar.nim(197) main
slivar.nim(105) expr_main
gnotate.nim(71) open
zipfiles.nim(194) extractFile
zipfiles.nim(182) extractFile
streams.nim(88) readData
SIGSEGV: Illegal storage access. (Attempt to read from nil?)
…On Fri, Apr 26, 2019 at 7:17 PM Brent Pedersen ***@***.***> wrote:
would you try runnng the command using this binary (just gunzip and chmod
+x) so we can see more debug info?
slivar_debug.gz
<https://github.com/brentp/slivar/files/3122143/slivar_debug.gz>
—
You are receiving this because you authored the thread.
Reply to this email directly, view it on GitHub
<#27 (comment)>, or mute
the thread
<https://github.com/notifications/unsubscribe-auth/ABBFINYPDWZPES4VD5P3AD3PSMTJVANCNFSM4HICPXHA>
.
--
Syed Najeeb Ashraf
|
ok. that's easy. You got a zip file from an old release. If you get one linked from the latest release here: https://github.com/brentp/slivar/releases you'll be fine. Sorry for the trouble. |
to clarify, by that, I mean grab the latest gnomad zip (GRCh37) |
Hi Brent, Now I am getting below error. slivar version: 0.1.2 54 samples matched in VCF and PED to be evaluated
|
Hi Najeeb, sorry for the trouble. Can you share a small part of the VCF that I can use to recreate the error? And do you get this with the original binary that you downloaded (or only the debug binary)? I can't see anything obvious to change just yet. |
@brentp - I also come across this error Not sure exactly why it happens, but seems to not like VCFs coming from GATK CombineVariants or after using slivar and using the |
if I remove all the |
This is needed to remove GATK.CommandLine from VCF header. slivar errors out when these are present. See: brentp/slivar#27 (comment)
can you give me a small VCF to recreate? |
Yup here here is an example VCF |
thank you. how can I see the problem. I tried:
without problem. |
PED file:
|
I used version |
if I put that in x.ped:
so, I'm not able to recreate. |
that's odd, I get:
|
what is your test.vcf? |
well, when it is uncompressed I get the error. I had to bgzip it when I uploaded to github haha. |
oh. wow. I get the error on the uncompressed, but not compressed. will fix. |
actually, this is already fixed in my dev branch will be out in next release. thanks for reporting. |
Thank you! Yea it was helpful to try to figure out why it was happening. What eventually fixed it for me was to compress VCF or to remove |
This is needed to remove GATK.CommandLine from VCF header. slivar errors out when these are present. See: brentp/slivar#27 (comment)
@brentp Hi Brent, Following up on this issue. We are using slivar v0.1.10 and see this issue.
But the thing is, I cannot recreate this issue since whenever I run it again, it passes through. We pass compressed vcf to slivar expr Do you have any suggestions?
|
can you share a small vcf and ped to recreate? you can try replacing "NA" with "0" or "-9" since those are accepted missing values in a pedigree file. |
@Tayazseven can you make sure you have latest release (version 0.1.12). |
@Phillip-a-richmond can you link a VCF with a RefCall variant that does have this problem? |
Okay so while trying to isolate the issue to the RefCall, I was unsuccessful. It turns out the issue that I solved with my zgrep for RefCall, was actually just because I uncompressed from vcf.gz to plain text, and then reran the slivar command on the uncompressed VCF. If I had to say where the issue is exactly, my guess is that it has something to do with the bgzipped VCF I'm using being read into slivar. So what I can reproduce with my files:
Error:
But this code does not give an error:
I'm attaching here the uncompressed and compressed VCFs, with least 10000 lines in each of them. Hopefully that can help your debugging process. Note, change the vcf.txt back to just vcf...github won't let me paste . |
thanks very much @Phillip-a-richmond , I will have a look. |
Not sure if it is for the same reasons, but I have a VCF that is generating the same error that I could point you to if that would be helpful. |
I can get this with the last release. With dev version, it seems to be fixed--though I'm not sure which change is to blame/praise for that. If so, I'll make a new release soon. Thanks. |
... and now with the binary attached. |
That binary didn't fix the problem for me.
The VCF is bgzip'd with tabix 1.10.2 and indexed. On the uncompressed VCF:
Is there a verbose/debug option you've got working on a binary you can share? Perhaps we can pinpoint where the complaint is on my end? Cheers, |
@Phillip-a-richmond , can you give the command that generates the error with the files that you posted above? |
Hi- Has this problem been solved? I am trying to implement slivar for patient data, so unfortunately I cannot share any VCFs. I have tried both with unzipped and bgzipped files, the most recent build, and the binary shared above. I get this error after passing my vcf through bcftools csq Thanks, |
I may have dropped the ball on this one. Raungar the error went away for me when I cut out the RefCall lines from my VCF. Below is a functional pipeline starting with a merged VCF from deeptrio-->glnexus-->bcftools (merged.vcf.gz)
I recommend grabbing an open source trio to share your errors. Plenty available, but I test with this one:
|
Hi @raungar , I'd like to get this resolved, but I have not been able to recreate. Can you try running the same command with this binary attached and let me know what the full output is? @Phillip-a-richmond if you can show me how to recreate, I'll look into it. Those header lines should not affect slivar. I am curious what a few of the RefCall lines look like so if you can paste those perhaps that will also give a clue. |
@raungar binary attached this time: (just download, gunzip, chmod +x slivar_debug, thne use |
I think I found the problem. Would you try this binary as well? |
Hi Brent and everyone, I believe I am having the same issue. First, those two binaries that you just posted did not solve the issue for me. But I do have some info that might help. Oddly this error will occur most of the time, but when repeatedly trying the command, it will rarely run successfully (maybe 10% of the time). For me at least, removing the RNC format field from glnexus VCFs seems to fix the error. From what I can tell the RNC format field does not follow the VCF spec. Attached are vcfs created using glnexus "get started page" and then removed the FORMAT/RNC |
@JakeHagen ok. this is a great start. Now, can you share a command that generates the sigsegv?
but it completes without issue (though it does give warning about multiallelics) |
The RNC field is valid VCF so something else is going on. If I can recreate the bug then I should be able to fix, but I don't see it yet. |
Hi- I figured out the problem for my VCF at least. This was purely a header issue for me -- I didn't have all the FORMAT tags in my header, which gave warnings but no errors. Adding in the proper FORMAT tags (ex: ##FORMAT=<ID=GQ,Number=1,Type=Integer,Description="Genotype Quality">) solved this problem for me. Thanks, |
@brentp Hi Brent, the command you used above works fine for me. Below is the command that errors for me.
The error usually only comes up when using |
Also, now that you confirmed RNC is valid (wasn't familiar with character fields), I tried removing other format fields. |
I can recreate with the command above. I will have a fix out in next release soon. Thanks for the example @JakeHagen |
Thats great, thank you |
OK. I have the fix in for this. Will add test and commit this weekend. Meanwhile, here is the binary. If you could verify that it works for you that would be great. I expect I can get a release out Sunday or Monday. Thanks so much for reporting everyone and for providing the test-case @JakeHagen , that is critical. This is a debug binary so will be slow, but should work for the cases this issue is referencing. |
Yep that seems to do the trick. Thanks a lot for this. |
this is finally, fully fixed in v0.2.4 |
Hey @brentp I'm running into this error again I think, only when using
and the output
And here's the input VCF that was annotated with Had the same error with multiple, now up to 3, VCFs. But works fine on a simple chr21 only test VCF |
Hi Matt, thanks for the test-case. You have this INFO:
at chr6:63280472 G/A. This is because you ran bcftools csq with |
Worked perfectly with |
Hi Brent,
I am trying to run the latest slivar version: 0.1.1 via a shell script. I am getting below error.
54 samples matched in VCF and PED to be evaluated
SIGSEGV: Illegal storage access. (Attempt to read from nil?)
slivar version: 0.1.1
When I am trying to run the same script with Version 0.0.5, I am not getting any error.
Do you have any idea that what I am getting this error?
Thanks
Najeeb
The text was updated successfully, but these errors were encountered: