This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Comparing calls from different VCFs

Hi all,

Does anyone know of a way to compare mutation calls for two vcf files? The vcfs should be derived from the same dataset but I may have used a different caller than the other vcf, so I want to compare if the calls were good or what do they filter in general.

vcf maf

That doesn't tell me what calls weren't included in the intersection that were in either vcf alone though, does it?

From the BEDtools 'intersect' documentation:

There will likely be cases where you’d like to know which “A” features do not overlap with any of the “B” features. Perhaps you’d like to know which SNPs don’t overlap with any gene annotations. The -v (an homage to “grep -v”) option will only report those “A” features that have no overlaps in “B”.

0 answers

No answers yet.

Log in to answer this question.