viral_seq 1.2.6 → 1.2.7

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: aee2536198f5579951f7c5f5f1b97ed27fb0bf5459e6dbfe79137a030951d443
-  data.tar.gz: 9ab9b69adf228d3429e02e5d47674ce0a0d137fa6231d940ef8d6cc9abfa1cdf
+  metadata.gz: 554845dba339d0e06b84c88bc117258516f391bdf58cce015c2669e7b2c6c0d5
+  data.tar.gz: 870280337c90d1f5b9ecbea6e6478d7e2dc22aa70917c6b2ecd94afaa185c1c6
 SHA512:
-  metadata.gz: d9b5def3e29f819b6d83902676e80fa2d81531cb44132942244125f6b63110d66389e7f9c61d6e9ea3d80ece57d22d95d5c3948c3f7f71d48835b2714794ffe0
-  data.tar.gz: 44b996d9caf9029f2ef522d9058410cdc3b4ec253ccc7bacca5d9ecf66f487624b0b6d18c960b6c7a0af7bfcd2bc72c8c02df59d3c4c7834ab80d4a9eed2a424
+  metadata.gz: 54db76e6fd8333ccebb19dee602378ec8dbe5d196ec7bd675e55f65db80cb06ac2ab51ce1f13ab7ea65c0a50ad49978bd3e9581074c497b298f0912858946fa8
+  data.tar.gz: 03d02329192465a9f278715c8a85e3a910e5c5c7252026980d29e669df823a5bdb4be323eeb56f7c9804b71fa8f1763c5a526227f3764315d6eb8e208934ce81

data/README.md

@@ -179,6 +179,12 @@ qc_seqhash.sdrm_hiv_pr(cut_off)
 
 ## Updates
 
+### Version 1.2.7-07152021
+
+  1. Optimzed the workflow of the `tcs` pipeline on raw data with uneven lengths. 
+  `tcs` version to v2.3.5.
+
+
 ### Version 1.2.6-07122021
 
   1. Optimized the workflow of the `tcs` pipeline in the "end-join/QC/Trimming" section.

data/bin/tcs

@@ -353,7 +353,8 @@ begin
       r1_seq_length = consensus_filtered.values[0][0].size
       r2_seq_length = consensus_filtered.values[0][1].size
     else
-      next
+      r1_seq_length = "n/a"
+      r2_seq_length = "n/a"
     end
     log.puts Time.now.to_s + "\t" + "R1 sequence #{r1_seq_length} bp"
     log.puts Time.now.to_s + "\t" + "R1 sequence #{r2_seq_length} bp"
@@ -426,7 +427,7 @@ begin
         if ref_end == 0
           ref_end = 0..(ViralSeq::RefSeq.get(ref_genome).size - 1)
         end
-        if primer[:end_join_option] == 1 and primer[:overlap] == 0
+        if primer[:end_join_option] == 1
           r1_sh = ViralSeq::SeqHash.fa(outfile_r1)
           r2_sh = ViralSeq::SeqHash.fa(outfile_r2)
           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

@@ -397,7 +397,9 @@ module ViralSeq
       (0..(seq_length - 1)).each do |position|
         all_base = []
         seq_array.each do |seq|
-          all_base << seq[position]
+          if seq[position]
+            all_base << seq[position]
+          end
         end
         base_count = all_base.count_freq
         max_base_list = []

data/lib/viral_seq/tcs_core.rb

@@ -305,7 +305,7 @@ module ViralSeq
       end
 
       def general_filter(seq)
-        if seq.size < $platform_sequencing_length
+        if seq.size < ($platform_sequencing_length - 1)
           return false
         elsif seq[1..-2] =~ /N/ # sequences with ambiguities except the 1st and last position removed
           return false

data/lib/viral_seq/version.rb

@@ -2,6 +2,6 @@
 # version info and histroy
 
 module ViralSeq
-  VERSION = "1.2.6"
-  TCS_VERSION = "2.3.5"
+  VERSION = "1.2.7"
+  TCS_VERSION = "2.3.6"
 end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: viral_seq
 version: !ruby/object:Gem::Version
-  version: 1.2.6
+  version: 1.2.7
 platform: ruby
 authors:
 - Shuntai Zhou
@@ -9,7 +9,7 @@ authors:
 autorequire:
 bindir: bin
 cert_chain: []
-date: 2021-07-12 00:00:00.000000000 Z
+date: 2021-07-15 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: bundler