The Stack Overflow podcast is back! Listen to an interview with our new CEO.

Questions tagged [gatk]

Genome Analysis Toolkit. Popular bioinformatics suite to call variants from second-generation data.

Filter by
Sorted by
Tagged with
1
vote
0answers
25 views

vcf tools Error: Require Genotypes in VCF file in order to output Fst statistics

I want to calculate Fst by vcf tools and GATK. I did this steps: for creating gvcf: gatk HaplotypeCaller -R ref.scf.fasta -I input.bam -ERC GVCF -O out.g.vcf ...
1
vote
1answer
35 views

Subset a multisample VCF file

I need some help with GATK. I have a multisample vcf file containing 18125 samples. After performing PCA I have around 303 samples as outliers and I wanted to remove these samples from my multisample ...
1
vote
0answers
46 views

How to use -csvStats option?

I'm trying to get CSV files from snpEff -csvStats option and got stuck. I tried the code below and only got the vcf file with no contents. How can I use -csvStats option properly? ...
2
votes
2answers
128 views

Variant calling without matched normal sample

I have WGS .bam files for 3 patients (tumour and its matched derived model namely organoid) but I don't matched normal sample. If I call variants of each patients (tumour and its matched organoid), ...
4
votes
2answers
77 views

Targeted NGS, up to 99% of reads have been marked as duplicates

Currently I'm performing whole analysis (pipeline from *.fastq to *.vcf) of 41 samples (targeted NGS). I rely on GATK best practices, however with some modifications. I decided to use the following ...
6
votes
1answer
71 views

GATK documentation for required depth to reliably call heterozygous mutation in diploid organism?

I'm looking for official GATK documentation (or a recent manuscript) that defines a general recommendation/requirement for sequencing depth to reliably call a heterozygous point mutation in a diploid ...
1
vote
0answers
116 views

Meaning of category in MutSigCV

I am trying to understand the working of MutSigCV. But I am not sure about some terms used in MAF file format which is category. According to MutSigCV, every mutation can be divided into the following ...
5
votes
2answers
517 views

No variant found using GATK 4.0 HaplotypeCaller

I am doing variant calling on RNA-seq datasets from wheat which is hexaploid,the binary alignment (BAM) files were created using STAR version 2.6.0c and variant calling was done using GATK 4.0 ...
3
votes
1answer
515 views

Increase number of threads for GATK 4.0 HaplotypeCaller

I am using GATK version 4.0, I tried to use multiple threads for calling variants using HaplotypeCaller using following command ...
3
votes
1answer
125 views

Human Genome Alignment and Variant Calling Benchmarks

I am doing a project that requires bench-marking the performance of computational tools for alignment and variant calling of human genome sequencing data. In particular I would like to establish ...
1
vote
0answers
45 views

Normalize by counting CNV within a pool of VCF files

I’m (trying) using the GATK4 germline CNV calling pipeline. I successfully got 57 VCFs from my sample batch, called with segments (obtained by merging the contiguous intervals), like in a classic VCF :...
4
votes
3answers
416 views

Issues performing variant calling with GATK

I am trying to perform variant calling on a BAM file generated through STAR version STAR_2.6.0b for wheat genome using GATK haplotypecaller as follows: ...
4
votes
2answers
320 views

GATK CombineVariants complains the contig order in the VCF files

I have called variants on two strains of C. elegans separately. I now want to merge the VCF files into one using the following code: Create a sequence dictionary of the reference sequence Sort the ...
6
votes
1answer
566 views

Merge 2 VCFs from different variant callers

I am working with WES data for detection of somatic variants and I have used two variant callers because no variant caller is complete in itself. I have used GATK Haplotypecaller for small variants ...
6
votes
2answers
5k views

Converting VCF file to PLINK bed/bim/fam files

I am trying to find the best way to convert VCF files to PLINK binary bed/bim/fam files, but it seems like there are many varied ways to do this. (For example, using Plink 1.9 --vcf tag, bcftools, ...
4
votes
1answer
161 views

Where can I get the population allele frequency vcf file?

I want to use GATK to estimate cross-sample contamination for Whole Genome Sequencing data. The specific tool is ContEst and it is run with: ...
5
votes
3answers
344 views

Is spark widely used in bioinformatics?

I learned that GATK 4 is using Spark for parallelization. I googled around, though I am still not quite sure how spark really works and how to use it in practice. Besides GATK 4, are any other ...
17
votes
2answers
1k views

Single-sample vs. joint genotyping

I am trying to understand the benefits of joint genotyping and would be grateful if someone could provide an argument (ideally mathematically) that would clearly demonstrate the benefit of joint vs. ...