If base quality scores set at 2 for clipped regions are restored to the original score with MergeBamAlignment, how will clipped regions not contribute to downstream analyses?
PlannedDear GATK Team,
In the tutorial (How to) Map and clean up short read sequence data efficiently, it states that MergeBamAlignment (Picard) 'restores original base quality scores'. My understanding is that this would include changing the base quality scores set at 2 for clipped regions (e.g. adapter sequences) marked with the XT tag back to original scores. Does this mean clipping status is lost and therefore, adapter sequences will be considered in downstream analyses?
If this is the case, BQSR, HaplotypeCaller, Mutect2 and other downstream tools could potentially consider problematic sequences (e.g. adapter sequences) in contrast to what has been described in previous posts. Therefore, clarification would be appreciated.
Thank you for your time and help.
Kind regards.
Please sign in to leave a comment.
0 comments