Edited wiki page PhyloBioRuby through web user interface.
authorcmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Sat, 26 Mar 2011 05:49:18 +0000 (05:49 +0000)
committercmzmasek@gmail.com <cmzmasek@gmail.com@ca865154-3058-d1c3-3e42-d8f55a55bdbd>
Sat, 26 Mar 2011 05:49:18 +0000 (05:49 +0000)
wiki/PhyloBioRuby.wiki

index 503fc75..ab6faf9 100644 (file)
@@ -43,6 +43,32 @@ msa.each do |entry|
 end
 }}}
 
+Blah
+{{{
+# Reads in a Fasta-formatted multiple sequence alignment (which does
+# not have to be aligned, though) and stores its sequences in
+# array 'seq_ary'.
+seq_ary = Array.new
+fasta_seqs = Bio::Alignment::MultiFastaFormat.new(File.open('bcl2.fasta').read)
+fasta_seqs.entries.each do |seq|
+  seq_ary.push( seq )
+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 }
+
+
+puts seqs.consensus
+
+# Writes multiple sequence alignment (possibly unaligned) 'seqs'
+# to a file in phylip format.
+File.open('out1.phylip', 'w') do |f|
+  f.write(seqs.output(:phylip))
+end
+}}}
+
 Relevant API documentation:
 
  * [http://bioruby.open-bio.org/rdoc/classes/Bio/ClustalW/Report.html Bio::ClustalW::Report]