http://netvouz.com/aunderwo/tag/report.alignment?feed=rss aunderwo's bookmarks tagged "report.alignment" on Netvouz http://code.google.com/p/forester/wiki/PhyloBioRuby #!/usr/bin/env ruby require 'bio' # Reads in a ClustalW-formatted multiple sequence alignment # from a file named "infile_clustalw.aln" and stores it in 'report'. report = Bio::ClustalW::Report.new(File.read('infile_clustalw.aln')) # Accesses the actual alignment. msa = report.alignment # Goes through all sequences in 'msa' and prints the # actual molecular sequence. msa.each do |entry|   puts entry.seq end aunderwo Fri, 18 Mar 2011 09:35:43 GMT