Frequency 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 ArticleSnps And Indels In Homozygous Genomes
I have a genomic sequence of an inbred mouse strain i.e. the genome is homozygous for a given position or same alleles are present on both of the strands. Now for variant calling, I have to compare it...
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 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 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 ArticleHow to call low frequency snp/indel with GATK's HC/UG
i want to call snp/indel with low frequency, like DP=100, ALT=10.how to set HC/UG ?with UG, i run with :--min_indel_count_for_genotyping 2--min_indel_fraction_per_sample 0.05get a nice result, but i...
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 ArticleBwa Aln For Snp/Indel Discovery With Paired End 100Bp Human Dna Reads
Hi, Can anyone show me a recommended command line input to achieve this? I have never used the program before and am afraid of entering inappropriate values for one of the many paramters.Or is there a...
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 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 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 ArticleHow To Retrieve Reads Supporting A Snp/Indel
I used samtools and varscan for variant calling like this:samtools mpileup \ -f genome.fa \ tku.sorted.rmdup.bam \ | java -jar VarScan.v2.3.5.jar pileup2snp \ --min-coverage 8 \ --min-reads2 2 \...
View ArticleSplit Vcf File Into Snps And Indels
Hi there,As recommended in the GATK best practices the Variant Quality Score Recalibration has to be done separately for SNPs and Indels. But, I didn't find the way to do this split in a clean way (for...
View ArticleAlternate Sources Of Reference Indel Vcfs Outside Of Dbsnp And 1000 Genomes...
Hi,I am in a midst of testing some software tools for detecting indel from single sample genome data, and was wondering whether there is any source of reference indel files outside of dbSNP and 1000...
View ArticleWhat Is The Difference Between Samtools Mpileup And Pileup
As of samtools 0.1.8 'mpileup' appears to have superseded 'pileup' for variant detection. The samtools website now has pileup as deprecated. So i am wondering whether i should repeat my recent...
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 ArticleClik here to view.
Indel Left/Right Alignment
When matching indels between different VCF files (generated by different callers), there is this issue with left / right indel alignment, For example: Here is a real example for one indel (from the...
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 Article