Function reverse_complement 

pub fn reverse_complement(seq: &str) -> String

Returns the reverse complement of a DNA sequence.

This helper preserves non-ACGT characters as-is because the reference and read sequences may contain N, and parity code expects those positions to survive.