Samtools: Local Realignment Around Indel
In samtools mpileup, adjust the option -m and -F to control when want to initiate indel realignment, any guideline to set the values for m and F for high coverage data (eg. >=20X)low coverage data...
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 ArticleWhat Do You Expect As The False Positive And Negative Rate For Snp'S And...
What do you expect as the false positive and negative rate for SNP's and INDELS in a WGS experiment?On which papers and data sets (inhouse or external) do you base this? Edit: Of course this depends on...
View ArticleDeletions, Translocations And Cnv From Rna-Seq
Is it possible to find out deletions tarnslocations and / or copy number variations between RNA-seq data of control vs treated samples. The data is PE with enough coverage ( say around 30X).Thanks
View Article1000 Genomes Indel/Snps File
Is there any place to download 1000 genome SNPs and Indels data in a single merged file for each pilot (pilot 1 , 2 , and 3) ? I am trying to use GATK and thought to use 1KG indels and SNPs vcf files...
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 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 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 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 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 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 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 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 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 ArticleReducing False Positive Indel Calls From Pindel
I am analyzing the specificity of Pindel's indel calls. I have matched tumor/normal WGS data, bwa aligned. The original calls were filtered to remove any that had supporting reads from the normal...
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 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 ArticleProgramming Challange: Pairwise Alignments To Multiple Alignment
I have a set of 10-12 very closely related chromosome sequences (from different strains) aligned to a "single" reference chromosome. Now I need to generate multiple sequence alignment of these without...
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