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 ArticleArlequin Input Frequency Data
Hi!I have frequency data for 50 biallelic loci and 500 samples in 3 populations and i can´t understand how to create an input file for arlequin with these data. If someone could help me it will be...
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 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 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 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 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 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 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 ArticleWhat Approach Would You Recommend For Large Indel Detection With Solid Data
I've been spending quite some time on following problem: I sequenced a bacterial genome using paired-end reads (SOLiD) and I have a quite good reference sequence. My goal is to detect changes in the...
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 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 ArticleIndel and missing in mach phased output
In mach phased output, how does it represent missing genotypes and indels?
View ArticleTool For Finding Indel In Which Part Of Genes
Are there any tools to help us determine whether an indel lies within a specific part of the gene ( Exon, Intron, ...)?
View ArticleNgs Tool Analysis For Snp&Indel Calling Of Pooling Data
What are the best and most used analysis tool actually, for SNP and indel calling of pooled data?I'm using: GATK Freebayes that have specified options for data in pool.Now I have found CRISP.Are there...
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 ArticleLarge Deletions With Bwa And Single End Reads
Hi, I am trying to detect deletions approximately 66bp in length by mapping 90 base reads to the human genome with BWA (and using GATK). By fine tuning the alignment parameters I have managed to...
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 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 Article