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 there are reads with indels either at/overlapping the position, or in the vicinity of the position, say 5 bp upstream or downstream. This is a read level query, it's not indel calling. Basically I want to for a given position consider the ±5 bp window around the position and ask if there are any indels within this window, at the level of reads, and if so count the number of indels in the window.
I think the CIGAR strings of the BAM file probably holds all the information needed, but the question is how to do this in a simple, practical way.
•
written
by
Anders Bergström •
1 answer
You could use the program I previously wrote for biostar. This program produces a XML description of the Reads overlapping a given genomic region. See the CIGAR elements under the 'align' tag.
•
written
by
Pierre Lindenbaum
Log in to answer this question.
-
Bowtie Vs Bwa For Indels
written by Varun Gupta
-
Read File Format .Bam
written by hudarembang •
-
How To Sort Sam File
written by Chen Sun
-
Sam To Bam - Loss Of Data Or Just Great Compression?
written by Clare
-
Filtration Of Reads With Length Lower Than 30 From Bam
written by filipzembol
-
Solution To Find Out Read Orientation As Well As Strand To Which Reads Align From Sam File
written by k2bhide
-
Rna-Seq Split Gap Size Cutoff?
written by Alex Reynolds
-
Library To Parse Absolute Reference Positions Of Indels From A Bam Cigar String
written by William
-
Merging Bam Files
written by rickyflintoff •
-
How To Interpret The Flag Field In The Sam Specification
written by deepthithomaskannan
