Questions tagged [gatk]
Genome Analysis Toolkit. Popular open-source bioinformatics suite to call variants from second-generation sequencing data.
59
questions
1
vote
1
answer
25
views
I am trying to create a subset of 10k variants from 25-30 unmapped contigs of a g.vcf file including the header
My objective is to take a g.vcf.gz file and from 25-30 unmapped contigs with titles like "NW_020192317.1", I want to make a subset of ~10k variants from each of the unmapped contigs and make ...
4
votes
0
answers
49
views
MergeBamAlignment error
I doing the alignment of samples following the GATK pipeline, and doing the MergeBamAlignment,like this:
...
0
votes
1
answer
73
views
A bash script for running on a bunch of bam files
I have some bam files in this directory
/data/Continuum/WES/results/
I want to run GATK mutation calling over bam files
I ...
1
vote
1
answer
45
views
Importing local files to a GATK docker image
I have pulled GATK docker on my computer but when running a GATK command, can not locate my input file
My OS version
...
1
vote
0
answers
14
views
GenomicsDB and recalibration recommended for cohort of mixed exome (GATK)
We have a cohort of around 3k exomes, most of them from different capture kits, such as Agilent v6, v8, Nextera, Twist, Clinical Research exome and many others. Right now we are creating genome db on ...
1
vote
1
answer
104
views
Getting VCF file that contain common SNPs from 6 VCF file using isec
I have 6 VCF files, where I would like to obtain the SNPs that are common (by position) in all the 6 files. I have tried this command
...
2
votes
1
answer
51
views
WES variant calling with DNBSEQ-T7: technical quality assessment
I recently received whole exome sequencing samples that were sequenced on an MGI sequencing instrument, DNBSEQ-T7. I am interested in somatic variant-calling on the paired tumor-normal samples. As ...
1
vote
1
answer
75
views
GATK SplitIntervals: A USER ERROR has occurred: but no positional argument is defined for this tool
I am trying to create a GATK g CNV pipeline and I need to make shards before running the GermlineCNVCaller therefore, I believe I need to use ...
2
votes
1
answer
180
views
Singularity & GATK 4.4.0.0. CreateSequenceDictionary
I ran the following code in bash (I added print statements for ease):
...
-1
votes
1
answer
51
views
GATK g CNV code Aborts Before Running Main Code
I have been struggling to get the code I wrote for running GATK g CNV on my Linux machine (Steps Here: https://gatk.broadinstitute.org/hc/en-us/articles/360035531152 ) to work. Specifically it seems ...
0
votes
0
answers
71
views
Combine gvcf for whole population or for different phenotypic groups?
I created 400 or so gvcf files (g.vcf.gz) using the gatk HaplotypeCaller. These files show populations for three distinct traits (HR, MR, HS).
Now, I want to use gatk GenomicDBImport to combine gvcfs ...
2
votes
2
answers
153
views
What should GC coverage bias plot of exome data look like?
I want to see if there is any GC coverage bias in some paired-end Illumina human exome data. Update This is for a variant calling project. I used the ...
2
votes
1
answer
90
views
Is it common to get different number of SNVs+Indels across samples from vcf files generated using GATK and DRAGEN (counts are higher for GATK)?
We have 2 vcf (Whole Exome Sequencing (WES) data; germline samples) files (e.g., vcf_1 and vcf_2). vcf_1 was generated (Ref. genome: hg38) using the GATK pipeline for 250 children and their parents ...
4
votes
0
answers
111
views
Write a bash script to run gatk, fix errors with input, and rerun until completion
I have a bam file that I want to run through GATK's SplitNCigarReads tool. Because of the way the bam file was generated, the program will often fail, with an error message stating:
...
1
vote
0
answers
47
views
Weird data in VCF file
I am using GATK HaploTypeCaller for SNP calling. I have 300 samples. As far as I understand that if a sample is missing a SNP it should be denoted by ./., however, I am getting a lot of values like 0/...
0
votes
0
answers
50
views
What is the best tools to find all the somatic mutations?
I am a beginner in this field. I would like to know which is the best way to get all the mutations from a certain sample. At this point, I am considering using GATK (here) and Maftools (here). I don't ...
3
votes
1
answer
212
views
Funcotator reference file error in GATK4
I am using a function called Funcotator (here) in GATK4. The code I am trying to run is:
...
6
votes
1
answer
885
views
Why does GATK produce both 0/1 and 1/0 genotypes in the same file? Are the two not equivalent?
I have always thought that 1/0 and 0/1 in VCF genotype fields are equivalent. And yet, GATK uses both. For example, these are ...
1
vote
1
answer
552
views
GATK variant filtration using "SelectVariants" and use of JEXL queries
How to do variants selection in some corner cases using GATK and JEXL expressions?
I am following the guidelines given in this links for variant selection for some specific cases are not as much ...
2
votes
2
answers
4k
views
How to convert a BAM file into a VCF file
I know that there is a lot of methods to do so, but most of them need to have a reference file. (For example, the GATK force me to provide a reference.)
However, I think the BAM file I am working with ...
-1
votes
2
answers
233
views
Error in running Mutect2
I am trying to run Mutect2 on a .bam file but I get an error by Googling can not be tackled
Have you seen this error?
...
-4
votes
1
answer
223
views
Getting error in running GATK
I am trying to making a fasta doc using GATK tools but I don't know what does this error say
Server2:/data2/RNASeq/Angel/gatk$ ./gatk CreateSequenceDictionary -R reference.fa
Using GATK jar /data2/...
2
votes
1
answer
226
views
MarkDuplicatesSpark failing with cryptic error message. MarkDuplicates succeeds
I have been trying to follow the GATK Best Practice Workflow for 'Data pre-processing for variant discovery' (https://gatk.broadinstitute.org/hc/en-us/articles/360035535912).
This has all been run on ...
2
votes
0
answers
49
views
About workflow input section in DRAGEN-GATK v1 whole genome analysis script
I am trying to understand the WDL script for the DRAGEN-GATK whole genome analysis and I am faced with this input section inside the workflow definition which I really don't get:
...
0
votes
1
answer
163
views
Convert VCF file to mpileup.txt
I am working on an iterative analysis that uses orthologous pipelines that require mpileup.txt files as input for a visualization step. This requires me to convert VCF files to mpileup.txt.
This ...
1
vote
1
answer
224
views
Dictionary and index of vcf for base recalibration step
I need suggestions in creating index and dictionary with vcf files. For the base recalibration step, I downloaded Homo_sapiens_assembly38.known_indels.vcf.gz from the given link: https://console.cloud....
2
votes
2
answers
743
views
snakemake multiple parameters for multiple input and single output in snakemake. ConbineGVCFs gatk problem
I have written a rule for CombineGVCFs in gatk4. The rule is as follow
...
0
votes
1
answer
208
views
Best practice for running GATK VQSR on X chromosome
According to GATK best practice, it is recommended that different VQSR models be built for SNPs and INDELs, because the annotations for high-quality SNPs and INDELs are systematically different (if I ...
2
votes
1
answer
84
views
Low pass sequencing has been reported to detect common variants. How low can one go and get reliable data? Is 2X pass sequencing analysis possible?
I would like to use low pass sequencing to replace a genotyping chip to be able to detect variants up to 0.1 % allele frequency in available population data. What is the minimum depth I can opt for to ...
1
vote
1
answer
490
views
GATK 4.2 using IndexFeatureFile with multiple files
I am trying to use GATK IndexFeatureFile to create ".vcf.gz.tbi" file. But I have 480 files.
For the pilot study I used:
...
2
votes
0
answers
41
views
Joint variant calling of two sample sets with different coverages (~5x vs ~25x)
Let's say 50 samples were sequenced at 25x. After a while, the group realizes that the experiment is under-powered and adds 100 samples more, but at a coverage of 5x. I'm talking about WGS.
The ...
0
votes
0
answers
41
views
How to control/normalize for number of reads when calling SNPs using RNA-Seq?
I used the GATK pipeline to call SNPs on males and females using RNA-Seq data. But the males have a higher read count (~43-46M reads) than the females (~40-42M reads). This causes SNP counts to be ...
2
votes
1
answer
41
views
Cromwell: can't find file ID
I am using Cromwell to test a CWL (common workflow language) implementation of the GATK somatic CNV (copy number variation) pipeline. I am currently getting the following error (full log at the end):
...
1
vote
1
answer
366
views
No MQ tags in VCF files
To call minority variants in my Mtb sequences I'm using a pipeline of
...
1
vote
1
answer
678
views
How to generate the json files for Cromwell workflow execution?
I am starting to use Cromwell WDL scripts for running GATK workflows on a slurm cluster. When using the paired-fastq-to-unmapped-bam script I need to define the fastq sepcifics in the appropriate json ...
2
votes
2
answers
1k
views
Marking optical or PCR duplicates with picard vs. samtools flagstat
I am trying to identify the best library prep method for noninvasive prenatal test samples, to be sequenced with Illumina Novaseq.
One metric that I am evaluating is the number/ percentage of PCR or ...
0
votes
1
answer
436
views
What samples can be used for a Mutect2 Panel of Normals
I'm working on calling somatic variants from solid tissue tumors. I have one normal and one tumor sample for each patient. I plan to use Mutect2 to call somatic variants after preprocessing the data ...
2
votes
3
answers
1k
views
Get start and end coordinates per chromosome
I will use GATK for SNP calling (HaplotypeCaller). I need to feed the interval file in the command, otherwise I get errors (even though I want to use the whole genome, not a subset - btw it's not the ...
0
votes
1
answer
352
views
Interpreting MultiQC plot of GATK BaseRecalibration data?
I'm using GATK4 to build a somatic variant calling pipeline. The pipeline uses MultiQC to aggregate quality control data, and one of the QC measures reported is base quality score recalibration from ...
2
votes
0
answers
47
views
NA12878 equivalent for somatic variant calling [duplicate]
I'm building a somatic variant calling pipeline with snakemake. The pipeline is hosted on Github, and I'd like to use Travis CI for continuous integration testing. ...
1
vote
1
answer
347
views
GATK 4.1.4.0 Mutect stats output
I have adapted my pipeline to the new filtering strategies implemented in gatk 4.1.4.0 and while looking at its output I noticed that the stat file generated by mutect2 has a negative number:
...
1
vote
0
answers
595
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
1
answer
311
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 ...
2
votes
1
answer
369
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
2
answers
170
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
2
answers
372
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
1
answer
138
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 ...
2
votes
0
answers
153
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
2
answers
1k
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
1
answer
4k
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
...