GATK filterExpression multiplication (*) mistaken for wildcard

[tsecogen]

Hi. I am using GATK for SNP identification and I have come across a problem when filtering the SNPs, where the * for multiplication (see code below) is being mistaken for a wildcard for all files in the folder. Does anyone have a solution for the * to be recognized as a multiplication and not a wildcard? Thank you!

System: Unix higher performance cluster

Code:

java -Xmx16g -jar GenomeAnalysisTK.jar \
-T VariantFiltration \
-R MtGenome_sections.fasta \
--variant GATK_4_Indel_filtered_Q30.vcf \
-o GATK_5_analysis_ready_Q30.vcf \
--clusterWindowSize 10 \
--filterExpression "MQ0 >= 4 && ((MQ0 / (1.0 * DP)) > 0.1)" \
--filterName "HARD_TO_VALIDATE" \
--filterExpression "SB >= -1.0" \
--filterName "StrandBiasFilter" \
--filterExpression "QUAL < 10" \
--filterName "QualFilter" \
--filterExpression "QUAL < 30.0 || QD < 5.0 || HRun > 5 || SB > -0.10" \
--filterName GATKStandard

[davidblaney]

Hi,

Have you tried backslash before it:

(1.0 \* DP)

or using single quotes around the expression?

[tsecogen]

Thanks for the suggestions but unfortunately neither of them work. I still had the same problem with * being recognized as a wild card.

[mastal]

GATK filterExpression multiplication (*) mistaken for wildcard

Have you tried MQ0 / DP > 0.1?

Obviously this will only work for cases when you want to multiply DP by 1.0.

[tsecogen]

Thanks for the thought, but DP (Depth of Coverage) changes for every SNP so it is not safe to assume it is 1 or any other number for all SNPs.

In the end I just transfered the input files to my imac to run that piece of the code. Not the most elegant work around but at least it allows me to move forward.

Thanks.