BWA 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 ArticleUnknown Indel Notation To Hgvs
I'm sorry about this caotic question, but I don't understand nothing at all... I have a file of Indels finded in genes, presumably from HGMD. Well, i don't know how to read this file, i don't have any...
View ArticleHow To Handle Reads Ending With Deletions In Gatk?
Hello,before asking my question, I should point out that I'm working with data that's not my own (publicly available), to learn and establish a proper workflow when real data wlll arrive in the...
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 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 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 ArticleNovel Sequence Insertions 1000 Genomes Project
I am looking for novel sequence insertions identified in the 1000 genomes project, and I found 3 files in this...
View ArticleWhat Kind Of Analysis Can Do You With Vcf Files?
What kinds of analyses can be completed if you have SNP, indel, and SV data in vcf files? I just need a few very broad and generic examples.
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 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 ArticleWhat Do You Do With Indels Data?
I have a list of insertion and deletion from a whole genome sequence obtained by next generation sequencer, chr, pos and sequence for each indels aligned to our reference genome.I'm not familiar with...
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 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 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 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 ArticleFunctional Consequence of In-frame Indel
Does anybody include in-frame indels in association studies? A quick pubmed / google / biostars search on this didn't come up with much. In-frame indels are often treated as synonymous, but since they...
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 ArticleMulti-Sample Indel Realignment Using Gatk
Is there a way to use the GATK Indel Realigner tool with a two-samples-input? The issue is that in my normal sample, several reads are mapped containing SNPs at specific positions, and ONE read having...
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 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 Article