Strange variant call it show a indel but looks like a SNP
Hello all, i have a variant from a ION PGM sequencing, called with the ION Torrent variant caller plugin CHR Â Â POS Â Â Â ID Â Â REF Â Â Â ALT chr22 Â 44324726Â . Â Â Â TÂ Â Â Â Â Â TG The flow cell...
View ArticleWhat Are The Next Steps For Processing Indels From Gatk
Hi, I have raw calls for indels called using GATK, I wonder, where to start? filtering these calls? Should I plot the quality, depth etc. and then decide threshold to filter the bad calls? Any...
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 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 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 ArticleMemory Allocation for VarScan
I am running the somatic variant caller in VarScan on several files. Â On average, the input .pileup files are ~50 GB each (with both tumor and normal samples).I've noticed that I've gotten 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 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 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 ArticleGatk 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 ArticleHow Best To Add Custom Variation Track To Ucsc Browser?
I have generated sequence variants (SNPs/INDELs) by NGS resequencing and would like to display them on the UCSC browser (I have resequenced against a UCSC reference genome). Obviously i could use BED...
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 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 Does Left-Alignment Work For Indels? Does It Address All Complicated...
I can vaguely guess its meaning but want to know, at least, how it works in different situations. Besides, after left alignment, can we write an indel in VCF file in a unique way?Thank you!
View ArticleHow to call low frequency snp/indel with GATK's HC/UG
i want to call snp/indel with low frequency, like DP=100, ALT=10.how to set HC/UG ?with UG, i run with :--min_indel_count_for_genotyping 2--min_indel_fraction_per_sample 0.05get a nice result, but i...
View Articlecalling 1000G SNPs/indels from ngs data (chipseq or other source)
I would like to call known (e.g. from 1000G project) SNPs and indels from sequence reads, in my case a chip-seq project with 20 samples, which i will later also use to query allelic imbalance. Is there...
View ArticleClassify Mutations Which Are Found Ny Using Exome Sequencing
Hi, everyone. I used exome sequencing to sequence tumor and blood, and after I remove mutations in blood and snp in dbsnp131, I found I still have about 2000 mutations. I read the exome sequencing...
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 ArticleIs there any way to quantify differences between neighbour-joining...
I am working on a project and I am currently only using the SNP data to make my trees. I would like to compare the trees made using SNP data to the trees made using indel data to ensure they are...
View ArticleHow To Call Variants (Snp, Indel, Sv) On A Bac Contig Aligned With Bwa Mem To...
I aligned a BAC contig (assembled from sanger sequences) to a reference genome using BWA-mem. The output alignments are very similar to the best end to end alignment I got from aligning the bac contig...
View Article