Indel 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 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 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 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 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 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 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 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 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 ArticleHow To Transform Adjacent Complex Variants In Vcf Format To Individuals'...
Hello, everyone. Recently, I have downloaded the latest variant data of 1000genomes project from ftp://ftp.ncbi.nlm.nih.gov/1000genomes/ftp/release/20110521/ and the mapped data from...
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 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 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 ArticleGenomic Alignment And Snp/Indel Calling - My First Ever "Pipeline"
So I have finally generated some reads and run it through what I guess could be called a very rudimentary 'pipeline'. I generated a million paired end reads with wgsim then aligned with bwa, and used...
View ArticleHow To Call Variants (Snp, Indel, Sv) On A Bac Contig Aligned With Bwa Mem To...
I aligned a BAC contig (assembled from sanger sequences) to a reference genome using BWA-mem. The output alignments are very similar to the best end to end alignment I got from aligning the bac contig...
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 ArticleTool: Phylosim - Monte Carlo Simulation Of Sequence Evolution In The R...
PhyloSimPhyloSim is an extensible object-oriented framework for the Monte Carlo simulation of sequence evolution written in 100 percent R. It is built on the top of the R.oo and ape packages and uses...
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 ArticleTools For Large Indels Discovery In Single-End Rna-Seq Reads
Hi all,Has anyone heard about any tool wich could deal with large indels discovery? I've found some for paired-end reads but never for single-end... My final goal is to find transposable elements...
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 Article