Deletions, 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 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 ArticleFrequency 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 ArticlePython - finding INDELS and SNPS
seq_a = 'GAGAGATTTTCCAATTCGACG-------CGGGGTCAGG--GAAATTT'Â Â seq_b = 'GAGAGATTGGCCTTAACTACCCAACCCACGGCCTGACCGAGGTCTTC'G A C T = Bases- = INDELPYTHON - I am very new to programming and need some help, I...
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 Indels -- Filtering Only Hits With Insertion In The Reference In The...
Can anyone suggest how to use samtools to filter only hits where there is a insertion in the reference that splits the sequence hit roughly by the middle? My sequences are in the range of 100-1000bp...
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 ArticleSamtools / Bcftools Missing Random Obvious Single Base Indels
I am using samtools to call variants in a haploid genome (yes I know it is designed for diploid). It finds SNPs easily, and most of the indels, but there are a few indels it can NOT seem to find no...
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 ArticleNext Gen Sequencing Indel Realignment With The Genome Analysis Toolkit (Gatk)
As I continue to add steps to my SNP/Indel discovery workflow, the latest recommendation is local realignment around Indels using GATK following the initial alignment step. I have just commenced the...
View ArticleIndel Detection For 454 Resequencing
I do some work for a small diagnostics company that has a requirement for small indel detection in 454 data. Those of you familiar with Roche's pipeline will be aware that AVA (Amplicon Variant...
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 ArticleInsertion Detection Tools
Hi,Please suggest me some tools which can detect large indels (esp. of size > 100bp). Some of the tools I have already tried include Breakdancer and Pindel and Hydra. Any other suggestions??
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 ArticleWhat Do You Do With Indels Data?
I have a list of insertion and deletion from a whole genome sequence obtained by next generation sequencer, chr, pos and sequence for each indels aligned to our reference genome.I'm not familiar with...
View ArticleLarge Deletions With Bwa And Single End Reads
Hi, I am trying to detect deletions approximately 66bp in length by mapping 90 base reads to the human genome with BWA (and using GATK). By fine tuning the alignment parameters I have managed to...
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 ArticleGenerate pileup for indels only using SAMtools mpileup
Is there a way to generate a pileup for indels only using the samtools mpileup command? The -I flag can be used to skip indel calling but there is no option to call indels only.Thanks.
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 ArticleLocating Indels In Gene
I have a excel file from windows (spreadsheet) contains start, end and indel sequences (-GGC means GGC deletion or +AAC means AAC insertion) in column 1, 2 and 3 respectively. I have a reference genome...
View Article