fastq-factory 0.1.6 → 0.1.7
Sign up to get free protection for your applications and to get access to all the features.
- data/VERSION +1 -1
- data/fastq-factory.gemspec +2 -2
- data/lib/trim_and_correct.rb +5 -3
- metadata +3 -3
data/VERSION
CHANGED
@@ -1 +1 @@
|
|
1
|
-
0.1.
|
1
|
+
0.1.7
|
data/fastq-factory.gemspec
CHANGED
@@ -5,11 +5,11 @@
|
|
5
5
|
|
6
6
|
Gem::Specification.new do |s|
|
7
7
|
s.name = "fastq-factory"
|
8
|
-
s.version = "0.1.
|
8
|
+
s.version = "0.1.7"
|
9
9
|
|
10
10
|
s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
|
11
11
|
s.authors = ["Anthony Underwood"]
|
12
|
-
s.date = "2012-
|
12
|
+
s.date = "2012-09-07"
|
13
13
|
s.description = "This tool can process fastq files, using fastq_quality_trimmer and quake to correct fastq files and then provide a quality asssessment of the data"
|
14
14
|
s.email = "anthony.underwood@hpa.org.uk"
|
15
15
|
s.executables = ["fastq-factory"]
|
data/lib/trim_and_correct.rb
CHANGED
@@ -2,6 +2,7 @@ def trim_and_correct_fastqs(sample_map, directory, forward_reads_suffix, forward
|
|
2
2
|
Dir.chdir(directory)
|
3
3
|
# trimming
|
4
4
|
sample_map.each do |sample_file_prefix, sample_name|
|
5
|
+
next if File.exists?("paired_#{sample_file_prefix}#{forward_reads_suffix}.trimmed.#{forward_reads_file_extension}")
|
5
6
|
puts "Trimming files for #{sample_name}"
|
6
7
|
#determine read length
|
7
8
|
read_length = calculate_read_length("#{directory}/#{sample_file_prefix}#{forward_reads_suffix}.#{forward_reads_file_extension}")
|
@@ -15,6 +16,7 @@ def trim_and_correct_fastqs(sample_map, directory, forward_reads_suffix, forward
|
|
15
16
|
# quake correction
|
16
17
|
# write file for quake
|
17
18
|
sample_map.each do |sample_file_prefix, sample_name|
|
19
|
+
next if File.exists?("#{sample_file_prefix}#{forward_reads_suffix}.trimmed.cor.#{forward_reads_file_extension}") || File.exists?(("paired_#{sample_file_prefix}#{reverse_reads_suffix}.trimmed.cor.#{forward_reads_file_extension}"))
|
18
20
|
puts "Error correcting files for #{sample_name}"
|
19
21
|
output_file = File.open("quake_file_list.txt","w")
|
20
22
|
output_file.puts "paired_#{sample_file_prefix}#{forward_reads_suffix}.trimmed.#{forward_reads_file_extension} paired_#{sample_file_prefix}#{reverse_reads_suffix}.trimmed.#{reverse_reads_file_extension}"
|
@@ -23,9 +25,9 @@ def trim_and_correct_fastqs(sample_map, directory, forward_reads_suffix, forward
|
|
23
25
|
`#{quake_path} -f quake_file_list.txt -k 15 -q #{quality_scale}`
|
24
26
|
end
|
25
27
|
sample_map.each do |sample_file_prefix, sample_name|
|
26
|
-
if File.exists?("
|
27
|
-
|
28
|
-
|
28
|
+
next if File.exists?("#{sample_file_prefix}#{forward_reads_suffix}.trimmed.cor.#{forward_reads_file_extension}")
|
29
|
+
# remove orphans
|
30
|
+
`perl /tmp/fastq-remove-orphans.pl -1 paired_#{sample_file_prefix}#{forward_reads_suffix}.trimmed.cor.#{forward_reads_file_extension} -2 paired_#{sample_file_prefix}#{reverse_reads_suffix}.trimmed.cor.#{reverse_reads_file_extension}`
|
29
31
|
end
|
30
32
|
|
31
33
|
# cleanup and rename files
|
metadata
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: fastq-factory
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
version: 0.1.
|
4
|
+
version: 0.1.7
|
5
5
|
prerelease:
|
6
6
|
platform: ruby
|
7
7
|
authors:
|
@@ -9,7 +9,7 @@ authors:
|
|
9
9
|
autorequire:
|
10
10
|
bindir: bin
|
11
11
|
cert_chain: []
|
12
|
-
date: 2012-
|
12
|
+
date: 2012-09-07 00:00:00.000000000 Z
|
13
13
|
dependencies:
|
14
14
|
- !ruby/object:Gem::Dependency
|
15
15
|
name: trollop
|
@@ -165,7 +165,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
|
|
165
165
|
version: '0'
|
166
166
|
segments:
|
167
167
|
- 0
|
168
|
-
hash: -
|
168
|
+
hash: -4361859413612648693
|
169
169
|
required_rubygems_version: !ruby/object:Gem::Requirement
|
170
170
|
none: false
|
171
171
|
requirements:
|