Counting Number Of Bam Reads Directly Within Set Of Intervals With Bedtools
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 not want reads overlapping the intervals, but ones that fall directly within them.
I tried the following:
intersectBed -abam reads.bam -b exons.gff -wb -f 1
this has redundancies, so I pipe it into coverageBed as follows:
intersectBed -abam reads.bam -b exons.gff -wb -f 1 | coverageBed -abam stdin -b exons.gff
Is this correct? Thanks.
• 6,529 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
User 9996
•
written
by
User 9996
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Question about number of reads within intervals
written by iraun
-
Bowtie2 Mapping Different Number Of Reads To Same Sequence When Ref-Seq Is Part Of Different Indexes
written by rohita.sinha •
-
Getting All Reads That Align To A Region In Compact Bed Format Using Bedtools?
written by user
-
Getting Number Of Reads In Intervals With Bedtools
written by user
-
Intersectbed/Coveragebed -Split Purify Exon?
written by Puriney
-
Extract Only Paired-End Reads That Map A Specific Interval
written by Nicolas Rosewick
-
Filtering Stream Bam Output Idealy From Pysam/Python
written by User 9996
-
Intersectbed Overlap
written by Nicolas Rosewick
-
Reporting The Bam Reads Overlapping A Set Of Intervals With Bedtools
written by User 9996
-
Bedtools Compare Multiple Bed Files?
written by Bioscientist
Yes, that is correct. I suggest using version 2.13.1 and use the -counts option in coverageBed. This will use less memory and run much more quickly than the default coverageBed.
is there a caveat about 0-based versus 1-based coordinates of SAM and BED format?
Both BAM (not SAM, which is 1-based) and BED are zero-based, half open. BEDTools handles this automatically anyway.
How can this be extended to get all mapped pairs that land in the boundaries of a certain intervals from BAM?