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

Invalid END position #497

Open
ACEnglish opened this issue Jul 18, 2024 · 3 comments
Open

Invalid END position #497

ACEnglish opened this issue Jul 18, 2024 · 3 comments
Assignees

Comments

@ACEnglish
Copy link

Using sniffles 2.4 in merge mode on chm13 produced an invalid variant.

This isn't a huge problem because most vcf parsers can handle the error plus SVs at the beginning of chromosomes should be a rare occurrence(?). But I just wanted to report it.

Error:
[W::vcf_parse_info] INFO/END=0 is smaller than POS at chr11:1

Offending Line:
chr11   1       Sniffles2.INS.8M0       t       tCCCAACCCCACCCGAACTCCGAACCCCAACTCCGAACCCCAACTCGAACCCGAACCCTAACCCGAACCTTGAATTCGAACT50      PASS    PRECISE;SVTYPE=INS;SVLEN=87;END=0;SUPPORT=3;COVERAGE=10,13,14,15,23;STRAND=+-;AC=22;STDEV_LEN=10.746;STDEV_POS=393.116;SUPP_VEC=0000000000000000000010000000000000000000000000000000000000000000000000000010000000000000010000100000000000000000000000000000000000001000000000010000000000000000001000000100001000000000000000000010000000000100100000000000000000000000000000000000000000000000000000000000000000000000010000000000000010000000000000000000000011000000000000000000000000000000000000010000010000000000000000000000100100000000100000000000000000000000000010000000000000000000000000000100000000000000000000000000011000000100000000000000001000000000010000000000000000100000000000000000000101000000100000010000000000000000000000000000010000000000000000000000000000000000000000000000001000000000000000000000000000000000000000000000011000000000000001000000000000001001000000000001000000000000000000000000000000000000000000000000000000000000000000000100000000000000000000000000000000000000000000000001000101000000000000010000000000000000000000000000000000000000000000000000000000000000000100000001010000000000001000000000000000000000000100000000000000000000000001001100000000000000010000000100000000000000100000000000000   GT:GQ:DR:DV:ID  0/0:0:25:0:NULL
@hermannromanek
Copy link
Collaborator

Hi @ACEnglish,

Which dataset did you see this on?

Thanks,
Hermann

@ACEnglish
Copy link
Author

ACEnglish commented Aug 12, 2024

This was in the 1kgp ONT data. Fritz has the snf files locally. https://www.biorxiv.org/content/10.1101/2024.04.18.590093v1.full.pdf

@lfpaulin
Copy link
Collaborator

@hermannromanek it looks like Sniffles is diminishing by 1 the end when outputting the END position and for INS its the same as start that when is 1 it becomes zero

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