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

Not able to merge insertions #199

Open
beantkapoor786 opened this issue Aug 30, 2023 · 5 comments
Open

Not able to merge insertions #199

beantkapoor786 opened this issue Aug 30, 2023 · 5 comments

Comments

@beantkapoor786
Copy link

Hello Fritz,

We have a tree genome pedigree in which we are comparing 5 different genomes to the reference genome. We aligned each genome using mummer4 and called structural variants. Now, we want to merge those structural variants using SURVIVOR's merge function. We were able to merge deletions and inversions but cannot merge insertions. Is it because the start and end positions in insertions are in the format END=START+1 and SVLEN=1? Please provide your thoughts. I am using the following command to do the merge but output vcf file is empty.

SURVIVOR merge insertions.txt 10000 1 1 1 0 1 mummer_insertions_merged.vcf

Output message -

merging entries: 51273
merging entries: 79656
merging entries: 75335
merging entries: 77095
merging entries: 79632

I have also attached the insertion vcf files here in case you want to look at them. Thanks!
insertions.zip

@sd880131
Copy link

I have the same problem

@fritzsedlazeck
Copy link
Owner

Hey there,
mmh the SVLen =1 is an issue. I wrote that wrapper script a while back for something... are you using this converter script from SURVIVOR for mummer ?

Btw. I would suggest this: (10k is a bit big)
SURVIVOR merge insertions.txt 10000 1 1 0 0 1 mummer_insertions_merged.vcf

The strand information might also be causing the issue...

Thanks
Fritz

@JiaB-He
Copy link

JiaB-He commented Dec 25, 2023

Hey there, mmh the SVLen =1 is an issue. I wrote that wrapper script a while back for something... are you using this converter script from SURVIVOR for mummer ?

Btw. I would suggest this: (10k is a bit big) SURVIVOR merge insertions.txt 10000 1 1 0 0 1 mummer_insertions_merged.vcf

The strand information might also be causing the issue...

Thanks Fritz

Hello, I did not find this conversion script

@zhiyongli1995
Copy link

I have the same problem. @beantkapoor786 Have you solved the problem? If so, please share the method. Thank you!

@beantkapoor786
Copy link
Author

Hello, apologies for the late response. I couldn't do it using survivor but found another tool to merge the SVs - https://github.com/DecodeGenetics/svimmer

Hope this helps.

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

5 participants