This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Merging consecutive non-variant calls

I'm looking for a tool that can merge non-variant segments from GATK together. Essentially what I'm looking for is the opposite of gvcftools' break_blocks command. Anyone know of anything?

gvcf vcf gatk gvcftools
• 1,942 views • link
updated by Chris Saunders Chris Saunders • written by gmcinnes • gmcinnes

1 answer

It has not been maintained for current GATK, but this tool is designed to merge any non-variants from GATK 1.x all-sites output:

https://sites.google.com/site/gvcftools/home/configuration-and-analysis/gatk_to_gvcf-usage

• 1 views • link
written by Chris Saunders Chris Saunders

Thanks, Chris. What if I already have gvcfs with non-variant blocks, but want to use less stringent criteria to merge the non-variant blocks together?

• 1 views • link
written by gmcinnes • gmcinnes

Sorry, don't know of anything that would do this right out of the box..

• 1 views • link
written by Chris Saunders Chris Saunders

Log in to answer this question.

More posts like this