bio-octopus 0.1.0 → 0.1.1
Sign up to get free protection for your applications and to get access to all the features.
- data/README.rdoc +9 -0
- data/VERSION +1 -1
- data/bin/biooctopus +21 -6
- data/bin/biooctopus.rb +21 -6
- data/bio-octopus.gemspec +2 -2
- metadata +4 -4
data/README.rdoc
CHANGED
@@ -2,6 +2,15 @@
|
|
2
2
|
|
3
3
|
A bioruby plugin for running and parsing octopus, available from http://octopus.cbr.su.se/ Currently, the plugin is not Ruby 1.9 (only 1.8) compatible, because of the dependency on the rio gem.
|
4
4
|
|
5
|
+
To run it as a script, a fasta file of sequences to predict, and a database formatted in legacy blast database format needs to be specified:
|
6
|
+
|
7
|
+
$ /usr/local/rvm/gems/ruby-1.8.7-p334/gems/bio-octopus-0.1.0/bin/biooctopus my.fasta /blastdb/UniProt2011_04/uniprot_sprot.fasta
|
8
|
+
|
9
|
+
Output are topology predictions, including signal peptide and transmembrane domain predictions. If that is not the best format for you, look at the script file and modify it - it is just a front-end. Also required are:
|
10
|
+
* a working install of SPOCTOPUS
|
11
|
+
* NCBI legacy blast
|
12
|
+
Binaries for these must be accessible in the PATH.
|
13
|
+
|
5
14
|
== Contributing to bio-octopus
|
6
15
|
|
7
16
|
* Check out the latest master to make sure the feature hasn't been implemented or the bug hasn't been fixed yet
|
data/VERSION
CHANGED
@@ -1 +1 @@
|
|
1
|
-
0.1.
|
1
|
+
0.1.1
|
data/bin/biooctopus
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
#!/usr/bin/env ruby
|
2
2
|
#
|
3
3
|
# BioRuby bio-octopus Plugin
|
4
|
-
# Version 0.
|
4
|
+
# Version 0.1.1
|
5
5
|
# Author:: Ben J Woodcroft
|
6
6
|
# Copyright:: 2011
|
7
7
|
# License:: The Ruby License
|
@@ -21,12 +21,12 @@ if $0 == __FILE__
|
|
21
21
|
].flatten
|
22
22
|
end
|
23
23
|
o.parse!
|
24
|
-
|
24
|
+
|
25
25
|
unless ARGV.length == 2
|
26
|
-
|
27
|
-
|
26
|
+
puts o.help
|
27
|
+
exit
|
28
28
|
end
|
29
|
-
|
29
|
+
|
30
30
|
fasta_path = ARGV[0]
|
31
31
|
blastdb_path = ARGV[1]
|
32
32
|
|
@@ -36,6 +36,21 @@ if $0 == __FILE__
|
|
36
36
|
result = runner.calculate(seq.seq, blastdb_path)
|
37
37
|
name = seq.definition
|
38
38
|
|
39
|
+
# Output singal peptide if it has one
|
40
|
+
if result.signal?
|
41
|
+
puts [
|
42
|
+
name,
|
43
|
+
'Signal peptide',
|
44
|
+
result.signal_peptide.start,
|
45
|
+
result.signal_peptide.stop
|
46
|
+
].join("\t")
|
47
|
+
else
|
48
|
+
puts [
|
49
|
+
name,
|
50
|
+
'No signal peptide found'
|
51
|
+
].join("\t")
|
52
|
+
end
|
53
|
+
|
39
54
|
if result.has_domain?
|
40
55
|
# At least one TMD found. Output each on a separate line
|
41
56
|
result.transmembrane_domains.each do |tmd|
|
@@ -50,7 +65,7 @@ if $0 == __FILE__
|
|
50
65
|
else
|
51
66
|
puts [
|
52
67
|
name,
|
53
|
-
'No
|
68
|
+
'No transmembrane domain found'
|
54
69
|
].join("\t")
|
55
70
|
end
|
56
71
|
end
|
data/bin/biooctopus.rb
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
#!/usr/bin/env ruby
|
2
2
|
#
|
3
3
|
# BioRuby bio-octopus Plugin
|
4
|
-
# Version 0.
|
4
|
+
# Version 0.1.1
|
5
5
|
# Author:: Ben J Woodcroft
|
6
6
|
# Copyright:: 2011
|
7
7
|
# License:: The Ruby License
|
@@ -21,12 +21,12 @@ if $0 == __FILE__
|
|
21
21
|
].flatten
|
22
22
|
end
|
23
23
|
o.parse!
|
24
|
-
|
24
|
+
|
25
25
|
unless ARGV.length == 2
|
26
|
-
|
27
|
-
|
26
|
+
puts o.help
|
27
|
+
exit
|
28
28
|
end
|
29
|
-
|
29
|
+
|
30
30
|
fasta_path = ARGV[0]
|
31
31
|
blastdb_path = ARGV[1]
|
32
32
|
|
@@ -36,6 +36,21 @@ if $0 == __FILE__
|
|
36
36
|
result = runner.calculate(seq.seq, blastdb_path)
|
37
37
|
name = seq.definition
|
38
38
|
|
39
|
+
# Output singal peptide if it has one
|
40
|
+
if result.signal?
|
41
|
+
puts [
|
42
|
+
name,
|
43
|
+
'Signal peptide',
|
44
|
+
result.signal_peptide.start,
|
45
|
+
result.signal_peptide.stop
|
46
|
+
].join("\t")
|
47
|
+
else
|
48
|
+
puts [
|
49
|
+
name,
|
50
|
+
'No signal peptide found'
|
51
|
+
].join("\t")
|
52
|
+
end
|
53
|
+
|
39
54
|
if result.has_domain?
|
40
55
|
# At least one TMD found. Output each on a separate line
|
41
56
|
result.transmembrane_domains.each do |tmd|
|
@@ -50,7 +65,7 @@ if $0 == __FILE__
|
|
50
65
|
else
|
51
66
|
puts [
|
52
67
|
name,
|
53
|
-
'No
|
68
|
+
'No transmembrane domain found'
|
54
69
|
].join("\t")
|
55
70
|
end
|
56
71
|
end
|
data/bio-octopus.gemspec
CHANGED
@@ -5,11 +5,11 @@
|
|
5
5
|
|
6
6
|
Gem::Specification.new do |s|
|
7
7
|
s.name = %q{bio-octopus}
|
8
|
-
s.version = "0.1.
|
8
|
+
s.version = "0.1.1"
|
9
9
|
|
10
10
|
s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
|
11
11
|
s.authors = ["Ben J Woodcroft"]
|
12
|
-
s.date = %q{2011-
|
12
|
+
s.date = %q{2011-05-02}
|
13
13
|
s.description = %q{Running and parsing of the protein transmembrane domain predictor octopus}
|
14
14
|
s.email = %q{gmail.com after donttrustben}
|
15
15
|
s.executables = ["biooctopus.rb", "biooctopus"]
|
metadata
CHANGED
@@ -1,13 +1,13 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: bio-octopus
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
hash:
|
4
|
+
hash: 25
|
5
5
|
prerelease:
|
6
6
|
segments:
|
7
7
|
- 0
|
8
8
|
- 1
|
9
|
-
-
|
10
|
-
version: 0.1.
|
9
|
+
- 1
|
10
|
+
version: 0.1.1
|
11
11
|
platform: ruby
|
12
12
|
authors:
|
13
13
|
- Ben J Woodcroft
|
@@ -15,7 +15,7 @@ autorequire:
|
|
15
15
|
bindir: bin
|
16
16
|
cert_chain: []
|
17
17
|
|
18
|
-
date: 2011-
|
18
|
+
date: 2011-05-02 00:00:00 +10:00
|
19
19
|
default_executable:
|
20
20
|
dependencies:
|
21
21
|
- !ruby/object:Gem::Dependency
|