Question

Filter Exon Aligned Reads That Are Mapped Against Human Genome.

0

Entering edit mode

10.7 years ago

Coryza ▴ 30

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.3k views

ADD COMMENT • link updated 3.6 years ago by Ram 44k • written 10.7 years ago by Coryza ▴ 30

Ram · Answer 1 · 2014-03-27

0

Entering edit mode

10.7 years ago

Pierre Lindenbaum 164k

duplicate of Samtools and Region List

ADD COMMENT • link 10.7 years ago by Pierre Lindenbaum 164k

0

Entering edit mode

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

ADD REPLY • link updated 5.1 years ago by Ram 44k • written 10.7 years ago by Coryza ▴ 30

0

Entering edit mode

The online manual is two versions behind.

ADD REPLY • link 10.7 years ago by Devon Ryan 104k

0

Entering edit mode

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

ADD REPLY • link 10.7 years ago by Coryza ▴ 30