編集距離k以内の塩基配列("ATGC"のみ許可)を深さ優先探索で列挙するRustコードを書いたので備忘録として残しておきます。もうちょっとコードをキレイにしたかったのとテストコードは宿題にします。 use bio::alphabets; #[derive(Clone, Copy)] enum Operations { Match, Subst(u8), Del, Ins(u8), } struct NeighborSeq { origin_seq: Vec<u8>, stack: Vec<StateNS>, curr: Vec<Operations>, alphabet: alphabets:…