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

Error loading vcf with some SV types #26

Open
cboursnell opened this issue Nov 21, 2019 · 1 comment
Open

Error loading vcf with some SV types #26

cboursnell opened this issue Nov 21, 2019 · 1 comment

Comments

@cboursnell
Copy link

cboursnell commented Nov 21, 2019

I have a GRIDSS vcf that I have updated the svtype field (which was all BND to start with) to be taken from the output of Clove. So these are BND, DEL, INV, TAN, TRA.
When I try to load this vcf into R and get the breakpointRanges I get an error.
My R is:

vcf <- readVcf(input, "hg38")
gr <- breakpointRanges(vcf)

and the error is

Error in names(cgr)[is.na(end), ] : incorrect number of dimensions

Update: I have tried just changing the deletion events from BND to DEL and this gives the same error. I can provide a minimal reproducible example if needed.

@d-cameron
Copy link
Member

Yes, please provide a minimal reproduction VCF with just one or two of the offending records.

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

2 participants