trim_hsp65_primers

Function trim_hsp65_primers 

Source
pub fn trim_hsp65_primers(seq: &[u8]) -> Vec<u8> 
Expand description

Trim hsp65 amplification primers from a sequencing read.

Unlike trim_start_end, which takes a single primer per boundary, hsp65 amplification uses two alternative primer pairs (Telenti and Patel sets), so each boundary is defined by a set of candidates — the earliest match wins for the start and the latest for the end.