no variant remained after GenotypeGVCFs gatk4
0
0
Entering edit mode
16 months ago
Mehmet ▴ 780

Dear all,

I am analyzing a DNA-seq data to identify variants in two genes of human using a single-sample based gatk4 pipeline (https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels-).

Haplotypecaller identified known variants that start with rs such as rs11490246. Later, I applied GenotypeGVCFs.

After applying GenotypeGVCFs, there was no known variants that were previously found at HaplotypeCaller step in output of GenotypeGVCFs.

Here is my all pipeline after generating analysis-ready bam files. This is for one bam of a single sample. I don't have multiple bam files from the same sample, which is why I followed single-sample approach at the website (https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels-).

Variant calling

HaplotypeCaller

gatk --java-options "-Xmx5g" HaplotypeCaller -R hg19.fa -I 21050.AnalyseReady.bam  -O 21050.HaplotypeCaller.output.vcf.gz -D dbsnp_138.hg19.vcf -ERC GVCF

GenotypeGVCFs

gatk --java-options "-Xmx4g" GenotypeGVCFs -R hg19.fa -V 21050.HaplotypeCaller.output.vcf.gz -O 21050.HaplotypeCaller.output.g.vcf.gz.GenotypeGVCFs.output.vcf.gz --tmp-dir=./Temp_files

CNNScoreVariants

gatk CNNScoreVariants -V 21050.HaplotypeCaller.output.g.vcf.gz.GenotypeGVCFs.output.vcf.gz  -R hg19.fa -O 21050.HaplotypeCaller.output.NEW.vcf.gz.CNNScoreVariants.OUT.annotated.vcf 

FilterVariantTranches

gatk --java-options "-Xmx7g" FilterVariantTranches -V 21050.HaplotypeCaller.output.NEW.vcf.gz.CNNScoreVariants.OUT.annotated.vcf  --resource Mills_and_1000G_gold_standard.indels.hg19.sites.vcf.EDITED.gz --info-key CNN_1D --snp-tranche 99.95 --indel-tranche 99.4 -O 21050.FilterVariantTranches.output.cnns.cnnfilter.vcf

Am I missing something? or Should I skip GenotypeGVCFs? I could not make sure about the GenotypeGVCFs step in single-sample variant call pipeline.

GenotypeGVCFs gatk4 variant • 447 views
ADD COMMENT

Login before adding your answer.

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