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 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 ArticleDifferences Running Dindel On Different Unix Environments
Hi all, I am hoping to try out dindel for indel-calling on my NGS data but coming up against some problems getting it to work on both institute's computing cluster and my desktop machine. On the...
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 ArticleGatk Threshold Parameter For Clipping Of Reads During Indel Calling Running...
GATK 1.6 claims to apply a Q20 threshold to clip ends of the reads for the indel caller:http://www.broadinstitute.org/gatk/guide/article?id=12373. Indel Calling with the Unified Genotyper [...] while...
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 ArticleSimulating A Read With A Deletion Of 50 Bp Into It.
Im trying to test some indel calling programs with large deletions in sequencing reads. I ve got this sequence in Chromosome 20 of hg19 : > CTAGCAAGGG GGCTGTATGG CTTGAGGCCA TAGTCCAGGA CATCATCGGG...
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 ArticleChip-Seq Peak Caller Able To Find Peaks Split By Indels?
Are there any chip-seq peak callers able to find cases where the region under a peak has been split by an indel between the reference and the chipped sample?For example, let's suppose one has done...
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 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 Is Difference Between Gatk And Dindel For Calling Indel
What is difference between GATK and Dindel for calling indel ,I see that in the GATK second step :discovery indel ,The genotype likelihoods calculation is inspired by Dindel.So I want to know the...
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 ArticleDeletions, Translocations And Cnv From Rna-Seq
Is it possible to find out deletions tarnslocations and / or copy number variations between RNA-seq data of control vs treated samples. The data is PE with enough coverage ( say around 30X).Thanks
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 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 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 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 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 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 Article