Entering edit mode
9.6 years ago
rmccloskey
▴
240
I have a number of BAM files with a very high number of reads flagged as duplicates (>50%). To confirm that these are indeed duplicates and the flags were not introduced erroneously during processing, I would like to clear all the duplicate read flags and run Picard tools MarkDuplicates to recheck these. Is there a better way to change the flags in a BAM file then to convert to SAM, edit the file with sed or something, and then convert back to BAM?