Gatk indel realignment error - mismatch in index files and dict file
Greetings,I am aligning pooled sequencing data to a new renferece genome. GATK won't generate intervals because not every scaffold in the reference is found in my bam index?What am I missing? It seems...
View ArticleSomaticIndelDetector vs UnifiedGenotyper
Hey Guys, I've been establishing my own first pipeline for calling variants: Indels and SNPs. For that, I used mainly GATK tools. I'm doing this for whole genome and whole exome sequences of mouse...
View ArticleWhat do you do with Indels data?
I have a list of insertion and deletion from a whole genome sequence obtained by next generation sequencer, chr, pos and sequence for each indels aligned to our reference genome. I'm not familiar with...
View ArticleGATK3 vs Dindel for small indel calling
I read some posts from Heng Li around two years ago that he suggested Dindel is the optimal software for small indel calling. I'm just wondering if this is still true today because GATK seems to...
View ArticleVarScan: variant detection in massively parallel sequencing data
VarScan is a platform-independent mutation caller for targeted, exome, and whole-genome resequencing data generated on Illumina, SOLiD, IonTorrent, Roche/454, and similar instruments. The newest...
View ArticleHow to handle reads ending with deletions in GATK?
Hello,before asking my question, I should point out that I'm working with data that's not my own (publicly available), to learn and establish a proper workflow when real data wlll arrive in the...
View ArticleGATK SomaticIndelDetector VCF info field missing
Hey everyone,I am a bit confused. I called my INDELs pairwise using the SomaticIndelDetector. But in the resulting VCF the info field seems to be missing, is this normal? Or am I getting something...
View ArticleSNPs and Indels in homozygous genomes
I have a genomic sequence of an inbred mouse strain i.e. the genome is homozygous for a given position or same alleles are present on both of the strands. Now for variant calling, I have to compare it...
View ArticleTrio: detecting a large de novo indel
I've been given a set of three BAM (father, mother, child) and I expect the child to contains some de-novo heterozygous variations . samtools mpileup have been used to find the small variations but...
View ArticleSplit VCF file into SNPs and Indels
Hi there,As recommended in the GATK best practices the Variant Quality Score Recalibration has to be done separately for SNPs and Indels. But, I didn't find the way to do this split in a clean way (for...
View ArticleHow are INDEL coded in VCF files
Hello,I have a question concerning how are INDEL coded in VCF files. I have for instance the following file:CHROM POS ID REF ALT QUAL FILTER8 18078835 . TTA T 46 PASS8 18078836 . TA T 138 PASSHow to...
View ArticleExtracting specific reads from VCF file
Hello,I am trying to extract indels with coverage ~1 and they are in a VCF file. How do I write a python script or a shell commnad to extract those indels with coverage 1 only? This is the data:chr1...
View ArticleMerging/concatenating vcf files
I have a vcf file of SNPs and another vcf file for my Indels. During SNP calling step by GATK, I created them separately instead of using -glm (calling together SNPs and Indels). Now to get the...
View ArticleFiltering vcf file
I was wondering how to filter the vcf file based on a few input arguments ( DP>10, MQ>30 and QD>20 or GT = "1/1" etc)? I m planning to use simple command on the command line to extract the...
View ArticleLocating indels in gene
I have a excel file from windows (spreadsheet) contains start, end and indel sequences (-GGC means GGC deletion or +AAC means AAC insertion) in column 1, 2 and 3 respectively. I have a reference genome...
View Articlestrelka indel allele counts
Hopefully there are some Strelka users out there who can help with this one. I'm looking for allelic counts of the ref and non-ref allele of the indels.I see these flags in the VCF header, do these...
View Articlesearching for deletion in a low number of reads
Hi everybody,I am having sort of a non-conform problem.we have a mutated mitochondrial genome (circular) of high coverage (~x10k). We know that there are a few deletions and insertions in there due to...
View ArticleInsertion Detection tools
Hi,Please suggest me some tools which can detect large indels (esp. of size > 100bp). Some of the tools I have already tried include Breakdancer and Pindel and Hydra. Any other suggestions??
View ArticleHow to detect overlaping indels in VCF file?
Hi all,I have a vcf file with indels from a whole genome analysis and I want to detect overlapping between indels... I tried to use BEDtools-intersect but it asks for 2 files, and I only have a single...
View ArticleMulti-sample Indel Realignment using GATK
Is there a way to use the GATK Indel Realigner tool with a two-samples-input? The issue is that in my normal sample, several reads are mapped containing SNPs at specific positions, and ONE read having...
View Article
/fit-in/580x348/filters:format(webp)/english-betterindia/media/media_files/2026/09/17/balaghat-4-2026-09-17-19-01-31.png)
