HaplotypeCaller GVCF on one bam files · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

HaplotypeCaller GVCF on one bam files

Hi Guys, I have two fastq files for a single patient. They have paired-end reads. I am using Haplotype caller to generate variants. So while generating the variants should I Just generate vcf directly from haplotype caller or a gvcf and then use genotype calling. From my understanding gvcf is effective while calling variants on a cohort of bam files. Is that correct?

Thank you

variant calling gatk haplotypecaller

• 1,951 views • link

written 5.7 years ago by rayzghimire •

gvcf is effective while calling variants on a cohort of bam files. Is that correct?

yes

• 0 views • link

written 5.7 years ago by Pierre Lindenbaum

Hi Pierre, So should I use genotypeCalling for a single gvcf file? How is it different from the vcf generated by HaplotypeCaller?

Thank you

• 0 views • link

written 5.7 years ago by rayzghimire •

you just call HaplotypeCaller without -ERC GVCF and you'll get a standard vcf

• 0 views • link

written 5.7 years ago by Pierre Lindenbaum

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Is this the correct order of steps in GATK germline cohort variant calling (VQSR workflow)?
written 10 months ago by Ramnaresh •

I’m performing germline variant calling using GATK (gVCF-based workflow) and would like to confirm whether the order of my steps is correct. What I have …
Making consensus sequence for each haplotype
written 4.7 years ago by michael.flower.14 •

I'm dealing with paired end amplicon sequencing data. I've produced a GVCF file with haplotype calls using: gatk HaplotypeCaller -R $REF -I "$BAM" -O "$OUT"/results/variants/${SN}_HaplotypeCallerPGT.vcf …
How to convert gvcf file to normal vcf file ?
written 5.4 years ago by DareDevil

Hi, I followed following steps to arrive at gvcf file: #Mark duplicates in the aligned file gatk MarkDuplicates -I s1.bam --METRICS_FILE s1_markdups_gatk_metrics.txt -O s1_sorted_gatkMarkDups.bam --VALIDATION_STRINGENCY …
bcftools - creating indels
written 5.8 years ago by rayzghimire •

Hi, So I am converting fastq files into vcf. I wanted to know while calling variants how to include both SNPs and indels using samtools …
[gatk 3.6.0] GVCF - haplotypecaller - cross referencing multiple VCF
written 6.1 years ago by Jalil Sharif •

Hi there, how do you make a GVCF file from bam files, using [gatk 3.6] haplotypecaller if wanting to compare patient with their two parents? …
Hard filtering vcf files
written 6.7 years ago by Kash

Hi all, I am trying to find a set of variants for a non-model organism that can be used as a known set (well-curated training/truth …
Query regarding somatic variant calling
written 7.0 years ago by iti.gupta •

Hello, I had a query regarding variant calling. I have three tissues from the same patient (Tumor, Non-Tumor, and Blood) I wanted to identify the …
GATK joint genotyping
written 8.4 years ago by win

Hi all, i am struggling a bit with preparing a cohort genome vcf file for joint genotyping using GATK. I tried with 30 BAMs from …
GATK - gene panel VCF - too many variants - how to filter out
written 9.3 years ago by stanedav •

Hello guys, I have question about filtering NGS data with GATK workflow. I have data from NGS Illumina (Sureselect kit), it is targeted gene panel …
Do freeBayes and platypus support gVCF (ie incremental variant calling for large sets of genomes)
written 11.8 years ago by William

Do freeBayes and platypus support gVCF? https://sites.google.com/site/gvcftools/home/about-gvcf With GATK you used to have the 1000 + 1 genome variant calling problem, where you had to …