Hi everyone, here's my question:
I have some BAM files resulting from mapping RNA-seq reads. I would like to know whether there is a way to extract the reads mapping within defined intronic regions.
When simply doing "samtools view file.bam chr1:Start-End > result.bam", I get reads in the output file even if there are no actual reads mapped within the intronic region of interest, this is a result of the fact that samtools outputs also read fragments spanning from one exon to another, and not mapping inside the intron; while I would like to know which reads map within a given intronic region.
Does anyone know how I may do this?
Thank you in advance.
•
written
by
gspirito •
1 answer
I quickly wrote something: http://lindenb.github.io/jvarkit/Biostar9469733.html
this tool looks at the cigar string and check if there any block has an overlap with the bed file.
java -jar dist/biostar9469733.jar --regions input.bed input.bam
Log in to answer this question.
-
How to extract intron counts from total RNA Sequencing?
written by PK •
-
Reads mapping to pseudogenes using RSubread
written by cstashko •
-
Histone marks enrichment analysis
written by gspirito •
-
Extract reads that are contained entirely within a give region
written by george.wiggins •
-
Extract only RNA-seq reads in a BAM file that extend into an intron
written by foxdie •
-
Identify reads that span intron-exon junctions in RNA-Seq
written by oriolebaltimore
-
Re-merge split BAM files
written by KT •
-
featurecounts vs samtools view
written by tonja.r
-
Extract reads within given region, and their mates
written by user230613
-
Count splice junction reads from SAM/BAM files
written by juncheng
