What Do You Expect As The False Positive And Negative Rate For Snp'S And...
What do you expect as the false positive and negative rate for SNP's and INDELS in a WGS experiment?On which papers and data sets (inhouse or external) do you base this? Edit: Of course this depends on...
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 ArticleComparing Genomic Assemblies To Find Indels
Hello everyone,We've sequenced a YAC containing a stretch of human DNA as an insert. We know what the insert sequence is supposed to be. However - as you could imagine - Illumina-based sequencing has...
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 ArticleGatk3 Vs Dindel For Small Indel Calling
I read some posts from Heng Li around two years ago that he suggested Dindel is the optimal software for small indel calling. I'm just wondering if this is still true today because GATK seems to...
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 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 ArticleSNP and In/Del Counting without Raw Sequence Data
Is there an easy way to do SNP calling on multiple sequences of a gene (from different members in a population) without the raw sequence data? Ā Most of the tutorials online involve starting from the...
View ArticleSequence Variation Detection For Prokaryotic Genomes
I have resequenced a bacterial genome using SOLiD paired end reads. After looking at SNPs and small Indels, I'm now interested in large Insertions/Deletions (Size of genes and gene clusters).Which...
View ArticleTrio: Detecting A Large De Novo Indel
I've been given a set of three BAM (father, mother, child) and I expect the child to contains some de-novo heterozygous variations . samtools mpileup have been used to find the small variations but...
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 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 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 ArticleDo You Guys Have Any Recommendation On Variant Calling Using Long Contigs Or...
I have assembled a human genome and used bwa mem aligning the contigs back to the reference genome. Do you guys have any recommendation of variant caller?
View ArticleHow To Detect Overlaping Indels In Vcf File?
Hi all,I have a vcf file with indels from a whole genome analysis and I want to detect overlapping between indels... I tried to use BEDtools-intersect but it asks for 2 files, and I only have a single...
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 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 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 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 ArticleFiltering Vcf File
I was wondering how to filter the vcf file based on a few input arguments ( DP>10, MQ>30 and QD>20 or GT = "1/1" etc)? I m planning to use simple command on the command line to extract the...
View Article