How 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 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 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 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 ArticleWhat Is The Best Tool For Mouse (Mm9 Or Mm10) Variant Annotations?
I am working with some mouse exome capture data, and I would like to provide functional annotations to variants (SNPs and small indels).For human data, I would typically use ANNOVAR. It looks like...
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 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 ArticleCalculating Tn (True Negatives) For An Indel Detection Method
This is a kind of follow-up inspired by the very good question/answers I read in "How to calculate sensitivity/selectivity of an algorithm that returns locations of possible matches?"My goal is to...
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 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 ArticleStrelka Indel Allele Counts
Hopefully there are some Strelka users out there who can help with this one. I'm looking for allelic counts of the ref and non-ref allele of the indels.I see these flags in the VCF header, do these...
View ArticleHow To Set Filter For Frequency Of Reads And Hapmap Exome Sample Results:
Hey All, I only used so far three filters for my whole exome pipeline (aligning to hg19) for a HapMap sample. I tried it on the NA19240 Hapmap sample from paper below (Table 3) which shows ~196...
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 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 ArticleBest Way To Compare Output Snp/Indels From Different Software?
I have a list of SNPs and indels from 4 different software. IndelGenotyper (broad tool)—for calling somatic indels Bambino—for calling SNP/Indels Somaticsniper—SNP calling Varscan – indel/snp/LOH Now...
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 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 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 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 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 Article