Benchmarking Read Alignment And Variant-Calling Algorithms (For Dummies)
Hi all, I am wondering if there is a good step by step guide of how to benchmark alignment and variant calling software. I do understand the premise e.g.Generate reads with known mutations Align to...
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 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 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 ArticleSomaticindeldetector 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 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 ArticleIs It Possible To Get Rs Numbers For Indels?
I am currently designing an Illumina custom GoldenGate platform (for the first time, so with lack of experience..). I have finished the selection of SNPs and I'm ready to submit the list for...
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 ArticleMerging/Concatenating Vcf Files
I have a vcf file of SNPs and another vcf file for my Indels. During SNP calling step by GATK, I created them separately instead of using -glm (calling together SNPs and Indels). Now to get the...
View ArticleHow To Identify/Flag The Presence Of Indels On The Read Level Nearby A Given...
Let's say I have list of genomic positions, e.g. SNP calls in a VCF file, and a BAM file with read alignments. What would be the best way to, for each of those positions, determine if in the BAM file...
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 ArticleMerging Vcf Files And Filling In Missing Genotypes
I have vcf files with genomic variants (e.g. SNPs and Indels) from several patients. I now want to merge them into one vcf file e.g. using vcftools' merge-vcf. Doing that I get a file that holds for...
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 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 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 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 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 Article1000 Genomes Indel/Snps File
Is there any place to download 1000 genome SNPs and Indels data in a single merged file for each pilot (pilot 1 , 2 , and 3) ? I am trying to use GATK and thought to use 1KG indels and SNPs vcf files...
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 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 Article