summaryrefslogtreecommitdiffstats
path: root/GRPH.rs
diff options
context:
space:
mode:
Diffstat (limited to 'GRPH.rs')
-rw-r--r--GRPH.rs22
1 files changed, 22 insertions, 0 deletions
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));
+ }
+}