cagnut_bwa 0.3.2 → 0.3.3

Sign up to get free protection for your applications and to get access to all the features.
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