cagnut_bwa 0.3.2 → 0.3.3

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
checksums.yaml CHANGED
@@ -1,7 +1,7 @@
1
1
  ---
2
2
  SHA1:
3
- metadata.gz: 97a5d23cc8316d1b59b094e87e31d8ee1aac2308
4
- data.tar.gz: 0b79d3b340ddb09a6c0d685e70d4f0a679cffecf
3
+ metadata.gz: 2e4b4f7800e674673c6c89535427a3338b1be1e4
4
+ data.tar.gz: b01d2b22d2987f7bd637338cba42cac186b04200
5
5
  SHA512:
6
- metadata.gz: ec0ef694d43675ba5e27548d63ab818e0fe57c1b574ffae58648a19da42fe84266fb2d0222ada1437c9f7b5ee59cac680e2ba4bde5c84cd4c96b4164cf28157e
7
- data.tar.gz: c36780081df3e223280c5ac07c917432df5b131f955cd5885d2c82a575850cafe7d2c0eb609b086158f18b5154f5a4d4df0a6b88f2c160998255e403e0071dda
6
+ metadata.gz: 4f4b424a1418f8a142723ba283d0fae36eee7e2f56354e50b7a3a6610982dd15ee4770304fe852a39686e0c95eb927ec8b0aed7d97031fff8c5e7f1f0be4d0d8
7
+ data.tar.gz: ed764a78b7f6e845908ec97bd39134c59ac37a2de02ac3d087fe8acb6f2ab6fa8ac93ee31b3ab287fa23bd90284335919cc622f548efb3418950dac7a4fa26b3
@@ -3,7 +3,7 @@ module CagnutBwa
3
3
  def check_tool tools_path, refs=nil
4
4
  super if defined?(super)
5
5
  ver = check_bwa tools_path['bwa'], refs['ref_fasta']
6
- check_bwa_index tools_path['bwa'], refs['ref_fasta'] if ver.blank?
6
+ check_bwa_index tools_path['bwa'], refs['ref_fasta'] if !ver.blank?
7
7
  end
8
8
 
9
9
  def check_bwa path, ref_path
@@ -26,14 +26,14 @@ else
26
26
  fi
27
27
 
28
28
  #check for correct number of sequences processed, based on fastq records
29
- PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o -P " \\d+ ")
29
+ PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o "[0-9]\+")
30
30
 
31
31
  echo "checking stdout file: " #{jobs_dir}/#{script_name}.err
32
32
  echo "bwa processed" $PROCESSED
33
33
 
34
34
  if [[ "#{input}" =~ gz$ ]]
35
35
  then
36
- LINESFASTQ1=$(zcat "#{input}" | wc -l)
36
+ LINESFASTQ1=$(gunzip -c "#{input}" | wc -l)
37
37
  else
38
38
  # non gz files
39
39
  LINESFASTQ1=$(wc -l "#{input}" | cut -d" " -f1 )
@@ -79,13 +79,13 @@ else
79
79
  fi
80
80
 
81
81
  #check for correct number of sequences processed, based on fastq records
82
- PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o -P " \\d+ ")
82
+ PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o "[0-9]\+")
83
83
  echo "checking stdout file: " #{jobs_dir}/#{script_name}.err
84
84
  echo "bwa processed" $PROCESSED
85
85
 
86
86
  if [[ "#{input2}" =~ gz$ ]]
87
87
  then
88
- LINESFASTQ2=$(zcat "#{input2}" | wc -l)
88
+ LINESFASTQ2=$(gunzip -c "#{input2}" | wc -l)
89
89
  else
90
90
  # non gz files
91
91
  LINESFASTQ2=$(wc -l "#{input2}" | cut -d" " -f1 )
@@ -26,7 +26,7 @@ else
26
26
  fi
27
27
 
28
28
  #check for correct number of sequences processed, based on fastq records
29
- PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o -P " \\d+ ")
29
+ PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o "[0-9]\+")
30
30
 
31
31
  echo "checking stdout file: " #{jobs_dir}/#{script_name}.err
32
32
  echo "bwa processed" $PROCESSED
@@ -34,8 +34,8 @@ echo "bwa processed" $PROCESSED
34
34
  if [[ "#{seq}" =~ gz$ ]]
35
35
  then
36
36
 
37
- LINESFASTQ1=$(zcat "#{seq}" | wc -l)
38
- LINESFASTQ2=$(zcat "#{seq2}" | wc -l)
37
+ LINESFASTQ1=$(gunzip -c "#{seq}" | wc -l)
38
+ LINESFASTQ2=$(gunzip -c "#{seq2}" | wc -l)
39
39
 
40
40
  else
41
41
  # non gz files
@@ -27,7 +27,7 @@ fi
27
27
 
28
28
  #check for correct number of sequences processed, based on fastq records
29
29
 
30
- PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o -P " \\d+ ")
30
+ PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o "[0-9]\+")
31
31
 
32
32
  echo "checking stdout file: " #{jobs_dir}/#{script_name}.err
33
33
  echo "bwa processed" $PROCESSED
@@ -35,15 +35,15 @@ else
35
35
  fi
36
36
 
37
37
  #check for correct number of sequences processed, based on fastq records
38
- PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o -P " \\d+ ")
38
+ PROCESSED=$(tail -5 #{jobs_dir}/#{script_name}.err | grep " sequences have been processed." | grep -o "[0-9]\+")
39
39
 
40
40
  echo "checking stdout file: " #{jobs_dir}/#{script_name}.err
41
41
  echo "bwa processed" $PROCESSED
42
42
 
43
43
  if [[ "#{fastq}" =~ gz$ ]]
44
44
  then
45
- LINESFASTQ1=$(zcat "#{fastq}" | wc -l)
46
- LINESFASTQ2=$(zcat "#{fastq2}.gz" | wc -l)
45
+ LINESFASTQ1=$(gunzip -c "#{fastq}" | wc -l)
46
+ LINESFASTQ2=$(gunzip -c "#{fastq2}" | wc -l)
47
47
  else
48
48
  # non gz files
49
49
  LINESFASTQ1=$(wc -l "#{fastq}" | cut -d" " -f1 )
@@ -1,3 +1,3 @@
1
1
  module CagnutBwa
2
- VERSION = "0.3.2"
2
+ VERSION = "0.3.3"
3
3
  end
metadata CHANGED
@@ -1,7 +1,7 @@
1
1
  --- !ruby/object:Gem::Specification
2
2
  name: cagnut_bwa
3
3
  version: !ruby/object:Gem::Version
4
- version: 0.3.2
4
+ version: 0.3.3
5
5
  platform: ruby
6
6
  authors:
7
7
  - Shi-Gang Wang
@@ -9,7 +9,7 @@ authors:
9
9
  autorequire:
10
10
  bindir: exe
11
11
  cert_chain: []
12
- date: 2016-11-07 00:00:00.000000000 Z
12
+ date: 2016-11-15 00:00:00.000000000 Z
13
13
  dependencies:
14
14
  - !ruby/object:Gem::Dependency
15
15
  name: cagnut_core