Gatk Indel Realignment Taking Forever! Help
Hi all, I am working on exome capture data for barley (1.3Gbp). I am interested in variant calling to find out SNPs in my sample. I have used SAMTools SNP calling and things get done in ~1 hr whereas...
View ArticleFrequency Of Homozygous Indels Vs Heterozygous Indels
I have been surprised to read in Genetic Variation in an Individual Human Exome that most of indels in coding regions are homozygous:The HuRef genome contains a total of 739 coding indels, which...
View ArticleDifferences Running Dindel On Different Unix Environments
Hi all, I am hoping to try out dindel for indel-calling on my NGS data but coming up against some problems getting it to work on both institute's computing cluster and my desktop machine. On the...
View ArticleGatk Indel Realignment Dies In A Memory Explosion... Any Fixes?
I got the following rather strange error while running GATK indel realignment on a considerably larger file (HiSeq) than I am used to:Exception in thread "main" java.lang.OutOfMemoryError: GC overhead...
View ArticleSamtools Indels -- Filtering Only Hits With Insertion In The Reference In The...
Can anyone suggest how to use samtools to filter only hits where there is a insertion in the reference that splits the sequence hit roughly by the middle? My sequences are in the range of 100-1000bp...
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 ArticleTool: Varscan: 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 ArticleClik here to view.
Indel Left/Right Alignment
When matching indels between different VCF files (generated by different callers), there is this issue with left / right indel alignment, For example: Here is a real example for one indel (from the...
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 ArticleExtracting Two Bam Files And Indexing Them To Compare Indels
I need to fetch 2 BAM files from an HTTP server and I only need 'chr1' so I'm attempting to use:samtools view -h http://cdna.eva.mpg.de/denisova/alignments/T_hg19_1000g.bam chr1 chr1.bamfor each file...
View ArticleChip-Seq Peak Caller Able To Find Peaks Split By Indels?
Are there any chip-seq peak callers able to find cases where the region under a peak has been split by an indel between the reference and the chipped sample?For example, let's suppose one has done...
View ArticleHow To Calculate Genotype Concordance Between Indel Call Sets?
How do people calculate genotype concordance between INDEL call sets? I have a NGS and BAC based INDEL call set (in vcf) but I get a higher false positive and false negative rate than mentioned in some...
View ArticleMerging Vcf Files And Filling In Missing Genotypes
I have vcf files with genomic variants (e.g. SNPs and Indels) from several patients. I now want to merge them into one vcf file e.g. using vcftools' merge-vcf. Doing that I get a file that holds for...
View ArticleAlternate Sources Of Reference Indel Vcfs Outside Of Dbsnp And 1000 Genomes...
Hi,I am in a midst of testing some software tools for detecting indel from single sample genome data, and was wondering whether there is any source of reference indel files outside of dbSNP and 1000...
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 ArticleTool: Phylosim - Monte Carlo Simulation Of Sequence Evolution In The R...
PhyloSimPhyloSim is an extensible object-oriented framework for the Monte Carlo simulation of sequence evolution written in 100 percent R. It is built on the top of the R.oo and ape packages and uses...
View ArticleSequence Variation Detection For Prokaryotic Genomes
I have resequenced a bacterial genome using SOLiD paired end reads. After looking at SNPs and small Indels, I'm now interested in large Insertions/Deletions (Size of genes and gene clusters).Which...
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 ArticleWhat Is The Best Tool For Mouse (Mm9 Or Mm10) Variant Annotations?
I am working with some mouse exome capture data, and I would like to provide functional annotations to variants (SNPs and small indels).For human data, I would typically use ANNOVAR. It looks like...
View ArticleBwa Indels
Comparing bwa 0.6.2 to bwa 0.5.9 I'm seeing very different behavior around indels. I'm using default options in both cases and as far as I can tell they haven't changed. I'm aligning Illumina paired...
View Article