Methods for Doing Consensus Variant Calling between multiple callers
1
1
Entering edit mode
2.2 years ago
djevanclark ▴ 40

We are currently retooling a pipeline to do consensus calling between 3 variant callers, Mutect2, strelka, and lancets. We would like some thoughts on the best way to do consensus calling on the VCF outputs. We have some older custom scripts that are very convoluted, and have also looked at selectvariants from GATK. Any input is appreciated.

WES DNA Consensus Variant Calling • 1.0k views
ADD COMMENT
0
Entering edit mode

This work might give some insipration: https://www.ncbi.nlm.nih.gov/pubmed/29945233

ADD REPLY
0
Entering edit mode

That is what we ended up going with. We are just looking to modify it now to support 3 other variant callers.

ADD REPLY
1
Entering edit mode

Hi! Have you used ATpoint's solution?

ADD REPLY
1
Entering edit mode

We are still trying to get it working, right now we are working to rewrite it into our CWL ecosystem.

ADD REPLY
0
Entering edit mode

As a side note to this, anyone interested in contributing to this project? i.e. wrapping in CWL. We are a group from Mount Sinai Icahn Medical school and are building a solid repo of CWL wrapped tools with some consistent standards. https://gitlab.com/iidsgt/biowdl

ADD REPLY
2
Entering edit mode
2.1 years ago
djevanclark ▴ 40

We ended up writing our own solution which can be pulled from here. https://hub.docker.com/r/sinaiiidgst/consensus_caller

It's simply invoked by running consensus_caller from bash.

ADD COMMENT

Login before adding your answer.

Traffic: 3290 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6