lederhosen 1.8.1 → 1.8.2

Sign up to get free protection for your applications and to get access to all the features.
data/lederhosen.gemspec CHANGED
@@ -5,7 +5,7 @@
5
5
 
6
6
  Gem::Specification.new do |s|
7
7
  s.name = "lederhosen"
8
- s.version = "1.8.1"
8
+ s.version = "1.8.2"
9
9
 
10
10
  s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
11
11
  s.authors = ["Austin G. Davis-Richardson"]
@@ -43,16 +43,10 @@ Gem::Specification.new do |s|
43
43
  "scripts/otu_ref_picking/readme.md",
44
44
  "scripts/readme.md",
45
45
  "spec/cli_spec.rb",
46
- "spec/data/ILT_L_9_B_001_1.txt.gz",
47
- "spec/data/ILT_L_9_B_001_3.txt.gz",
48
- "spec/data/ILT_L_9_B_002_1.txt.gz",
49
- "spec/data/ILT_L_9_B_002_3.txt.gz",
50
- "spec/data/example.fastq",
51
46
  "spec/data/test.uc",
52
- "spec/data/trimmed_sizes.txt",
47
+ "spec/data/trimmed/ILT_L_9_B_001.fasta",
53
48
  "spec/no_tasks_spec.rb",
54
- "spec/spec_helper.rb",
55
- "spec/trimmer_spec.rb"
49
+ "spec/spec_helper.rb"
56
50
  ]
57
51
  s.homepage = "http://audy.github.com/lederhosen"
58
52
  s.licenses = ["MIT"]
@@ -54,6 +54,9 @@ module Lederhosen
54
54
  else
55
55
  dat[level] || 'unparsed_name'
56
56
  end
57
+
58
+ # remove commas from name
59
+ name = name.tr(',', '_')
57
60
 
58
61
  # the next two lines are what is slow
59
62
  level_sample_cluster_count[level][input_file][name] += 1
@@ -3,7 +3,7 @@ module Lederhosen
3
3
  MAJOR = 1
4
4
  MINOR = 8
5
5
  CODENAME = 'Karottensaft' # changes for minor versions
6
- PATCH = 1
6
+ PATCH = 2
7
7
 
8
8
  STRING = [MAJOR, MINOR, PATCH].join('.')
9
9
  end
data/spec/cli_spec.rb CHANGED
@@ -12,30 +12,20 @@ describe Lederhosen::CLI do
12
12
  $?.success?.should be_true
13
13
  end
14
14
 
15
- it 'should trim paired QSEQ reads' do
16
- `./bin/lederhosen trim --reads-dir=spec/data/IL*.txt.gz --out-dir=#{$test_dir}/trimmed`
17
- $?.success?.should be_true
18
- end
19
-
20
- it 'should trim paired, interleaved FASTQ reads' do
21
- `./bin/lederhosen trim --reads-dir=spec/data/example.fastq --out-dir=#{$test_dir}/trimmed --read-type fastq`
22
- $?.success?.should be_true
23
- end
24
-
25
15
  it 'can create a usearch udb using usearch' do
26
- `./bin/lederhosen make_udb --input #{$test_dir}/trimmed/ILT_L_9_B_001.fasta --output #{$test_dir}/test_db.udb`
16
+ `./bin/lederhosen make_udb --input spec/data/trimmed/ILT_L_9_B_001.fasta --output #{$test_dir}/test_db.udb`
27
17
  $?.success?.should be_true
28
18
  end
29
19
 
30
20
  it 'can (dry run) cluster reads using usearch and output usearch cmd to stdout' do
31
- stdout = `./bin/lederhosen cluster --dry-run --input #{$test_dir}/trimmed/ILT_L_9_B_001.fasta --database #{$test_dir}/test_db.udb --identity 0.95 --output #{$test_dir}/clusters.uc`
21
+ stdout = `./bin/lederhosen cluster --dry-run --input spec/data/trimmed/ILT_L_9_B_001.fasta --database #{$test_dir}/test_db.udb --identity 0.95 --output #{$test_dir}/clusters.uc`
32
22
  stdout.should match /^usearch/
33
23
  $?.success?.should be_true
34
24
  File.exists?(File.join($test_dir, 'clusters.uc')).should be_false
35
25
  end
36
26
 
37
27
  it 'can cluster reads using usearch' do
38
- `./bin/lederhosen cluster --input #{$test_dir}/trimmed/ILT_L_9_B_001.fasta --database #{$test_dir}/test_db.udb --identity 0.95 --output #{$test_dir}/clusters.uc`
28
+ `./bin/lederhosen cluster --input spec/data/trimmed/ILT_L_9_B_001.fasta --database #{$test_dir}/test_db.udb --identity 0.95 --output #{$test_dir}/clusters.uc`
39
29
  $?.success?.should be_true
40
30
  File.exists?(File.join($test_dir, 'clusters.uc')).should be_true
41
31
  end
@@ -57,12 +47,12 @@ describe Lederhosen::CLI do
57
47
  end
58
48
 
59
49
  it 'should split a fasta file into smaller fasta files (optionally gzipped)' do
60
- `./bin/lederhosen split_fasta --input=#{$test_dir}/trimmed/ILT_L_9_B_001.fasta --out-dir=#{$test_dir}/split/ --gzip true -n 100`
50
+ `./bin/lederhosen split_fasta --input=spec/data/trimmed/ILT_L_9_B_001.fasta --out-dir=#{$test_dir}/split/ --gzip true -n 100`
61
51
  $?.success?.should be_true
62
52
  end
63
53
 
64
54
  it 'should print representative sequences from uc files' do
65
- `./bin/lederhosen get_reps --input=#{$test_dir}/clusters.uc --database=#{$test_dir}/trimmed/ILT_L_9_B_001.fasta --output=#{$test_dir}/representatives.fasta`
55
+ `./bin/lederhosen get_reps --input=#{$test_dir}/clusters.uc --database=spec/data/trimmed/ILT_L_9_B_001.fasta --output=#{$test_dir}/representatives.fasta`
66
56
  end
67
57
 
68
58
  it 'should create a fasta file containing representative reads for each cluster'