1000 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 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 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 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 ArticleIs It Possible To Get Rs Numbers For Indels?
I am currently designing an Illumina custom GoldenGate platform (for the first time, so with lack of experience..). I have finished the selection of SNPs and I'm ready to submit the list for...
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 ArticleSamtools: 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 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 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 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 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 Difference Between Gatk And Dindel For Calling Indel
What is difference between GATK and Dindel for calling indel ,I see that in the GATK second step :discovery indel ,The genotype likelihoods calculation is inspired by Dindel.So I want to know the...
View ArticleAnyone Know A *General* Variant (Snp/Indel) Effect Predictor / Annotator?
I'm looking for a general variant annotator ... one that would call a variant as synonymous, non-synonymous, non-genic, UTR, frame-shift, truncation, etc. By general I mean that we should be able to...
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 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 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 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 ArticleCalculating Tn (True Negatives) For An Indel Detection Method
This is a kind of follow-up inspired by the very good question/answers I read in "How to calculate sensitivity/selectivity of an algorithm that returns locations of possible matches?"My goal is to...
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 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 Article
