Filter BAM to keep only alignments at an interval, ignoring reads spliced over
Hi,
I am looking at a rare splicing event, and would like to keep only reads in a specific interval that have an alignment at that location.
With a simple samtools view alignment.cram chr1:12345-12456 I also get all reads that are spliced across that interval (screenshot below).
How do I get rid of reads that do not actually align at this location?
Bonus points if I can still keep the mates of the reads that align at this location...
• 1,272 views
•
link
written
by
WouterDeCoster
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Help manually processing strand in paired-end reads when fishing for lariats in bulk RNA-seq
written by txema.heredia •
-
Tutorial: Yes .. BBMap can do that! - Part I : bbmap (aligner), bbduk (scan/trim), repair (fix PE r…
written by GenoMax
-
Event- vs transcript-based approaches for analyzing alternative splicing
written by fmerkal •
-
Unexpected hump in GC content distribution from FasQC report in ChIP-seq
written by msimmer92
-
De novo RNA-seq
written by blawney •
-
News: Metabarcoding: Transforming the way in which biodiversity can be surveyed
written by Physalia-courses
-
Tutorial: InsideDNA: RNA-Seq de novo assembly using Trinity
written by shilparaopradeep •
-
Unable To Replicate Splice Junction In Tophat
written by Dan D
-
Tutorial: Mapping Reads With Bwa And Bowtie
written by MSU course 2011 •
-
Intron Retention And Alternative 5'3' Splice Site Identification
written by Arun
what's the gray section ? it's the read you want ? do you only want the reads that end at the splice event (single point location) ?
I want reads in a specific short interval, but if I do that with samtools view I also get reads that are spliced over that interval.
see also : Extract paired-end reads with one end falls within a given region from a bam file ; Remove reads from bam file that partially covers a region ; Is it possible to remove reads from BAM which match the reference at a position