This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Individual or joint variant call ing

Hello

I need a suggestion. I have exome data of 6 subjects. As the sample size is very less, should we do variant calling on every sample (https://gatk.broadinstitute.org/hc/en-us/articles/360037226672-CNNScoreVariants) rather than joint variant calling?.

P.S These 6 subjects belong to same generation or a different generation in a family with common ancestor.

Thanks

Priya

exome filtering analysis variants
• 977 views • link
updated by brunobsouzaa brunobsouzaa • written by priya.bmg • priya.bmg

Personally, I prefer performing variant call on every sample. But doing joint genotyping is also valid... So, see what best fit your purposes!

• 0 views • link
written by brunobsouzaa brunobsouzaa

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Sample Level Filtering During GATK Germline Short Variant Discovery
    written by jon.klonowski jon.klonowski
    I was looking over GATK4's pipeline for Germline Short Variant Discovery (https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels-) and what was lost on me is where sample level filtering occurs. I …
  • Variant calling for matched samples
    written by shivangi.agarwal800 shivangi.agarwal800
    Hi Guys, I want to call variants for 36 samples having matched normal. Can you suggest me how to run (what are the steps) GATK …
  • What is GenotypeGVCFs?
    written by wormball • wormball
    Hello! This article https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels- says i should use HaplotypeCaller in GVCF mode and GenotypeGVCFs then, and this article https://gatk.broadinstitute.org/hc/en-us/articles/360035531192-RNAseq-short-variant-discovery-SNPs-Indels- advises to use HaplotypeCaller without GenotypeGVCFs. …
  • Using GATK in conda
    written by michael.flower.14 • michael.flower.14
    I'm new to this, but I've successfully dabbled with alignment and variant calling. I installed conda as per the biostar handbook. It didn't have GATK, …
  • no variant remained after GenotypeGVCFs gatk4
    written by Mehmet Mehmet
    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 …
  • Contig Ploidy Model
    written by akshitiz09 • akshitiz09
    Hi everyone, I have some WGS data for humans to process and I'm using gatk for CNV calling. [https://gatk.broadinstitute.org/hc/en-us/articles/360035531152--How-to-Call-common-and-rare-germline-copy-number-variant] I'm running DetermineGermlineContigPloidy in CASE MODE. …
  • monozygotic twins ped file specification
    written by cocchi.e89 cocchi.e89
    Hi, I got some monozygotic twins samples (and a technical replicate of the same sample sequences twice) and I am running a `GATK GenotypeGVCFs` on …
  • GATK GenomicDBImport missing argument --intervals for Illumina HiSeq Exome??
    written by yussab • yussab
    Hi everyone, I'm using the GATK germline short variant discovery (SNPs + Indels) pipeline. https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels- After using HaplotypeCaller to call the GVCF, I'm not consolidating …
  • GATK somatic CNV results
    written by cocchi.e89 cocchi.e89
    I am using [GATK somatic-CNV pipeline][1] for some hg38 human sample (normal-tumor matched). In CollectAllelicCounts interval_list ([here][2] and [here][3] the discussions) I have tried to …
  • Mutect2 somatic hemizygote variant on chrX
    written by cocchi.e89 cocchi.e89
    I am working on calling some somatic variants from WES samples, I followed step-by-step the [GATK guidelines][1]. This is not a tumor sample so I …