Suggest vcf filtered parameter for gatk VariantFiltration
Using GATK we generated a vcf file and now using the following command, we have been trying to filter those reads. However, we observed little to no change in number of variants. I am unable to understand the problem Is there any error in the script below? Could you kindly suggest any other script or parameters to extract high quality snps
gatk VariantFiltration -R ref.fna -V raw.vcf -O filtered.vcf --genotype-filter-expression "QD < 2.0 || FS > 60.0 || MQ < 40.0 || MQRankSum < -12.5 || ReadPosRankSum < -8.0" --genotype-filter-name "my_filters"
• 3,725 views
•
link
updated
by
kanika.151
•
written
by
Mbillah
•
written
by
Mbillah
1 answer
Hello,
These are not genotype filters. The correct tag would be: --filter-expression "FS > 60.0" --filter-name "FSgreaterthan60"
Please refer to this discussion:
gatk_variant_filtration_issue
• 0 views
•
link
written
by
kanika.151
Log in to answer this question.
More posts like this
-
gatk Hardfilter Error
written by bestone •
-
No Passing Variants GATK
written by brandnewatthis •
-
Critical criteria for filtering variants by VariantFiltration
written by Alireza •
-
GATK 3.8.0 VariantFiltration ERROR input string
written by rturba •
-
VariantFiltration warning "undefined variable"
written by felipead66
-
GATK filtering console error
written by rbioinfo •
-
GATK VariantFiltration invalid expression
written by ciemanek
-
Question about GATK hardfiltering conditions
written by bitpir
-
Applying hard filters for variants
written by bioinforesearchquestions
-
Data_Type_Specific_Filters For Gatk Variantfiltration?
written by interfect •
