extract the number of reads around TSS
2
0
Entering edit mode
5.0 years ago
Ben ▴ 50

I have coded oen python script to extract the amount of reads around TSS from ChIP-seq sam file, but the speed is too slow. Does anyone have any suggestions to solve this question? Thanks in advance!

ChIP-Seq • 1.7k views
ADD COMMENT
3
Entering edit mode
5.0 years ago

Via BEDOPS:

$ bedmap --echo --count TSSs.bed <(sam2bed < reads.sam) > answer.bed

If you want to pad the TSSs, say by 1k bases on both sides, add the --range N option:

$ bedmap --echo --count --range 1000 TSSs.bed <(sam2bed < reads.sam) > answer.bed

This will likely increase the number of overlaps between reads and (padded) TSSs.

ADD COMMENT
1
Entering edit mode
5.0 years ago
h.mon 34k

bam-readcount.

ADD COMMENT

Login before adding your answer.

Traffic: 1353 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6