Help with Piranha bin size for BED files
0
0
Entering edit mode
7.8 years ago
anu014 ▴ 190

I am using Piranha for Clip-seq analysis. I'm not having score column for my BED file. How much my 'bin-size/-b' should be? Can I take difference between 'start' & 'end' coordinates?

The Piranha manual says that : The response file (always the first argument, and required) may contain: Raw read locations, in which case Piranha will bin the reads for you. If you provide your input like this, you MUST set the bin size option, or Piranha will treat your input as being already binned. If Piranha is creating bins from raw reads, it will always start at the first index of each chromosome (i.e. index 1) and move at a step size equal to the bin size. Bins with no reads in them will not be retained. link : https://github.com/smithlabcode/piranha

Command used : ./piranha-1.2.1/bin/Piranha -b 43 -s ./G/SRR1690969.bed > piranha_peaks_output

My BED file looks like this :

chr1 630917 630960 . . + chr1 632324 632367 . . + chr1 632327 632370 . . + chr1 632328 632371 . . + chr1 632329 632372 . . + chr1 632331 632374 . . + chr1 632332 632375 . . + chr1 632333 632376 . . + chr1 632334 632377 . . + chr1 632336 632379 . . + chr1 632340 632383 . . + chr1 632342 632385 . . + chr1 632343 632386 . . + chr1 632344 632387 . . + chr1 632346 632389 . . +

Thank you in advance.

sequencing alignment ChIP-Seq next-gen • 2.5k views
ADD COMMENT
0
Entering edit mode

I have the same question. I get "ERROR: binning reads failed. Reason: size of required bins doesn't match size of actual bins" for the following command: ./Piranha -s -i 100 SMUG1.bed \ IgG.bed \ -o peaks.bed where 100 is bin size = read length.

ADD REPLY

Login before adding your answer.

Traffic: 1983 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