Question: multi-sample VCF file filtering using 'PASS' for multi-sample VCF files.
0
gravatar for fbsja
2.3 years ago by
fbsja0
United Kingdom
fbsja0 wrote:

hiya,

Wanted to enquire about the filtering using the 'PASS' tag under the FILTER column in the VCF file for multi-sample VCF files.

I used the VariantFiltration GATK command below to generate different tags including the PASS filter.

########################################################################################
java -XX:-DoEscapeAnalysis -jar /home/software/GATK-3.5-0/GenomeAnalysisTK.jar \
-R hg19_withouM.fa \
-T VariantFiltration \
--variant /home/john/CFTR/multi_sample.vcf \
-o /home/john/CFTR/multi_sample.filtered.vcf \
--clusterWindowSize 10 \
--filterExpression "MQ0 >= 4 && ((MQ0 / (1.0 * DP)) > 0.1)" \
--filterName "HARD_TO_VALIDATE" \
--filterExpression "DP < 5 " \
--filterName "LowCoverage" \
--filterExpression "QUAL < 30.0 " \
--filterName "VeryLowQual" \
--filterExpression "QUAL > 30.0 && QUAL < 50.0 " \
--filterName "LowQual" \
--filterExpression "QD < 1.5 " \
--filterName "LowQD" \
--filterExpression "SB > -10.0 " \
--filterName "StrandBias"
########################################################################################

Does the PASS filter - correspond only to the first samples or all samples? I was planning to filter for just PASS sites as a general rule of thumb but this drastically reduces my total number of SNPs.

Wondering if I'd be better off keeping all variants and apply general QC filters like genotype rate? maf? individual missingness?

thank you very much.

cheers, J

sequencing multi-sample vcf • 1.2k views
ADD COMMENTlink modified 2.3 years ago by Pierre Lindenbaum117k • written 2.3 years ago by fbsja0

Does the PASS filter - correspond only to the first samples or all samples?

ALL samples.

for sample/ genotypes, use

genotypeFilterExpression and genotypeFilterName https://software.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_gatk_tools_walkers_filters_VariantFiltration.php#--genotypeFilterExpression

ADD REPLYlink modified 2.3 years ago • written 2.3 years ago by Pierre Lindenbaum117k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2276 users visited in the last hour