viral_seq 1.2.6 → 1.2.7
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/README.md +6 -0
- data/bin/tcs +3 -2
- data/lib/viral_seq/seq_hash.rb +3 -1
- data/lib/viral_seq/tcs_core.rb +1 -1
- data/lib/viral_seq/version.rb +2 -2
- metadata +2 -2
checksums.yaml
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
---
|
2
2
|
SHA256:
|
3
|
-
metadata.gz:
|
4
|
-
data.tar.gz:
|
3
|
+
metadata.gz: 554845dba339d0e06b84c88bc117258516f391bdf58cce015c2669e7b2c6c0d5
|
4
|
+
data.tar.gz: 870280337c90d1f5b9ecbea6e6478d7e2dc22aa70917c6b2ecd94afaa185c1c6
|
5
5
|
SHA512:
|
6
|
-
metadata.gz:
|
7
|
-
data.tar.gz:
|
6
|
+
metadata.gz: 54db76e6fd8333ccebb19dee602378ec8dbe5d196ec7bd675e55f65db80cb06ac2ab51ce1f13ab7ea65c0a50ad49978bd3e9581074c497b298f0912858946fa8
|
7
|
+
data.tar.gz: 03d02329192465a9f278715c8a85e3a910e5c5c7252026980d29e669df823a5bdb4be323eeb56f7c9804b71fa8f1763c5a526227f3764315d6eb8e208934ce81
|
data/README.md
CHANGED
@@ -179,6 +179,12 @@ qc_seqhash.sdrm_hiv_pr(cut_off)
|
|
179
179
|
|
180
180
|
## Updates
|
181
181
|
|
182
|
+
### Version 1.2.7-07152021
|
183
|
+
|
184
|
+
1. Optimzed the workflow of the `tcs` pipeline on raw data with uneven lengths.
|
185
|
+
`tcs` version to v2.3.5.
|
186
|
+
|
187
|
+
|
182
188
|
### Version 1.2.6-07122021
|
183
189
|
|
184
190
|
1. Optimized the workflow of the `tcs` pipeline in the "end-join/QC/Trimming" section.
|
data/bin/tcs
CHANGED
@@ -353,7 +353,8 @@ begin
|
|
353
353
|
r1_seq_length = consensus_filtered.values[0][0].size
|
354
354
|
r2_seq_length = consensus_filtered.values[0][1].size
|
355
355
|
else
|
356
|
-
|
356
|
+
r1_seq_length = "n/a"
|
357
|
+
r2_seq_length = "n/a"
|
357
358
|
end
|
358
359
|
log.puts Time.now.to_s + "\t" + "R1 sequence #{r1_seq_length} bp"
|
359
360
|
log.puts Time.now.to_s + "\t" + "R1 sequence #{r2_seq_length} bp"
|
@@ -426,7 +427,7 @@ begin
|
|
426
427
|
if ref_end == 0
|
427
428
|
ref_end = 0..(ViralSeq::RefSeq.get(ref_genome).size - 1)
|
428
429
|
end
|
429
|
-
if primer[:end_join_option] == 1
|
430
|
+
if primer[:end_join_option] == 1
|
430
431
|
r1_sh = ViralSeq::SeqHash.fa(outfile_r1)
|
431
432
|
r2_sh = ViralSeq::SeqHash.fa(outfile_r2)
|
432
433
|
r1_sh = r1_sh.hiv_seq_qc(ref_start, (0..(ViralSeq::RefSeq.get(ref_genome).size - 1)), indel, ref_genome)
|
data/lib/viral_seq/seq_hash.rb
CHANGED
@@ -397,7 +397,9 @@ module ViralSeq
|
|
397
397
|
(0..(seq_length - 1)).each do |position|
|
398
398
|
all_base = []
|
399
399
|
seq_array.each do |seq|
|
400
|
-
|
400
|
+
if seq[position]
|
401
|
+
all_base << seq[position]
|
402
|
+
end
|
401
403
|
end
|
402
404
|
base_count = all_base.count_freq
|
403
405
|
max_base_list = []
|
data/lib/viral_seq/tcs_core.rb
CHANGED
@@ -305,7 +305,7 @@ module ViralSeq
|
|
305
305
|
end
|
306
306
|
|
307
307
|
def general_filter(seq)
|
308
|
-
if seq.size < $platform_sequencing_length
|
308
|
+
if seq.size < ($platform_sequencing_length - 1)
|
309
309
|
return false
|
310
310
|
elsif seq[1..-2] =~ /N/ # sequences with ambiguities except the 1st and last position removed
|
311
311
|
return false
|
data/lib/viral_seq/version.rb
CHANGED
metadata
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: viral_seq
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
version: 1.2.
|
4
|
+
version: 1.2.7
|
5
5
|
platform: ruby
|
6
6
|
authors:
|
7
7
|
- Shuntai Zhou
|
@@ -9,7 +9,7 @@ authors:
|
|
9
9
|
autorequire:
|
10
10
|
bindir: bin
|
11
11
|
cert_chain: []
|
12
|
-
date: 2021-07-
|
12
|
+
date: 2021-07-15 00:00:00.000000000 Z
|
13
13
|
dependencies:
|
14
14
|
- !ruby/object:Gem::Dependency
|
15
15
|
name: bundler
|