Using -"-allow-old-rms-mapping-quality-annotation-data" to force genotyping gvcf created with old version
AnsweredHello,
I created sample gvcfs using GATK v.3 and want to joint genotype those files.
I used GATK v.4 to create a Genomics DB using GenomicsDBimport and it worked well.
But when using GenotypeGVCFs to genotype the database, I got the error:
A USER ERROR has occurred: Bad input: Presence of '-RAW_MQ' annotation is detected. This GATK version expects key RAW_MQandDP with a tuple of sum of squared MQ values and total reads over variant genotypes as the value. This could indicate that the provided input was produced with an older version of GATK. Use the argument '--allow-old-rms-mapping-quality-annotation-data' to override and attempt the deprecated MQ calculation. There may be differences in how newer GATK versions calculate DP and MQ that may result in worse MQ results. Use at your own risk
I thus tried to use GenotypeGVCFs with the argument '--allow-old-rms-mapping-quality-annotation-data' and it worked.
Nonetheless, I am wondring about the impact of using this argument and I could not find any information about it.
Is using this argument making GATK v.4 using the same way to calculate MQ as GATK v.3?
Does this argument amount to using GATK v.3 for all other calculations of statistics or solely for the calculation of MQ?
For the downstream filtering, should I use GATK v.3 or v.4?
Thanks!
If not an error, choose a category for your question(REQUIRED):
a)How do I (......)?
b) What does (......) mean?
c) Why do I see (......)?
d) Where do I find (......)?
e) Will (......) be in future releases?
-
There have been a few previous discussions about this argument on the forum. Please take a look and if you have more questions let me know.
- https://gatk.broadinstitute.org/hc/en-us/community/posts/360074172971-Missing-MQ-scores-in-gvcf
- https://gatk.broadinstitute.org/hc/en-us/community/posts/360075376311-GenotypeGVCF-by-intervals
Best,
Genevieve
Please sign in to leave a comment.
1 comment