transrate 1.0.0 → 1.0.1
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/.travis.yml +17 -0
- data/README.md +1 -1
- data/Rakefile +1 -1
- data/bin/transrate +2 -375
- data/docs/transrate_logo_full.png +0 -0
- data/files.txt +78 -0
- data/lib/transrate.rb +1 -0
- data/lib/transrate/assembly.rb +4 -0
- data/lib/transrate/cmdline.rb +543 -0
- data/lib/transrate/comparative_metrics.rb +3 -0
- data/lib/transrate/snap.rb +10 -0
- data/lib/transrate/version.rb +1 -1
- data/packaging/build_deps_linux.sh +57 -0
- data/packaging/build_deps_macosx.sh +55 -0
- data/packaging/bundler-config +3 -0
- data/packaging/minify.sh +43 -0
- data/packaging/transrate +20 -0
- data/test/data/test_contig_nc1.fa +2 -0
- data/test/data/test_contig_nc2.fa +4 -0
- data/test/data/test_contig_nc3.fa +6 -0
- data/test/data/test_contig_nc4.fa +4 -0
- data/test/data/test_contig_nc5.fa +6 -0
- data/test/data/test_contig_nc6.fa +2 -0
- data/test/data/test_contig_nc7.fa +2 -0
- data/test/data/test_reference_aa1.fa +2 -0
- data/test/data/test_reference_nc1.fa +2 -0
- data/test/helper.rb +35 -6
- data/test/test_assembly.rb +5 -3
- data/test/test_cmd.rb +1 -1
- data/test/test_cmdline.rb +114 -0
- data/test/test_comp_metrics.rb +2 -2
- data/test/test_contig.rb +1 -1
- data/test/test_contig_metrics.rb +1 -1
- data/test/test_inline.rb +1 -1
- data/test/test_optimiser.rb +16 -16
- data/test/test_read_metrics.rb +1 -1
- data/test/test_salmon.rb +1 -1
- data/test/test_snap.rb +1 -1
- data/test/test_transrate.rb +1 -1
- data/test/test_transrater.rb +2 -2
- data/test/vagrant/centos_6.5_64/Vagrantfile +122 -0
- data/test/vagrant/debian_7.4_64/Vagrantfile +126 -0
- data/test/vagrant/debian_7.4_64/provision.sh +28 -0
- data/test/vagrant/fedora_20_64/Vagrantfile +122 -0
- data/test/vagrant/fedora_20_64/provision.sh +16 -0
- data/test/vagrant/linux-x86_64/Vagrantfile +10 -0
- data/test/vagrant/osx/Vagrantfile +18 -0
- data/test/vagrant/ubuntu_12.04_64/Vagrantfile +126 -0
- data/test/vagrant/ubuntu_12.04_64/provision.sh +24 -0
- data/transrate.gemspec +2 -2
- metadata +38 -17
- data/test/test_bin.rb +0 -139
@@ -0,0 +1,16 @@
|
|
1
|
+
# Setup
|
2
|
+
sudo yum install curl
|
3
|
+
|
4
|
+
# Install a multi-user RVM install with the latest ruby`
|
5
|
+
\curl -sSL https://get.rvm.io | sudo bash -s stable --ruby
|
6
|
+
source /etc/profile.d/rvm.sh
|
7
|
+
|
8
|
+
# Add user to RVM group
|
9
|
+
sudo usermod -a -G rvm vagrant
|
10
|
+
su vagrant
|
11
|
+
|
12
|
+
# Install transrate
|
13
|
+
gem install transrate --pre
|
14
|
+
|
15
|
+
# Install dependencies
|
16
|
+
transrate --install-deps
|
@@ -0,0 +1,10 @@
|
|
1
|
+
# -*- mode: ruby -*-
|
2
|
+
# vi: set ft=ruby :
|
3
|
+
|
4
|
+
# Vagrantfile API/syntax version. Don't touch unless you know what you're doing!
|
5
|
+
VAGRANTFILE_API_VERSION = "2"
|
6
|
+
|
7
|
+
Vagrant.configure(VAGRANTFILE_API_VERSION) do |config|
|
8
|
+
config.vm.box = "andy_b_84/debian-squeeze64-vanilla"
|
9
|
+
config.vm.network :forwarded_port, guest: 22, host: 18723
|
10
|
+
end
|
@@ -0,0 +1,18 @@
|
|
1
|
+
# -*- mode: ruby -*-
|
2
|
+
# vi: set ft=ruby :
|
3
|
+
|
4
|
+
# Vagrantfile API/syntax version. Don't touch unless you know what you're doing!
|
5
|
+
VAGRANTFILE_API_VERSION = "2"
|
6
|
+
|
7
|
+
Vagrant.configure(VAGRANTFILE_API_VERSION) do |config|
|
8
|
+
# Every Vagrant virtual environment requires a box to build off of.
|
9
|
+
config.vm.box = "http://files.dryga.com/boxes/osx-yosemite-0.2.0.box"
|
10
|
+
|
11
|
+
# config.vm.synced_folder "../data", "/vagrant_data"
|
12
|
+
config.vm.provider "virtualbox" do |vb|
|
13
|
+
config.vm.synced_folder ".", "/vagrant", type: "nfs"
|
14
|
+
end
|
15
|
+
|
16
|
+
config.vm.provision "shell", path: "../../../packaging/build_deps_macosx.sh"
|
17
|
+
|
18
|
+
end
|
@@ -0,0 +1,126 @@
|
|
1
|
+
# -*- mode: ruby -*-
|
2
|
+
# vi: set ft=ruby :
|
3
|
+
|
4
|
+
# Vagrantfile API/syntax version. Don't touch unless you know what you're doing!
|
5
|
+
VAGRANTFILE_API_VERSION = "2"
|
6
|
+
|
7
|
+
Vagrant.configure(VAGRANTFILE_API_VERSION) do |config|
|
8
|
+
# All Vagrant configuration is done here. The most common configuration
|
9
|
+
# options are documented and commented below. For a complete reference,
|
10
|
+
# please see the online documentation at vagrantup.com.
|
11
|
+
|
12
|
+
# Every Vagrant virtual environment requires a box to build off of.
|
13
|
+
config.vm.box = "hashicorp/precise64"
|
14
|
+
|
15
|
+
# Provision the virtual machine with a local script, provision.sh,
|
16
|
+
# in the same directory as the Vagrantfile
|
17
|
+
config.vm.provision "shell", path: "provision.sh"
|
18
|
+
|
19
|
+
# Disable automatic box update checking. If you disable this, then
|
20
|
+
# boxes will only be checked for updates when the user runs
|
21
|
+
# `vagrant box outdated`. This is not recommended.
|
22
|
+
# config.vm.box_check_update = false
|
23
|
+
|
24
|
+
# Create a forwarded port mapping which allows access to a specific port
|
25
|
+
# within the machine from a port on the host machine. In the example below,
|
26
|
+
# accessing "localhost:8080" will access port 80 on the guest machine.
|
27
|
+
# config.vm.network "forwarded_port", guest: 80, host: 8080
|
28
|
+
|
29
|
+
# Create a private network, which allows host-only access to the machine
|
30
|
+
# using a specific IP.
|
31
|
+
# config.vm.network "private_network", ip: "192.168.33.10"
|
32
|
+
|
33
|
+
# Create a public network, which generally matched to bridged network.
|
34
|
+
# Bridged networks make the machine appear as another physical device on
|
35
|
+
# your network.
|
36
|
+
# config.vm.network "public_network"
|
37
|
+
|
38
|
+
# If true, then any SSH connections made will enable agent forwarding.
|
39
|
+
# Default value: false
|
40
|
+
# config.ssh.forward_agent = true
|
41
|
+
|
42
|
+
# Share an additional folder to the guest VM. The first argument is
|
43
|
+
# the path on the host to the actual folder. The second argument is
|
44
|
+
# the path on the guest to mount the folder. And the optional third
|
45
|
+
# argument is a set of non-required options.
|
46
|
+
# config.vm.synced_folder "../data", "/vagrant_data"
|
47
|
+
|
48
|
+
# Provider-specific configuration so you can fine-tune various
|
49
|
+
# backing providers for Vagrant. These expose provider-specific options.
|
50
|
+
# Example for VirtualBox:
|
51
|
+
#
|
52
|
+
# config.vm.provider "virtualbox" do |vb|
|
53
|
+
# # Don't boot with headless mode
|
54
|
+
# vb.gui = true
|
55
|
+
#
|
56
|
+
# # Use VBoxManage to customize the VM. For example to change memory:
|
57
|
+
# vb.customize ["modifyvm", :id, "--memory", "1024"]
|
58
|
+
# end
|
59
|
+
#
|
60
|
+
# View the documentation for the provider you're using for more
|
61
|
+
# information on available options.
|
62
|
+
|
63
|
+
# Enable provisioning with CFEngine. CFEngine Community packages are
|
64
|
+
# automatically installed. For example, configure the host as a
|
65
|
+
# policy server and optionally a policy file to run:
|
66
|
+
#
|
67
|
+
# config.vm.provision "cfengine" do |cf|
|
68
|
+
# cf.am_policy_hub = true
|
69
|
+
# # cf.run_file = "motd.cf"
|
70
|
+
# end
|
71
|
+
#
|
72
|
+
# You can also configure and bootstrap a client to an existing
|
73
|
+
# policy server:
|
74
|
+
#
|
75
|
+
# config.vm.provision "cfengine" do |cf|
|
76
|
+
# cf.policy_server_address = "10.0.2.15"
|
77
|
+
# end
|
78
|
+
|
79
|
+
# Enable provisioning with Puppet stand alone. Puppet manifests
|
80
|
+
# are contained in a directory path relative to this Vagrantfile.
|
81
|
+
# You will need to create the manifests directory and a manifest in
|
82
|
+
# the file default.pp in the manifests_path directory.
|
83
|
+
#
|
84
|
+
# config.vm.provision "puppet" do |puppet|
|
85
|
+
# puppet.manifests_path = "manifests"
|
86
|
+
# puppet.manifest_file = "site.pp"
|
87
|
+
# end
|
88
|
+
|
89
|
+
# Enable provisioning with chef solo, specifying a cookbooks path, roles
|
90
|
+
# path, and data_bags path (all relative to this Vagrantfile), and adding
|
91
|
+
# some recipes and/or roles.
|
92
|
+
#
|
93
|
+
# config.vm.provision "chef_solo" do |chef|
|
94
|
+
# chef.cookbooks_path = "../my-recipes/cookbooks"
|
95
|
+
# chef.roles_path = "../my-recipes/roles"
|
96
|
+
# chef.data_bags_path = "../my-recipes/data_bags"
|
97
|
+
# chef.add_recipe "mysql"
|
98
|
+
# chef.add_role "web"
|
99
|
+
#
|
100
|
+
# # You may also specify custom JSON attributes:
|
101
|
+
# chef.json = { mysql_password: "foo" }
|
102
|
+
# end
|
103
|
+
|
104
|
+
# Enable provisioning with chef server, specifying the chef server URL,
|
105
|
+
# and the path to the validation key (relative to this Vagrantfile).
|
106
|
+
#
|
107
|
+
# The Opscode Platform uses HTTPS. Substitute your organization for
|
108
|
+
# ORGNAME in the URL and validation key.
|
109
|
+
#
|
110
|
+
# If you have your own Chef Server, use the appropriate URL, which may be
|
111
|
+
# HTTP instead of HTTPS depending on your configuration. Also change the
|
112
|
+
# validation key to validation.pem.
|
113
|
+
#
|
114
|
+
# config.vm.provision "chef_client" do |chef|
|
115
|
+
# chef.chef_server_url = "https://api.opscode.com/organizations/ORGNAME"
|
116
|
+
# chef.validation_key_path = "ORGNAME-validator.pem"
|
117
|
+
# end
|
118
|
+
#
|
119
|
+
# If you're using the Opscode platform, your validator client is
|
120
|
+
# ORGNAME-validator, replacing ORGNAME with your organization name.
|
121
|
+
#
|
122
|
+
# If you have your own Chef Server, the default validation client name is
|
123
|
+
# chef-validator, unless you changed the configuration.
|
124
|
+
#
|
125
|
+
# chef.validation_client_name = "ORGNAME-validator"
|
126
|
+
end
|
@@ -0,0 +1,24 @@
|
|
1
|
+
su vagrant
|
2
|
+
|
3
|
+
# Setup
|
4
|
+
sudo apt-get update -y
|
5
|
+
sudo apt-get install -y curl build-essential git
|
6
|
+
|
7
|
+
# Install a multi-user RVM install with the latest ruby
|
8
|
+
sudo gpg --keyserver hkp://keys.gnupg.net --recv-keys D39DC0E3
|
9
|
+
\curl -sSL https://get.rvm.io | bash -s stable --ruby
|
10
|
+
source /usr/local/rvm/scripts/rvm
|
11
|
+
|
12
|
+
# Add user to RVM group
|
13
|
+
sudo usermod -a -G rvm vagrant
|
14
|
+
|
15
|
+
# Install transrate
|
16
|
+
git clone https://github.com/Blahah/transrate.git
|
17
|
+
cd transrate
|
18
|
+
sudo chmod -R 777 .
|
19
|
+
|
20
|
+
# Install dependencies
|
21
|
+
gem install bundler rake
|
22
|
+
bundle install
|
23
|
+
bundle exec rake compile
|
24
|
+
bundle exec bin/transrate --install-deps
|
data/transrate.gemspec
CHANGED
@@ -28,8 +28,8 @@ Gem::Specification.new do |gem|
|
|
28
28
|
gem.add_development_dependency 'test-unit', '~> 3.0'
|
29
29
|
gem.add_development_dependency 'rake', '~> 10.3', '>= 10.3.2'
|
30
30
|
gem.add_development_dependency 'rake-compiler', '~> 0.9', '>= 0.9.2'
|
31
|
-
gem.add_development_dependency '
|
32
|
-
gem.add_development_dependency 'minitest', '~>
|
31
|
+
gem.add_development_dependency 'minitest', '~> 5'
|
32
|
+
gem.add_development_dependency 'minitest-reporters', '~> 1', '>= 1.0.17'
|
33
33
|
gem.add_development_dependency 'simplecov', '~> 0.8', '>= 0.8.2'
|
34
34
|
gem.add_development_dependency 'shoulda', '~> 3.5', '>= 3.5.0'
|
35
35
|
gem.add_development_dependency 'coveralls', '~> 0.7', '>= 0.7.2'
|
metadata
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: transrate
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
version: 1.0.
|
4
|
+
version: 1.0.1
|
5
5
|
platform: ruby
|
6
6
|
authors:
|
7
7
|
- Richard Smith-Unna
|
@@ -12,7 +12,7 @@ authors:
|
|
12
12
|
autorequire:
|
13
13
|
bindir: bin
|
14
14
|
cert_chain: []
|
15
|
-
date: 2015-
|
15
|
+
date: 2015-07-07 00:00:00.000000000 Z
|
16
16
|
dependencies:
|
17
17
|
- !ruby/object:Gem::Dependency
|
18
18
|
name: yell
|
@@ -209,45 +209,39 @@ dependencies:
|
|
209
209
|
- !ruby/object:Gem::Version
|
210
210
|
version: 0.9.2
|
211
211
|
- !ruby/object:Gem::Dependency
|
212
|
-
name:
|
212
|
+
name: minitest
|
213
213
|
requirement: !ruby/object:Gem::Requirement
|
214
214
|
requirements:
|
215
215
|
- - "~>"
|
216
216
|
- !ruby/object:Gem::Version
|
217
|
-
version: '
|
218
|
-
- - ">="
|
219
|
-
- !ruby/object:Gem::Version
|
220
|
-
version: 0.9.7
|
217
|
+
version: '5'
|
221
218
|
type: :development
|
222
219
|
prerelease: false
|
223
220
|
version_requirements: !ruby/object:Gem::Requirement
|
224
221
|
requirements:
|
225
222
|
- - "~>"
|
226
223
|
- !ruby/object:Gem::Version
|
227
|
-
version: '
|
228
|
-
- - ">="
|
229
|
-
- !ruby/object:Gem::Version
|
230
|
-
version: 0.9.7
|
224
|
+
version: '5'
|
231
225
|
- !ruby/object:Gem::Dependency
|
232
|
-
name: minitest
|
226
|
+
name: minitest-reporters
|
233
227
|
requirement: !ruby/object:Gem::Requirement
|
234
228
|
requirements:
|
235
229
|
- - "~>"
|
236
230
|
- !ruby/object:Gem::Version
|
237
|
-
version: '
|
231
|
+
version: '1'
|
238
232
|
- - ">="
|
239
233
|
- !ruby/object:Gem::Version
|
240
|
-
version:
|
234
|
+
version: 1.0.17
|
241
235
|
type: :development
|
242
236
|
prerelease: false
|
243
237
|
version_requirements: !ruby/object:Gem::Requirement
|
244
238
|
requirements:
|
245
239
|
- - "~>"
|
246
240
|
- !ruby/object:Gem::Version
|
247
|
-
version: '
|
241
|
+
version: '1'
|
248
242
|
- - ">="
|
249
243
|
- !ruby/object:Gem::Version
|
250
|
-
version:
|
244
|
+
version: 1.0.17
|
251
245
|
- !ruby/object:Gem::Dependency
|
252
246
|
name: simplecov
|
253
247
|
requirement: !ruby/object:Gem::Requirement
|
@@ -318,6 +312,7 @@ extensions:
|
|
318
312
|
extra_rdoc_files: []
|
319
313
|
files:
|
320
314
|
- ".gitignore"
|
315
|
+
- ".travis.yml"
|
321
316
|
- CITATION
|
322
317
|
- Gemfile
|
323
318
|
- LICENSE
|
@@ -326,11 +321,14 @@ files:
|
|
326
321
|
- bin/transrate
|
327
322
|
- deps/blast.yaml
|
328
323
|
- deps/deps.yaml
|
324
|
+
- docs/transrate_logo_full.png
|
329
325
|
- ext/transrate/extconf.rb
|
330
326
|
- ext/transrate/transrate.c
|
327
|
+
- files.txt
|
331
328
|
- lib/transrate.rb
|
332
329
|
- lib/transrate/assembly.rb
|
333
330
|
- lib/transrate/cmd.rb
|
331
|
+
- lib/transrate/cmdline.rb
|
334
332
|
- lib/transrate/comparative_metrics.rb
|
335
333
|
- lib/transrate/contig.rb
|
336
334
|
- lib/transrate/contig_metrics.rb
|
@@ -341,6 +339,11 @@ files:
|
|
341
339
|
- lib/transrate/transrater.rb
|
342
340
|
- lib/transrate/version.rb
|
343
341
|
- lib/transrate/writer.rb
|
342
|
+
- packaging/build_deps_linux.sh
|
343
|
+
- packaging/build_deps_macosx.sh
|
344
|
+
- packaging/bundler-config
|
345
|
+
- packaging/minify.sh
|
346
|
+
- packaging/transrate
|
344
347
|
- test/data/150uncovered.l.fq
|
345
348
|
- test/data/150uncovered.r.fq
|
346
349
|
- test/data/Os.protein.2.fa
|
@@ -352,11 +355,20 @@ files:
|
|
352
355
|
- test/data/sorghum_100.fa
|
353
356
|
- test/data/sorghum_transcript.fa
|
354
357
|
- test/data/test.sf
|
358
|
+
- test/data/test_contig_nc1.fa
|
359
|
+
- test/data/test_contig_nc2.fa
|
360
|
+
- test/data/test_contig_nc3.fa
|
361
|
+
- test/data/test_contig_nc4.fa
|
362
|
+
- test/data/test_contig_nc5.fa
|
363
|
+
- test/data/test_contig_nc6.fa
|
364
|
+
- test/data/test_contig_nc7.fa
|
365
|
+
- test/data/test_reference_aa1.fa
|
366
|
+
- test/data/test_reference_nc1.fa
|
355
367
|
- test/data/tiny.sam
|
356
368
|
- test/helper.rb
|
357
369
|
- test/test_assembly.rb
|
358
|
-
- test/test_bin.rb
|
359
370
|
- test/test_cmd.rb
|
371
|
+
- test/test_cmdline.rb
|
360
372
|
- test/test_comp_metrics.rb
|
361
373
|
- test/test_contig.rb
|
362
374
|
- test/test_contig_metrics.rb
|
@@ -367,6 +379,15 @@ files:
|
|
367
379
|
- test/test_snap.rb
|
368
380
|
- test/test_transrate.rb
|
369
381
|
- test/test_transrater.rb
|
382
|
+
- test/vagrant/centos_6.5_64/Vagrantfile
|
383
|
+
- test/vagrant/debian_7.4_64/Vagrantfile
|
384
|
+
- test/vagrant/debian_7.4_64/provision.sh
|
385
|
+
- test/vagrant/fedora_20_64/Vagrantfile
|
386
|
+
- test/vagrant/fedora_20_64/provision.sh
|
387
|
+
- test/vagrant/linux-x86_64/Vagrantfile
|
388
|
+
- test/vagrant/osx/Vagrantfile
|
389
|
+
- test/vagrant/ubuntu_12.04_64/Vagrantfile
|
390
|
+
- test/vagrant/ubuntu_12.04_64/provision.sh
|
370
391
|
- transrate.gemspec
|
371
392
|
homepage: https://github.com/Blahah/transrate
|
372
393
|
licenses:
|
data/test/test_bin.rb
DELETED
@@ -1,139 +0,0 @@
|
|
1
|
-
require 'helper'
|
2
|
-
require 'csv'
|
3
|
-
|
4
|
-
class TestTransrateBin < Test::Unit::TestCase
|
5
|
-
|
6
|
-
context "Transrate" do
|
7
|
-
|
8
|
-
setup do
|
9
|
-
|
10
|
-
end
|
11
|
-
|
12
|
-
teardown do
|
13
|
-
files = ["sorghum_100.1.fastq.sorghum_100.2.fastq.sorghum_100.bam.bai",
|
14
|
-
"sorghum_100.1.fastq.sorghum_100.2.fastq.sorghum_100.bam",
|
15
|
-
"sorghum_100_into_Os.protein.2.1.blast",
|
16
|
-
"assembly.2_into_Os.protein.2.1.blast",
|
17
|
-
"sorghum_100.nhr", "sorghum_100.nin", "sorghum_100.nsq",
|
18
|
-
"Os.protein.2_into_sorghum_100.2.blast",
|
19
|
-
"Os.protein.2_into_assembly.2.2.blast",
|
20
|
-
"assembly.2.nhr", "assembly.2.nin", "assembly.2.nsq",
|
21
|
-
"Os.protein.2.phr", "Os.protein.2.pin", "Os.protein.2.psq",
|
22
|
-
"transrate_assemblies.csv",
|
23
|
-
"sorghum_100.fa_quant.sf",
|
24
|
-
"sorghum_100.fa_bam_info.csv",
|
25
|
-
"transrate_sorghum_100.fa_contigs.csv",
|
26
|
-
"sorghum_100.1.fastq-sorghum_100.2.fastq-read_count.txt",
|
27
|
-
"sorghum_100.1.fastq.sorghum_100.2.fastq.sorghum_100.assigned.bam",
|
28
|
-
"bad.sorghum_100.fa", "chimeric.sorghum_100.fa",
|
29
|
-
"fragmented.sorghum_100.fa", "good.sorghum_100.fa"]
|
30
|
-
files.each do |file|
|
31
|
-
File.delete(file) if File.exist?(file)
|
32
|
-
end
|
33
|
-
`rm -rf sorghum_100`
|
34
|
-
`rm -rf logs`
|
35
|
-
`rm -rf libParams`
|
36
|
-
end
|
37
|
-
|
38
|
-
should "run help" do
|
39
|
-
c=Transrate::Cmd.new("bundle exec bin/transrate --help")
|
40
|
-
c.run
|
41
|
-
assert c.stdout =~ /DESCRIPTION/
|
42
|
-
assert_equal true, c.status.success?, "exit status"
|
43
|
-
end
|
44
|
-
|
45
|
-
should "fail on non existent assembly files" do
|
46
|
-
c=Transrate::Cmd.new("bundle exec bin/transrate --assembly foo.fasta")
|
47
|
-
c.run
|
48
|
-
assert_equal false, c.status.success?, "exit success"
|
49
|
-
end
|
50
|
-
|
51
|
-
should "fail on non existent reference files" do
|
52
|
-
c=Transrate::Cmd.new("bundle exec bin/transrate --reference foo.fasta")
|
53
|
-
c.run
|
54
|
-
assert_equal false, c.status.success?, "exit status"
|
55
|
-
end
|
56
|
-
|
57
|
-
should "run on test data" do
|
58
|
-
assembly = File.join(File.dirname(__FILE__), 'data', 'sorghum_100.fa')
|
59
|
-
reference = File.join(File.dirname(__FILE__), 'data', 'Os.protein.2.fa')
|
60
|
-
left = File.join(File.dirname(__FILE__), 'data', 'sorghum_100.1.fastq')
|
61
|
-
right = File.join(File.dirname(__FILE__), 'data', 'sorghum_100.2.fastq')
|
62
|
-
cmd = "bundle exec bin/transrate --assembly #{assembly}"
|
63
|
-
cmd << " --reference #{reference}"
|
64
|
-
cmd << " --left #{left}"
|
65
|
-
cmd << " --right #{right}"
|
66
|
-
c = Transrate::Cmd.new("#{cmd}")
|
67
|
-
c.run
|
68
|
-
assert_equal true, c.status.success?, "exit status"
|
69
|
-
assert File.exist?("transrate_assemblies.csv"), "csv file doesn't exist"
|
70
|
-
assert File.exist?("transrate_sorghum_100.fa_contigs.csv"),
|
71
|
-
"contig csv file doesn't exist"
|
72
|
-
hash = {}
|
73
|
-
CSV.foreach("transrate_assemblies.csv", :headers => true,
|
74
|
-
:header_converters => :symbol,
|
75
|
-
:converters => :all) do |row|
|
76
|
-
row.headers
|
77
|
-
row.fields
|
78
|
-
row.headers.zip(row.fields).each do |header, field|
|
79
|
-
hash[header]=field
|
80
|
-
end
|
81
|
-
end
|
82
|
-
assert_in_delta 137748, hash[:n_bases], 1000, "number of bases"
|
83
|
-
assert_equal 1692, hash[:n50], "n50"
|
84
|
-
assert_equal 2, hash[:n_refs_with_crbb], "number of crb hits"
|
85
|
-
assert_equal 2, hash[:n_contigs_with_crbb], "number of contigs with hits"
|
86
|
-
end
|
87
|
-
|
88
|
-
should "run on test data with comma separated list of fastq files" do
|
89
|
-
assembly = File.join(File.dirname(__FILE__), 'data', 'sorghum_100.fa')
|
90
|
-
left = []
|
91
|
-
right = []
|
92
|
-
left << File.join(File.dirname(__FILE__), 'data', 'sorghum_100.1.fastq')
|
93
|
-
left << File.join(File.dirname(__FILE__), 'data', '150uncovered.l.fq')
|
94
|
-
right << File.join(File.dirname(__FILE__), 'data', 'sorghum_100.2.fastq')
|
95
|
-
right << File.join(File.dirname(__FILE__), 'data', '150uncovered.r.fq')
|
96
|
-
cmd = "bundle exec bin/transrate --assembly #{assembly}"
|
97
|
-
cmd << " --left #{left.join(",")}"
|
98
|
-
cmd << " --right #{right.join(",")}"
|
99
|
-
c = Transrate::Cmd.new("#{cmd}")
|
100
|
-
c.run
|
101
|
-
if !(c.status.success?)
|
102
|
-
puts c.stdout
|
103
|
-
puts c.stderr
|
104
|
-
end
|
105
|
-
assert_equal true, c.status.success?, "exit status"
|
106
|
-
assert File.exist?("transrate_assemblies.csv")
|
107
|
-
assert File.exist?("transrate_sorghum_100.fa_contigs.csv"),
|
108
|
-
"contig csv file doesn't exist"
|
109
|
-
hash = {}
|
110
|
-
CSV.foreach("transrate_assemblies.csv", :headers => true,
|
111
|
-
:header_converters => :symbol,
|
112
|
-
:converters => :all) do |row|
|
113
|
-
row.headers.zip(row.fields).each do |header, field|
|
114
|
-
hash[header]=field
|
115
|
-
end
|
116
|
-
end
|
117
|
-
assert_in_delta 137748, hash[:n_bases], 1000, "number of bases"
|
118
|
-
assert_equal 1692, hash[:n50], "n50"
|
119
|
-
assert_equal 25006 + 223, hash[:fragments], "number of reads"
|
120
|
-
end
|
121
|
-
|
122
|
-
should "fail when one of multiple assemblies is missing" do
|
123
|
-
assembly = File.join(File.dirname(__FILE__), 'data', 'assembly.2.fa')
|
124
|
-
assembly2 = File.join(File.dirname(__FILE__), 'data', 'foo.fa')
|
125
|
-
reference = File.join(File.dirname(__FILE__), 'data', 'Os.protein.2.fa')
|
126
|
-
left = File.join(File.dirname(__FILE__), 'data', '150uncovered.l.fq')
|
127
|
-
right = File.join(File.dirname(__FILE__), 'data', '150uncovered.r.fq')
|
128
|
-
cmd = "bundle exec bin/transrate "
|
129
|
-
cmd << " --assembly #{assembly},#{assembly2}"
|
130
|
-
cmd << " --reference #{reference}"
|
131
|
-
cmd << " --left #{left}"
|
132
|
-
cmd << " --right #{right}"
|
133
|
-
c = Transrate::Cmd.new("#{cmd}")
|
134
|
-
c.run
|
135
|
-
assert_equal false, c.status.success?, "exit status"
|
136
|
-
end
|
137
|
-
|
138
|
-
end
|
139
|
-
end
|