This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!
Filtering does not work due to missing MQ values
I'm trying to filter my VCF on MQ value. I've done it for other annotations using VariantFiltration with no problems, but MQ is just not collaborating.
So first I had it in an expression with the other filters, and realized that the filter on MQ was not working when checking the output as there were lots of call with values that should have been filtered.
Then I did it individually (--filterExpression "MQ<40.0") and it starts running fine until it collapses when it encounters a call with missing data for MQ. Error message:
Message: for input string "nan"
I've checked with grep and indeed there're some missing MQ values (why?) but it seems to me the tool should be able to ignore these...? I can't seem to find an option to tell it to ignore them - is there a way to do it?
I also tried with SelectVariants and got the exact same error.
I'd appreciate any help you can give me here !
PS - I'm on GATK 3.7