How Many Is Too Many? Germline And Somatic Coding Indels And Snps In Cancer...
Would anyone care to share their experience with variant calling in cancer genomics using tumor - normal pair to find somatic vs germline variants especially indels?I have been getting an unbelievably...
View ArticleInsertion Detection Tools
Hi,Please suggest me some tools which can detect large indels (esp. of size > 100bp). Some of the tools I have already tried include Breakdancer and Pindel and Hydra. Any other suggestions??
View ArticleChip-Seq Peak Caller Able To Find Peaks Split By Indels?
Are there any chip-seq peak callers able to find cases where the region under a peak has been split by an indel between the reference and the chipped sample?For example, let's suppose one has done...
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 ArticleBWA Small Indels vs Roche Newbler Assembler
I've been experimenting with using BWA to assemble 454 (1kb reads), comparing previous assemblies performed with the proprietary Roche assembler and I've noticed that BWA seems to be able to handle...
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 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 ArticlePython - finding INDELS and SNPS
seq_a = 'GAGAGATTTTCCAATTCGACG-------CGGGGTCAGG--GAAATTT' seq_b = 'GAGAGATTGGCCTTAACTACCCAACCCACGGCCTGACCGAGGTCTTC'G A C T = Bases- = INDELPYTHON - I am very new to programming and need some help, I...
View ArticleBenchmarking Read Alignment And Variant-Calling Algorithms (For Dummies)
Hi all, I am wondering if there is a good step by step guide of how to benchmark alignment and variant calling software. I do understand the premise e.g.Generate reads with known mutations Align to...
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 ArticleGatk Indel Realignment Taking Forever! Help
Hi all, I am working on exome capture data for barley (1.3Gbp). I am interested in variant calling to find out SNPs in my sample. I have used SAMTools SNP calling and things get done in ~1 hr whereas...
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 ArticleIs There A Convention For Representing Indels In Diploid Genome Sequences?
If an individual has a heterozygous indel how should that be represented in a sequence?I am looking for some examples along the lines ofA[C]GT or A[C/-]GT
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 ArticleNext Gen Sequencing Indel Realignment With The Genome Analysis Toolkit (Gatk)
As I continue to add steps to my SNP/Indel discovery workflow, the latest recommendation is local realignment around Indels using GATK following the initial alignment step. I have just commenced the...
View ArticleHow To Calculate Genotype Concordance Between Indel Call Sets?
How do people calculate genotype concordance between INDEL call sets? I have a NGS and BAC based INDEL call set (in vcf) but I get a higher false positive and false negative rate than mentioned in some...
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 ArticleDoes Anyone Know Of A Publication Reviewing Snp/Indel Finding Algorithms For...
Hi all, Is anyone aware of a good publication that reviews current methods for detecting SNPs/Indels from next gen sequencing data? I have done some searching but so far in vain. It would be great to...
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 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 Article