Critical criteria for filtering variants by VariantFiltration
0
0
Entering edit mode
21 months ago
Alireza ▴ 30

Hi all.

I run the GATK VariantFiltration with the following parameters (according to GATK recommendation) to find robust variants. Then I used them as input for annotating variants.

Do you have any suggestions for better filtration?

Is it recommended to run VariantFiltration on different variants (somatic, germline, CNV)?

For SNP:

gatk VariantFiltration \
        -R genome.fa \
        -V germline_raw_SNP.vcf \
        -O germline_filtered_SNP.vcf \
            -filter-name "QD_filter" -filter "QD < 2.0" \
        -filter-name "FS_filter" -filter "FS > 60.0" \
        -filter-name "MQ_filter" -filter "MQ < 40.0" \
        -filter-name "SOR_filter" -filter "SOR > 4.0" \
        -filter-name "MQRankSum_filter" -filter "MQRankSum < -12.5" \
        -filter-name "ReadPosRankSum_filter" -filter "ReadPosRankSum < -8.0" \
        -genotype-filter-expression "DP < 10" \
        -genotype-filter-name "DP_filter" \
        -genotype-filter-expression "GQ < 10" \
        -genotype-filter-name "GQ_filter"

For INDELS:

gatk VariantFiltration \
        -R genome.fa \
        -V germline_raw_INDEL.vcf \
        -O germline_filtered_INDEL.vcf \
            -filter-name "QD_filter" -filter "QD < 2.0" \
        -filter-name "FS_filter" -filter "FS > 200.0" \
        -filter-name "SOR_filter" -filter "SOR > 10.0" \
        -genotype-filter-expression "DP < 10" \
        -genotype-filter-name "DP_filter" \
        -genotype-filter-expression "GQ < 10" \
        -genotype-filter-name "GQ_filter"
WGS Filtration gatk WES • 576 views
ADD COMMENT

Login before adding your answer.

Traffic: 816 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