bio-gadget 0.1.0
Sign up to get free protection for your applications and to get access to all the features.
- data/.gitignore +21 -0
- data/Gemfile +4 -0
- data/Gthorfile +46 -0
- data/LICENSE +22 -0
- data/README.org +20 -0
- data/Rakefile +2 -0
- data/bio-gadget.gemspec +19 -0
- data/lib/bio-gadget.rb +7 -0
- data/lib/bio-gadget/version.rb +5 -0
- metadata +70 -0
data/.gitignore
ADDED
data/Gemfile
ADDED
data/Gthorfile
ADDED
@@ -0,0 +1,46 @@
|
|
1
|
+
module Bio
|
2
|
+
class Gadget < Thor
|
3
|
+
|
4
|
+
namespace :bio
|
5
|
+
|
6
|
+
desc 'qvstat QUAL', 'statistics of quality values in *.qual file'
|
7
|
+
def qvstat(qualfile)
|
8
|
+
stat = Hash.new
|
9
|
+
myopen(qualfile) { |fp|
|
10
|
+
fp.each { |line|
|
11
|
+
next if /^[\#\>]/ =~ line
|
12
|
+
qvs = line.rstrip.split
|
13
|
+
qvs.each_index { |i|
|
14
|
+
qv = qvs[i]
|
15
|
+
stat[qv] = Array.new unless stat.key?(qv)
|
16
|
+
if stat[qv][i].nil?
|
17
|
+
stat[qv][i] = 1
|
18
|
+
else
|
19
|
+
stat[qv][i] += 1
|
20
|
+
end
|
21
|
+
}
|
22
|
+
}
|
23
|
+
}
|
24
|
+
statfile = qualfile.sub(/.qual(.gz)?$/, '.qvstat')
|
25
|
+
open(statfile, 'w') { |out|
|
26
|
+
qvs = stat.keys.sort { |a, b| a.to_i <=> b.to_i }
|
27
|
+
qvs.each { |qv|
|
28
|
+
out.puts "#{qv} #{stat[qv].join(' ')}"
|
29
|
+
}
|
30
|
+
}
|
31
|
+
end
|
32
|
+
|
33
|
+
private
|
34
|
+
|
35
|
+
def myopen(file, &block)
|
36
|
+
# how to write?
|
37
|
+
f = (/\|/ !~ file && /\.gz$/ =~ file) ? "| gunzip -c #{file}" : file
|
38
|
+
unless block.nil?
|
39
|
+
o = open(f); block.call(o); o.close
|
40
|
+
else
|
41
|
+
open(f)
|
42
|
+
end
|
43
|
+
end
|
44
|
+
|
45
|
+
end
|
46
|
+
end
|
data/LICENSE
ADDED
@@ -0,0 +1,22 @@
|
|
1
|
+
Copyright (c) 2012 Shintaro Katayama
|
2
|
+
|
3
|
+
MIT License
|
4
|
+
|
5
|
+
Permission is hereby granted, free of charge, to any person obtaining
|
6
|
+
a copy of this software and associated documentation files (the
|
7
|
+
"Software"), to deal in the Software without restriction, including
|
8
|
+
without limitation the rights to use, copy, modify, merge, publish,
|
9
|
+
distribute, sublicense, and/or sell copies of the Software, and to
|
10
|
+
permit persons to whom the Software is furnished to do so, subject to
|
11
|
+
the following conditions:
|
12
|
+
|
13
|
+
The above copyright notice and this permission notice shall be
|
14
|
+
included in all copies or substantial portions of the Software.
|
15
|
+
|
16
|
+
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
|
17
|
+
EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
|
18
|
+
MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
|
19
|
+
NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
|
20
|
+
LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
|
21
|
+
OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
|
22
|
+
WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
|
data/README.org
ADDED
@@ -0,0 +1,20 @@
|
|
1
|
+
* Bio::Gadget
|
2
|
+
|
3
|
+
*Bio::Gadget* - Gadgets for bioinformatics
|
4
|
+
|
5
|
+
* Installation
|
6
|
+
|
7
|
+
: gem install bio-gadget
|
8
|
+
|
9
|
+
* Usage
|
10
|
+
|
11
|
+
: gthor help bio
|
12
|
+
: gthor list bio
|
13
|
+
|
14
|
+
* Contributing
|
15
|
+
|
16
|
+
1. Fork it
|
17
|
+
2. Create your feature branch (`git checkout -b my-new-feature`)
|
18
|
+
3. Commit your changes (`git commit -am 'Added some feature'`)
|
19
|
+
4. Push to the branch (`git push origin my-new-feature`)
|
20
|
+
5. Create new Pull Request
|
data/Rakefile
ADDED
data/bio-gadget.gemspec
ADDED
@@ -0,0 +1,19 @@
|
|
1
|
+
# -*- encoding: utf-8 -*-
|
2
|
+
require File.expand_path('../lib/bio-gadget/version', __FILE__)
|
3
|
+
|
4
|
+
Gem::Specification.new do |gem|
|
5
|
+
gem.authors = ["Shintaro Katayama"]
|
6
|
+
gem.email = ["shintaro.katayama@gmail.com"]
|
7
|
+
gem.description = %q{Gadgets for bioinformatics}
|
8
|
+
gem.summary = gem.description
|
9
|
+
gem.homepage = "https://github.com/shka/ruby-bio-gadget"
|
10
|
+
|
11
|
+
gem.files = `git ls-files`.split($\)
|
12
|
+
gem.executables = gem.files.grep(%r{^bin/}).map{ |f| File.basename(f) }
|
13
|
+
gem.test_files = gem.files.grep(%r{^(test|spec|features)/})
|
14
|
+
gem.name = "bio-gadget"
|
15
|
+
gem.require_paths = ["lib"]
|
16
|
+
gem.version = Bio::Gadget::VERSION
|
17
|
+
|
18
|
+
gem.add_dependency 'gthor'
|
19
|
+
end
|
data/lib/bio-gadget.rb
ADDED
metadata
ADDED
@@ -0,0 +1,70 @@
|
|
1
|
+
--- !ruby/object:Gem::Specification
|
2
|
+
name: bio-gadget
|
3
|
+
version: !ruby/object:Gem::Version
|
4
|
+
version: 0.1.0
|
5
|
+
prerelease:
|
6
|
+
platform: ruby
|
7
|
+
authors:
|
8
|
+
- Shintaro Katayama
|
9
|
+
autorequire:
|
10
|
+
bindir: bin
|
11
|
+
cert_chain: []
|
12
|
+
date: 2012-08-20 00:00:00.000000000 Z
|
13
|
+
dependencies:
|
14
|
+
- !ruby/object:Gem::Dependency
|
15
|
+
name: gthor
|
16
|
+
requirement: !ruby/object:Gem::Requirement
|
17
|
+
none: false
|
18
|
+
requirements:
|
19
|
+
- - ! '>='
|
20
|
+
- !ruby/object:Gem::Version
|
21
|
+
version: '0'
|
22
|
+
type: :runtime
|
23
|
+
prerelease: false
|
24
|
+
version_requirements: !ruby/object:Gem::Requirement
|
25
|
+
none: false
|
26
|
+
requirements:
|
27
|
+
- - ! '>='
|
28
|
+
- !ruby/object:Gem::Version
|
29
|
+
version: '0'
|
30
|
+
description: Gadgets for bioinformatics
|
31
|
+
email:
|
32
|
+
- shintaro.katayama@gmail.com
|
33
|
+
executables: []
|
34
|
+
extensions: []
|
35
|
+
extra_rdoc_files: []
|
36
|
+
files:
|
37
|
+
- .gitignore
|
38
|
+
- Gemfile
|
39
|
+
- Gthorfile
|
40
|
+
- LICENSE
|
41
|
+
- README.org
|
42
|
+
- Rakefile
|
43
|
+
- bio-gadget.gemspec
|
44
|
+
- lib/bio-gadget.rb
|
45
|
+
- lib/bio-gadget/version.rb
|
46
|
+
homepage: https://github.com/shka/ruby-bio-gadget
|
47
|
+
licenses: []
|
48
|
+
post_install_message:
|
49
|
+
rdoc_options: []
|
50
|
+
require_paths:
|
51
|
+
- lib
|
52
|
+
required_ruby_version: !ruby/object:Gem::Requirement
|
53
|
+
none: false
|
54
|
+
requirements:
|
55
|
+
- - ! '>='
|
56
|
+
- !ruby/object:Gem::Version
|
57
|
+
version: '0'
|
58
|
+
required_rubygems_version: !ruby/object:Gem::Requirement
|
59
|
+
none: false
|
60
|
+
requirements:
|
61
|
+
- - ! '>='
|
62
|
+
- !ruby/object:Gem::Version
|
63
|
+
version: '0'
|
64
|
+
requirements: []
|
65
|
+
rubyforge_project:
|
66
|
+
rubygems_version: 1.8.23
|
67
|
+
signing_key:
|
68
|
+
specification_version: 3
|
69
|
+
summary: Gadgets for bioinformatics
|
70
|
+
test_files: []
|