Strange variant call it show a indel but looks like a SNP
Hello all, i have a variant from a ION PGM sequencing, called with the ION Torrent variant caller plugin CHR Â Â POS Â Â Â ID Â Â REF Â Â Â ALT chr22 Â 44324726Â . Â Â Â TÂ Â Â Â Â Â TG The flow cell...
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 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 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 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 and missing in mach phased output
In mach phased output, how does it represent missing genotypes and indels?Â
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 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 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 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 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 ArticleVcf Snps And Indels Mouse Genomes
Hi, I am using SNPs and INDELs data from the mouse genomes project and the data seem to have conflicting information. I mean, seems they report a SNP in a deleted position. Thus my question, Is it...
View ArticleDeletions, Translocations And Cnv From Rna-Seq
Is it possible to find out deletions tarnslocations and / or copy number variations between RNA-seq data of control vs treated samples. The data is PE with enough coverage ( say around 30X).Thanks
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 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 ArticleSimulating A Read With A Deletion Of 50 Bp Into It.
Im trying to test some indel calling programs with large deletions in sequencing reads. I ve got this sequence in Chromosome 20 of hg19 : > CTAGCAAGGG GGCTGTATGG CTTGAGGCCA TAGTCCAGGA CATCATCGGG...
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 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 ArticleComparative Genomic Analysis Of Bacterial Field Strains
Hi, I have a question on comparative genomics of bacterial field strains.If there is a very low number (only 30) of variants (base substitutions, indels etc) in a comparative genomics study of several...
View ArticleGatk Threshold Parameter For Clipping Of Reads During Indel Calling Running...
GATK 1.6 claims to apply a Q20 threshold to clip ends of the reads for the indel caller:http://www.broadinstitute.org/gatk/guide/article?id=12373. Indel Calling with the Unified Genotyper [...] while...
View Article