Suggest vcf filtered parameter for gatk VariantFiltration
1
0
Entering edit mode
2.5 years ago
Mbillah ▴ 120

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"
gatk4 snps • 2.0k views
ADD COMMENT
1
Entering edit mode
8 months ago
kanika.151 ▴ 120

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

ADD COMMENT

Login before adding your answer.

Traffic: 1747 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6