sigalign

Similarity Guided Alignment Algorithm

Caveat

This library is currently under development. Some key features are missing, and functions can be changed anytime without notification.

What is Sigalign?

• Sigalign is pairwise alignment algorithm using gap-affine penalty for nucleotide and amino-acid sequences.

  Features

• Sigalign has two values that evaluate the similarity between sequences.

  • (1) Aligned length
  • (2) Alignment penalty per aligned length
• Sigalign returns the result of sequences satisfying the similarity criterion above without exception.
• Sigalign has only five parameters that affect the alignment result.
  • Three are penalty values.
    1. Mismatch penalty
    2. Gap open penalty
    3. Gap extend penalty
  • Two are cutoff for similarity values, which are aligned length and alignment penalty per aligned length.
    1. Minimum aligned length
    2. Maximum penalty per length
• A reference, which is an alignment target, has two main features.
  • A search range can be specified after reference generation.
  • A method of storing sequence data can be selected in various ways depending on the purpose.
    • The SequenceProvider storing the sequence is defined as an interface (trait in Rust).
• Sigalign provides two alignment algorithms
  • Semi-global alignment
    • Return the results of alignment that either the target or query sequence extended to the end from the left and right sides of the sequence respectively.
  • Local alignment
    • Return the longest alignment among the alignment satisfying the similarity cutoff is output as a result.

Quick Start

```rust use sigalign::{Reference, Aligner}; use sigalign::basicsequenceprovider::InMemoryProvider;

// (1) Make Reference let mut sequenceprovider = InMemoryProvider::newempty(); sequenceprovider.addlabeledsequence( "record1".tostring(), b"ATCAAACTCACAATTGTATTTCTTTGCCAGCTGGGCATATACTTTTTCCGCACCCTCATTTAACTTCTTGGATAACGGAAGCACACCGATCTTAACCGGAGCAAGTGCCGGATGAAAATGGAAAACGGTTCTTACGTCCGGCTTTTCCTCTGTTCCGATATTTTCCTCATCGTATGCAGCACATAAAAATGCCAGAACCA".tovec(), ); sequenceprovider.addlabeledsequence( "record2".tostring(), b"TTCCATCAAACTCACAATTGTATTTCTTTGCCAGCTGGGCATATACTTTTTCCGCACCCTCATTTAACTTCTTGGATAACGGAAGCACACCGATCTTAACCGGAGCGTATGCAGCACATAAAAAT".tovec(), ); let mut reference = Reference::newwithdefaultconfig(sequenceprovider).unwrap();

// (2) Make Aligner let aligner = Aligner::new(4, 6, 2, 100, 0.1).unwrap();

// (3) Alignment with query let query = b"TTCCTCTGTCATCAAACTCACAATTGTATTTCTTTGCCAGCTGGGCATATACTTTTTCCGCCCCCTCATTTAACTTCTTGGATAACGGAAGCACACCGATCTTAACCGGAGGTGCCGGATGAAAATGGAAAACGGTTCTTACGTCCGGCTTTTCCTCTGTTCCGATATTTTCCTCAT"; // - Semi-global alignment let resultsemiglobal: String = aligner.semiglobalalignmentlabeled(&mut reference, query).unwrap(); // - Local alignment let resultlocal: String = aligner.localalignmentlabeled(&mut reference, query).unwrap(); ```