YeoLab / eCLIP

Other
38 stars 26 forks source link

Why only read2 is used for peak calling? #43

Open FionaMoon opened 10 months ago

FionaMoon commented 10 months ago

Hi developer,

As I ran eclip pipeline with pare-end data, I noticed that this pipeline only outputs the second read in each pair for use with peak caller.

According to eclip experiment, the crosslink site should be in read1, like this:

image

So, I wonder why we call peaks based on read2 but not read1?

Thank you in advance!