Hello, is there any way to remove low quality score bases from bam files, leaving the space blank or replaced with a filler character (i.e. Ns)?
I have mapped my sequencing reads, and adjusted the base quality scores (using mapDamage, which is conceptually similar to GATK's BQSR). I would now like to remove bases that have quality score less than 20 prior to visualizing the misincorporation rate, regardless of what the read's mapping score was (understanding that there can be poorly called bases on a confidently mapped reads). Is there a tool to do this? Thanks in advance!
Awesome, works like a charm. Thank you!