This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Merge vcf files of structural variants

Hi,

I have many single sample vcf files containing structural variants. I'm trying to merge all the vcfs to one multi-sample vcf, however as the breakpoints of structural variants are only estimates, the position of the same structural variant in different samples would be slightly different. Is there a way to merge individual vcfs containing structural variants into one multi sample vcf? taking account for example a percentage of overlap to call structural variants in different samples as the same variant.

Thanks

structural_variants vcf merge
• 8,228 views • link
updated by Benbird • Benbird • written by ma17 • ma17

2 answers

For that purpose I use SURVIVOR merge, which has a distance setting specifying the maximal difference between pairs of breakpoints (start1 vs start2 and end1 vs end2). As such it doesn't use a percentage of overlap, but it is also not requiring exact breakpoints.

About the percentage overlap the author of SURVIVOR said that he doesn't think it is correct that large SVs easier reach that minimal percentage, and therefore he just uses a distance. I believe he uses a distance of 1kb, but I commonly use 500bp.

• 1 views • link
written by WouterDeCoster WouterDeCoster

Thanks! Exactly what I was looking for

• 1 views • link
written by ma17 • ma17

When an overlap between vcf files A and B is found, which SV representation (i.e. which breakpoints) is reported by SURVIROR?

• 0 views • link
written by bi_Scholar • bi_Scholar

There is one recently published method that could address this question: SV-Pop https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2718-4

• 0 views • link
written by Benbird • Benbird

Log in to answer this question.

More posts like this