Is there a tool to convert soft clipping to hard clipping in bam files?
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?
• 3,414 views
•
link
written
by
Ryan Thompson
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Bamutils/ClipOverlap
written by chrisk •
-
How to properly mock a (Pysam) read
written by ManuelDB
-
How to mark as QC fail reads with specific CIGARs
written by ManuelDB
-
Samtools ampliconclip hard clip not giving same results as soft clip
written by Kevin •
-
converting Bam to fastq while removing clipping(hard/soft clip bases)
written by eli_bayat •
-
Tool: Trimviz, a lightweight QC tool to assess impact of read-trimming and aligner soft-clipping
written by stuart archer •
-
which annotation to use and is it advisable to use bam to fastq
written by novicebioinforesearcher •
-
Influence of RNA-seq read splitting on the BAQ calculation in samtools mpileup
written by Manuel Landesfeind
-
Fusion with strand orientation
written by mangfu100
FYI, PileOMeth has a
--methylKitoption 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.