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 (http://genome.sph.umich.edu/wiki/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?
Question: Is there a tool to convert soft clipping to hard clipping in bam files?
18 months ago by
Ryan Thompson • 3.3k
TSRI, La Jolla, CA
Ryan Thompson • 3.3k wrote:
ADD COMMENT • link •
Please log in to add an answer.
Powered by Biostar version 2.3.0
Traffic: 1339 users visited in the last hour