Filter Exon Aligned Reads That Are Mapped Against Human Genome. · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Filter Exon Aligned Reads That Are Mapped Against Human Genome.

Hi,

I am wondering if it is possible via samtools/bamtools to extract all the paired-end mapped reads (against the human hg19 genome) that fall in exon regions. I have a BED file that defines all these regions. I already know that it's possible to calculate this (coverageBed -abam reads.bam -b exons.bed > exons.bed.coverage), but I have no clue how to extract the reads from the .bam file). I tried using the samtools/bamtools manual but couldn't figure it out.

Anyone who can help me?

Thanks

exons genome mapping hg19 bowtie2

• 3,849 views • link

updated 5.3 years ago by Ram

• written 12.4 years ago by Coryza •

1 answer

duplicate of Samtools and Region List

• 0 views • link

written 12.4 years ago by Pierre Lindenbaum

Hmm wonder why I could not find that in the samtools manual..

• 0 views • link

updated 6.8 years ago by Ram

• written 12.4 years ago by Coryza •

The online manual is two versions behind.

• 0 views • link

written 12.4 years ago by Devon Ryan

Ah well, thanks for the help both of you ;)

• 0 views • link

written 12.4 years ago by Coryza •

Log in to answer this question.

More posts like this

exon.bed file for coverage
written 5.3 years ago by Sabeen •

Hello, I am trying to find coverage for my exon.bed file. for this i have used bedtools by following code **bedtools bamtobed -i reads.bam | …
How to properly get BED regions coverage
written 6.7 years ago by MAPK

I have a sorted bam file and a bed file with chromosome regions. I would like to: > 1. get mean coverage for the regions …
Filtering paired and single end reads from a BAM file
written 11.3 years ago by James Ashmore

<p>I have a paired-end BAM file which has gone through numerous filtering steps (MAPQ, MarkDuplicates, Mitochondrial, Unmapped). I would like to extract all properly paired …
Question about number of reads within intervals
written 11.8 years ago by iraun

<p>Hi there,</p> <p>This question is very basic but I need to ensure that I'm going on the right way. I need to calculate the number …
Remove Intronic Regions in .BAM
written 12.2 years ago by K.Nijbroek •

Hi I have a .BAM file which contains discordantly and concordantly mapped mate-pairs. I used bedtools Pairtobed to extract the mate-pairs which both show overlap …
Calculating Exome Coverage
written 12.4 years ago by Coryza •

<p>*// Edit to make the post more clear (Mapping done via Bowtie2). My problem is that when counting Exome Coverage via coverageBed gives different results …
Getting All Reads That Align To A Region In Compact Bed Format Using Bedtools?
written 13.6 years ago by user

<p>I'm trying to find all the reads (by name) from a BAM file that align to various regions in a bed file. Right now I …
Intersectbed/Coveragebed -Split Purify Exon?
written 13.9 years ago by Puriney

<p>all.reads.bam file records mapped RNA-seq reads data, including: </p> <ol> <li>exon:exon junction</li> <li>exon body</li> <li>intron body</li> <li>exon:intron junction</li> </ol> <p><strong>Q1:</strong> When calculating RPKM for given …
Reporting The Bam Reads Overlapping A Set Of Intervals With Bedtools
written 14.8 years ago by User 9996

<p>I am trying to use <a href='https://code.google.com/p/bedtools/'>bedtools</a> to pull out the reads falling directly within a set of BED coordinates. While this command does it …
Counting Number Of Bam Reads Directly Within Set Of Intervals With Bedtools
written 14.9 years ago by User 9996

<p>how can I count the number of BAM reads falling directly within a set of intervals, given in a GFF format? Note that I do …