From 00e7efaee4bb250758563289ef509101b43715f1 Mon Sep 17 00:00:00 2001 From: Jesse Luehrs Date: Sat, 9 Mar 2013 00:09:38 -0600 Subject: grph --- GRPH.rs | 22 ++++++++++++++++++++++ 1 file changed, 22 insertions(+) create mode 100644 GRPH.rs diff --git a/GRPH.rs b/GRPH.rs new file mode 100644 index 0000000..867aade --- /dev/null +++ b/GRPH.rs @@ -0,0 +1,22 @@ +use str::{len,view}; + +extern mod rosalind; +use rosalind::fasta::FASTAReader; + +fn main() { + let mut dna: ~[(~str, ~str)] = ~[]; + let reader = FASTAReader::new(); + for reader.each_sequence |cur_name, cur_seq| { + let cur_len = len(cur_seq); + for dna.each |&(prev_name, prev_seq)| { + let prev_len = len(prev_seq); + if view(cur_seq, 0, 3) == view(prev_seq, prev_len - 3, prev_len) { + io::println(fmt!("%s %s", prev_name, cur_name)); + } + if view(prev_seq, 0, 3) == view(cur_seq, cur_len - 3, cur_len) { + io::println(fmt!("%s %s", cur_name, prev_name)); + } + } + dna.push((cur_name, cur_seq)); + } +} -- cgit v1.2.3