Programming 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 ArticleHow To Identify/Flag The Presence Of Indels On The Read Level Nearby A Given...
Let's say I have list of genomic positions, e.g. SNP calls in a VCF file, and a BAM file with read alignments. What would be the best way to, for each of those positions, determine if in the BAM file...
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 ArticlePindel Error: Shiftedvector Out-Of-Range Error
Hello,I’ve been running pindel like so: pindel --config-file bam_config.txt --fasta /path/to/hg19/fasta/ --chromosome ALL --output-prefix /path/to/outputand I’ve noticed a lot of my pindel jobs fail...
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 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 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 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 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 ArticleHow To Transform Adjacent Complex Variants In Vcf Format To Individuals'...
Hello, everyone. Recently, I have downloaded the latest variant data of 1000genomes project from ftp://ftp.ncbi.nlm.nih.gov/1000genomes/ftp/release/20110521/ and the mapped data from...
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 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 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 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 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 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 ArticleBest Way To Compare Output Snp/Indels From Different Software?
I have a list of SNPs and indels from 4 different software. IndelGenotyper (broad tool)—for calling somatic indels Bambino—for calling SNP/Indels Somaticsniper—SNP calling Varscan – indel/snp/LOH Now...
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 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 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 Article