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

Bug found: SV LEN value cutted #213

Open
leone93 opened this issue Mar 26, 2024 · 0 comments
Open

Bug found: SV LEN value cutted #213

leone93 opened this issue Mar 26, 2024 · 0 comments

Comments

@leone93
Copy link

leone93 commented Mar 26, 2024

Hy Fritz,
I just found his bug trying to merge the SV coming from 5 different callers (pbsv,svim,sniffles,cutesv,svim-asm) for the same sample.

This are only two example (just a part of the vcf line) but the files are full of that:
Chr11 5514126 svim_asm.INS.11708,pbsv.INS.27728 G GAATACATATATCCAAGGAAAAAGTACGAATTACACCCC> SVLEN=84
Chr11 5514126 svim.INS.17070,cuteSV.INS.2923 G GAATACATATATCCAAGGAAAAAGTACGAATTACACCCCTGAACTAT> SVLEN=1584

Chr11 5710087 svim_asm.INS.11719,pbsv.INS.27758 G GCTATTTGTAACACTCTGAAAATTCGACCGACAAAATCA> SVLEN=60
Chr11 5710087 svim.INS.17086,Sniffles2.INS.162SA,cuteSV.INS.2935 G GGCTATTTGTAACACTCTGAAAA> SVLEN=3760

Two notes: the original file are parsed/simplified to make the merging work of survivor more easy reducing most of the field.
Where the error is present, the input files are perfect. There the original value is as for the other 3 caller, 1584 and 3760.
It seems only that for the first two, pbsv and svm-asm, (the last 2 in the sample list that I give to survivor) it cut the the first two number of the SVLEN, so from 1584 to 84, from 3760 to 60; to me it is very weird but it is everywhere. Any idea?

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

1 participant