Question: Is there a tool to convert soft clipping to hard clipping in bam files?
gravatar for Ryan Thompson
3.8 years ago by
Ryan Thompson3.4k
TSRI, La Jolla, CA
Ryan Thompson3.4k wrote:

My lab is developing a pipeline to analyze some paired-end methyl-seq data. One of the steps is to soft-clip overlaps between read pairs to prevent downstream tools from double-counting the overlapped region, using bamutil clipOverlap ( However, one of the downstream tools, methylKit, will not accept soft-clipped reads, so I would like to convert the soft clipping to hard clipping in order to input it into methylKit. Is there a tool for converting soft clipping into hard clipping? Or is the conversion simple enough to implement myself?

bam clipping • 1.6k views
ADD COMMENTlink modified 3.8 years ago • written 3.8 years ago by Ryan Thompson3.4k

FYI, PileOMeth has a --methylKit option that will create methlyKit-compatible output and not require you to do any mucking around with artificially hard/soft clipping bases. Perhaps that already does what you need.

ADD REPLYlink written 3.8 years ago by Devon Ryan94k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1890 users visited in the last hour