INNER CODE UNIT · Rust

check_sequence_consistency

informationsea/transanno · liftover-rs/src/chain/mod.rs:125

    pub fn check_sequence_consistency<G1: GenomeSequence, G2: GenomeSequence>(
        &self,
        original_sequence: &mut G1,
        new_sequence: &mut G2,
    ) -> Result<(), LiftOverError> {
        // Check chromosome length
        if let Some(expected_len) = original_sequence
            .get_contig_list()
            .iter()
            .filter(|x| x.0 == self.original_chromosome.name)
            .map(|x| x.1)
            .next()
        {
            if expected_len != self.original_chromosome.length {
                error!(
                    "Length of {} in chain file is {}, but {} in original FASTA is {}",
                    self.original_chromosome.name,
                    self.original_chromosome.length,

View source record →