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 +4 -4
- data/lib/cagnut_bwa/check_tools.rb +1 -1
- data/lib/cagnut_bwa/functions/templates/aln.sh +4 -4
- data/lib/cagnut_bwa/functions/templates/bwa_samp_one_fastq.sh +3 -3
- data/lib/cagnut_bwa/functions/templates/bwa_samse_one_fastq.sh +1 -1
- data/lib/cagnut_bwa/functions/templates/samp.sh +3 -3
- data/lib/cagnut_bwa/version.rb +1 -1
- metadata +2 -2
checksums.yaml
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
---
|
2
2
|
SHA1:
|
3
|
-
metadata.gz:
|
4
|
-
data.tar.gz:
|
3
|
+
metadata.gz: 2e4b4f7800e674673c6c89535427a3338b1be1e4
|
4
|
+
data.tar.gz: b01d2b22d2987f7bd637338cba42cac186b04200
|
5
5
|
SHA512:
|
6
|
-
metadata.gz:
|
7
|
-
data.tar.gz:
|
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 -
|
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=$(
|
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 -
|
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=$(
|
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 -
|
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=$(
|
38
|
-
LINESFASTQ2=$(
|
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 -
|
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 -
|
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=$(
|
46
|
-
LINESFASTQ2=$(
|
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 )
|
data/lib/cagnut_bwa/version.rb
CHANGED
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.
|
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-
|
12
|
+
date: 2016-11-15 00:00:00.000000000 Z
|
13
13
|
dependencies:
|
14
14
|
- !ruby/object:Gem::Dependency
|
15
15
|
name: cagnut_core
|