Check a bam file if it was hard clipped?
0
0
Entering edit mode
3.0 years ago

Is there a programmatic way to verify that a BAM file was hard clipped? The BAM file header does not show the -Y was used. Does this require manual verification in IGV?

Hardclipping BAM IGV • 979 views
ADD COMMENT
2
Entering edit mode

Does this require manual verification in IGV?

I guess you would not know that unless you compare to original reads.

How about checking CIGAR strings for H e.g. 6H37M?

ADD REPLY
0
Entering edit mode

The BAM file header does not show the -Y was used

what is that -Y please ?

ADD REPLY
1
Entering edit mode

From bwa site (https://www.mankier.com/1/bwa)

"-Y Use soft clipping CIGAR operation for supplementary alignments. By default, BWA-MEM uses soft clipping for the primary alignment and hard clipping for supplementary alignments."

ADD REPLY

Login before adding your answer.

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