Quantcast
Channel: Recent Discussions — GATK-Forum
Viewing all articles
Browse latest Browse all 12345

Combining RAW_MQ and MQ

$
0
0

Hi,

I have two batches of gvcfs :
- one was processed using GATK-3.3
- one was processed using GATK-3.5

I want to perform a joint calling + VQSR using both batches together. The problem is that the batch from GATK-3.3 has the MQ annotation while the other from GATK-3.5 has the RAW_MQ annotation.

Can the GenotypeGVCFs of GATK3.5 handle both MQ and RAW_MQ to compute a unique MQ value that can be used by VQSR ?

It seems that VariantAnnotator and CombineGVCFs cannot add the MQ annotation to the gvcfs from GATK-3.5 (using the -A parameter).

Best


Viewing all articles
Browse latest Browse all 12345

Trending Articles



<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>