Function simulate_ligation 

pub fn simulate_ligation(
    inputs: &[(&SequenceRecord, &[String])],
    selected: &[FragmentSelection],
    topology: Topology,
) -> Result<LigationReport, LigationError>