Next Gen Sequencing Indel Realignment With The Genome Analysis Toolkit (Gatk)
As I continue to add steps to my SNP/Indel discovery workflow, the latest recommendation is local realignment around Indels using GATK following the initial alignment step. I have just commenced the...
View ArticleHow To Set Filter For Frequency Of Reads And Hapmap Exome Sample Results:
Hey All, I only used so far three filters for my whole exome pipeline (aligning to hg19) for a HapMap sample. I tried it on the NA19240 Hapmap sample from paper below (Table 3) which shows ~196...
View ArticleIndel Detection For 454 Resequencing
I do some work for a small diagnostics company that has a requirement for small indel detection in 454 data. Those of you familiar with Roche's pipeline will be aware that AVA (Amplicon Variant...
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 ArticleIndel Discovery Delly, Pindel, Samtools, Gatk
I'm testing out samtools vs GATK for snp and indel calling, and looking at using pindel for SV in particular focusing on insertions and Delly for the other SV. What experience do people have of SNP...
View ArticleFunctional Consequence of In-frame Indel
Does anybody include in-frame indels in association studies? A quick pubmed / google / biostars search on this didn't come up with much. In-frame indels are often treated as synonymous, but since they...
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 ArticleGenomic Alignment And Snp/Indel Calling - My First Ever "Pipeline"
So I have finally generated some reads and run it through what I guess could be called a very rudimentary 'pipeline'. I generated a million paired end reads with wgsim then aligned with bwa, and used...
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 ArticleWhole Genome Dna Amplification Before Applying Wgs
Dear allWe have some DNA samples that I would like to run whole genome sequencing on. Unfortunately these samples only contain very little DNA (too little for WGS). It has then been suggested to do...
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 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 ArticleWhat Is The Difference Between Samtools Mpileup And Pileup
As of samtools 0.1.8 'mpileup' appears to have superseded 'pileup' for variant detection. The samtools website now has pileup as deprecated. So i am wondering whether i should repeat my recent...
View ArticleStructural Variation Detection On A Short Reads Simulation
Hi,I am trying to do the following exercise (still playing with indels):simulate short reads (and indels) out of a human genome chromosome,align the short reads simulated in (1),call indels using a SV...
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 ArticleDetecting Cancer Somatic Mutations From Multi-Sample Low-Pass Wgs Data
Hi All,Q1. I'm trying to find recurrent mutations (SNPs and indels) from multi-sample (~100) low-pass (3~4x) whole genome sequencing data for a certain cancer type. They are all paired with matched...
View ArticleGatk Threshold Parameter For Clipping Of Reads During Indel Calling Running...
GATK 1.6 claims to apply a Q20 threshold to clip ends of the reads for the indel caller:http://www.broadinstitute.org/gatk/guide/article?id=12373. Indel Calling with the Unified Genotyper [...] while...
View ArticleDo You Guys Have Any Recommendation On Variant Calling Using Long Contigs Or...
I have assembled a human genome and used bwa mem aligning the contigs back to the reference genome. Do you guys have any recommendation of variant caller?
View ArticleSamtools / Bcftools Missing Random Obvious Single Base Indels
I am using samtools to call variants in a haploid genome (yes I know it is designed for diploid). It finds SNPs easily, and most of the indels, but there are a few indels it can NOT seem to find no...
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 Article
