1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
#!/usr/bin/env ruby
#
# gb2fasta.rb - convert GenBank entry into FASTA format (nuc)
#
# Copyright (C) 2001 KATAYAMA Toshiaki <k@bioruby.org>
# Copyright (C) 2002 Yoshinori K. Okuji <o@bioruby.org>
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# $Id: gb2fasta.rb,v 0.6 2008/02/05 12:11:11 pjotr Exp $
#
require 'bio'
include Bio
ff = FlatFile.new(GenBank, ARGF)
while gb = ff.next_entry
print gb.seq.to_fasta("gb:#{gb.entry_id} #{gb.definition}", 70)
end
|
