X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=wiki%2FPhyloBioRuby.wiki;h=4b6b45ffa63c3c20de0a76453f618aa7b8714b94;hb=248d79b589b1f461b2ca2aa77f2d6457d3f3cadf;hp=595e20f52408c567a385981e5b97652f021b5e6b;hpb=bb1e65895bc0d93d9e90b511b78a29cd3573810f;p=jalview.git

diff --git a/wiki/PhyloBioRuby.wiki b/wiki/PhyloBioRuby.wiki
index 595e20f..4b6b45f 100644
--- a/wiki/PhyloBioRuby.wiki
+++ b/wiki/PhyloBioRuby.wiki
@@ -10,7 +10,7 @@ Tutorial for multiple sequence alignments and phylogenetic methods in [http://bi
 
 Eventually, this is expected to be placed on the official !BioRuby page.
 
-Author: [http://www.cmzmasek.net/ Christian M Zmasek], Sanford-Burnham Medical Research Institute
+Author: [https://sites.google.com/site/cmzmasek/ Christian Zmasek], Sanford-Burnham Medical Research Institute
 
  
 Copyright (C) 2011 Christian M Zmasek. All rights reserved.
@@ -23,6 +23,41 @@ Copyright (C) 2011 Christian M Zmasek. All rights reserved.
 
 === Reading in a Multiple Sequence Alignment from a File ===
 
+This automatically determines the format
+{{{
+#!/usr/bin/env ruby
+require 'bio'
+
+seq_ary = Array.new
+ff = Bio::FlatFile.auto('bcl2.fasta')
+ff.each_entry do |entry|
+  seq_ary.push(entry)
+  puts entry.entry_id          # prints the identifier of the entry
+  puts entry.definition        # prints the definition of the entry
+  puts entry.seq               # prints the sequence data of the entry
+end
+
+# Creates a multiple sequence alignment (possibly unaligned) named
+# 'seqs' from array 'seq_ary'.
+seqs = Bio::Alignment.new(seq_ary)
+seqs.each { |seq| puts seq.to_s }
+
+# Writes multiple sequence alignment (possibly unaligned) 'seqs'
+# to a file in PHYLIP format.
+File.open('out0.phylip', 'w') do |f|
+  f.write(seqs.output(:phylip))
+end
+
+# Writes multiple sequence alignment (possibly unaligned) 'seqs'
+# to a file in FASTA format.
+File.open('out0.fasta', 'w') do |f|
+  f.write(seqs.output(:fasta))
+end
+}}}
+
+
+==== ClustalW Format ====
+
 The following example shows how to read in a *ClustalW*-formatted multiple sequence alignment.
 
 {{{
@@ -43,6 +78,7 @@ msa.each do |entry|
 end
 }}}
 
+==== FASTA Format ====
 
 The following example shows how to read in a *FASTA*-formatted multiple sequence file. (_This seems a little clumsy, I wonder if there is a more direct way, avoiding the creation of an array.)
 {{{
@@ -67,7 +103,7 @@ seqs.each { |seq| puts seq.to_s }
 
 # Writes multiple sequence alignment (possibly unaligned) 'seqs'
 # to a file in PHYLIP format.
-File.open('out1.phylip', 'w') do |f|
+File.open('outfile.phylip', 'w') do |f|
   f.write(seqs.output(:phylip))
 end
 }}}
@@ -78,6 +114,13 @@ Relevant API documentation:
  * [http://bioruby.open-bio.org/rdoc/classes/Bio/Alignment.html Bio::Alignment]
  * [http://bioruby.open-bio.org/rdoc/classes/Bio/Sequence.html Bio::Sequence]
 
+=== Creating a Multiple Sequence Alignment ===
+
+
+=== Creating a Multiple Sequence Alignment from a Database ===
+
+?
+
 === Writing a Multiple Sequence Alignment to a File ===