Question: GATK DP filter for multiple samples
0
gravatar for sun.nation
2.1 years ago by
sun.nation120
United States
sun.nation120 wrote:

I used GATK to call SNPs and filter based on DP < 10. I did get SNPs but having problems to get filtered SNPs.

Command to filter SNPs:

$java -jar $gatk \
   -T VariantFiltration \
   -R Corca1_AssemblyScaffolds.fa \
   -o snpsDPlt10.vcf \
   --variant snps.vcf \
   --genotypeFilterExpression "DP < 10" \
   --genotypeFilterName "DP-lt10" \
   --setFilteredGtToNocall

$java -jar $gatk \
   -T SelectVariants \
   -R Corca1_AssemblyScaffolds.fa \
   -V snpsDPlt10.vcf  \
   -o snpsDPlt10-noCall.vcf \
   --setFilteredGtToNocall

$java -jar $gatk \
   -T SelectVariants \
   -R Corca1_AssemblyScaffolds.fa \
   -V snpsDPlt10-noCall.vcf  \
   -o snpsDPlt10-final-lt10DP.vcf \
   --excludeNonVariants

Final VCF file has annotated SNPs, how can I get SNPs with DP>10:

#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  cotton3.fastq_1.paired.trimmed  cotton4.fastq_1.paired.trimmed
scaffold_1      67871   .       G       A       485.77  PASS    AC=1;AF=1.00;AN=1;DP=13;Dels=0.00;FS=0.000;HaplotypeScore=0.0000;MLEAC=2;MLEAF=1.00;MQ=
60.00;MQ0=0;QD=32.77;SOR=1.609  GT:AD:DP:FT:GQ:PL       ./.:0,3:3:DP-lt10:99:118,0      1:0,10:10:PASS:99:396,0
scaffold_1      67878   .       T       A       485.77  PASS    AC=1;AF=1.00;AN=1;DP=13;Dels=0.00;FS=0.000;HaplotypeScore=0.0000;MLEAC=2;MLEAF=1.00;MQ=
60.00;MQ0=0;QD=30.77;SOR=1.609  GT:AD:DP:FT:GQ:PL       ./.:0,3:3:DP-lt10:99:122,0      1:0,10:10:PASS:99:392,0
snps gatk • 1.2k views
ADD COMMENTlink written 2.1 years ago by sun.nation120
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: 930 users visited in the last hour