Somaticindeldetector Vs Unifiedgenotyper
Hey Guys, I've been establishing my own first pipeline for calling variants: Indels and SNPs. For that, I used mainly GATK tools. I'm doing this for whole genome and whole exome sequences of mouse...
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 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 ArticleHow 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 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 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 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 ArticleClassify Mutations Which Are Found Ny Using Exome Sequencing
Hi, everyone. I used exome sequencing to sequence tumor and blood, and after I remove mutations in blood and snp in dbsnp131, I found I still have about 2000 mutations. I read the exome sequencing...
View ArticleHow To Make A Circular Plot Of Snps And Indels ??
Hello everybody,I am analysing SNPs and INDELs in several prokaryotic genomes of the same strain. My goal is to compare them with a reference genome . I would like to learn how to plot the SNPs in a...
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 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 ArticleHow Are Indel Coded In Vcf Files
Hello,I have a question concerning how are INDEL coded in VCF files. I have for instance the following file:CHROM POS ID REF ALT QUAL FILTER8 18078835 . TTA T 46 PASS8 18078836 . TA T 138 PASSHow to...
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 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 ArticleIndel Notation In Variant Calling
HelloI am sorry for the basic question but I am struggling to find any details of the nomenclature of indels by variant calling software. Unfortunately i am unable to access the details of the software...
View ArticleGatk Somaticindeldetector Vcf Info Field Missing
Hey everyone,I am a bit confused. I called my INDELs pairwise using the SomaticIndelDetector. But in the resulting VCF the info field seems to be missing, is this normal? Or am I getting something...
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 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 ArticleLocating Indels In Gene
I have a excel file from windows (spreadsheet) contains start, end and indel sequences (-GGC means GGC deletion or +AAC means AAC insertion) in column 1, 2 and 3 respectively. I have a reference genome...
View Article