viral_seq 1.2.9 → 1.6.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 2ccdd3f318deb458540b703ae00afcbdd195a92145c22dfb6c95ad14ac366146
-  data.tar.gz: 4f131cd69fef8548c1a466c2e23be008b08bed50ca78ef12195e4d4b39e052eb
+  metadata.gz: 6bc062e2f207fd9e33d97885125b656c4b24ae8056221020bb611c191104109b
+  data.tar.gz: a0086f02da866821c134a9433317b7016a429e37e775de65d434dbae1b7efd23
 SHA512:
-  metadata.gz: 4cb44299c0061794b8c32c23c850cd9c53eb9245d7124ed75d35367225e600e2c630df9853d13a8003a8a1c6c793edea58a188871f328866043132585846e825
-  data.tar.gz: 92594266725b7691c5203cd478e84a0fea7fc6d08124efad6311a4e2bdcc7dbd0d6218a54e5b63d19e3c3a8ad789c49906388bac7ff59611bc5810e1c614d8ef
+  metadata.gz: dc29bb6914196e4aa0fa6d7c2db25d4db1179f40a551dc0784534fe2225c857ef22663815d483412db1a715977574a01ab457a5d4913e38a7a9ac439b38795ab
+  data.tar.gz: f42559bc07b75b2d9a4023af0313cc40e037697970c0ac4a6c34fd04938b5ab879ece19fd95d05237075292c38a93e315018c7e0c03bd0ee9ab8cf1716b894ed

data/Gemfile.lock

@@ -1,12 +1,12 @@
 PATH
   remote: .
   specs:
-    viral_seq (1.1.1)
-      colorize (>= 0.1)
-      combine_pdf (>= 1.0.0)
-      muscle_bio (>= 0.4)
-      prawn (>= 2.3.0)
-      prawn-table (>= 0.2.0)
+    viral_seq (1.3.0)
+      colorize (~> 0.1)
+      combine_pdf (~> 1.0, >= 1.0.0)
+      muscle_bio (~> 0.4)
+      prawn (~> 2.3, >= 2.3.0)
+      prawn-table (~> 0.2, >= 0.2.0)
 
 GEM
   remote: https://rubygems.org/

data/README.md

@@ -55,7 +55,7 @@ Example commands:
 ---
 ### `tcs_log`
 
-Use `tcs_log` script to pool run logs and TCS fasta files after one batch of `tcs` jobs.
+Use `tcs_log` script to pool run logs and TCS fasta files after one batch of `tcs` jobs. This command generates log.html to visualize the sequencing runs.
 
 
 Example file structure:  
@@ -179,10 +179,32 @@ qc_seqhash.sdrm_hiv_pr(cut_off)
 
 ## Updates
 
+### Version 1.6.0-01042022
+
+  1. Update the `ViralSeq::TcsCore::detection_limit` with pre-calculated values to save processing time.
+  2. Update `tcs` pipeline to v2.5.0. HTML report will generated after running `tcs_log` script after `tcs` pipeline.
+
+### Version 1.5.0-01042022
+
+  1. Added a function to calcute detection limit/sensitivity for minority variants (R required). `ViralSeq::TcsCore::detection_limit`
+  2. Added a function to get a sub SeqHash object given a range of nt positions. `ViralSeq::SeqHash#nt_range`
+  3. Added a function to quality check dna sequences comparing with sample consensus for indels. `ViralSeq::SeqHash#qc_indel`
+  4. Added a function for DNA variant analysis. Return a Hash object that can output as a JSON file. `ViralSeq::SeqHash#nt_variants`
+  5. Added a function to check the size of sequences of a SeqHash object. `ViralSeq::SeqHash#check_nt_size`
+
+### Version 1.4.0-10132021
+
+  1. Added a function to calculate false detectionr rate (FDR, aka, Benjamini-Hochberg correction) for minority mutations detected in the sequences. `ViralSeq::SeqHash#fdr`
+  2. Updated `bin\tcs_sdrm` script to add FDR value to each DRMs detected.
+
+### Version 1.3.0-08302021
+
+  1. Fixed a bug in the `tcs` pipeline.
+
 ### Version 1.2.9-08022021
 
   1. Fixed a bug when reading the input primer sequences in lowercases.
-  2. Fix a bug in the method ViralSeq::Math::RandomGaussian
+  2. Fixed a bug in the method ViralSeq::Math::RandomGaussian
 
 ### Version 1.2.8-07292021
 

data/bin/tcs

@@ -32,9 +32,18 @@ require 'optparse'
 
 options = {}
 
-banner = '-'*50 + "\n" +
-        '| The TCS Pipeline ' + "Version #{ViralSeq::TCS_VERSION}".red.bold + " by " + "Shuntai Zhou".blue.bold + ' |' + "\n" +
-        '-'*50 + "\n"
+# banner = '-'*50 + "\n" +
+#         '| The TCS Pipeline ' + "Version #{ViralSeq::TCS_VERSION}".red.bold + " by " + "Shuntai Zhou".blue.bold + ' |' + "\n" +
+#         '-'*50 + "\n"
+
+banner = "\n" +
+"████████  ██████ ███████     ██████  ██ ██████  ███████ ██      ██ ███    ██ ███████\n".light_red +
+"   ██    ██      ██          ██   ██ ██ ██   ██ ██      ██      ██ ████   ██ ██\n".light_yellow +
+"   ██    ██      ███████     ██████  ██ ██████  █████   ██      ██ ██ ██  ██ █████\n".light_green +
+"   ██    ██           ██     ██      ██ ██      ██      ██      ██ ██  ██ ██ ██\n".light_cyan +
+"   ██     ██████ ███████     ██      ██ ██      ███████ ███████ ██ ██   ████ ███████\n".light_magenta
+
+banner += "\nVersion #{ViralSeq::TCS_VERSION}".red.bold + " by " + "Shuntai Zhou".blue.bold + "\n\n"
 
 OptionParser.new do |opts|
   opts.banner = banner + "Usage: tcs -j"
@@ -200,7 +209,7 @@ begin
     summary_json[:paired_raw_sequence] = paired_seq_number
     if paired_seq_number < raw_sequence_number * 0.001
       summary_json[:warnings] <<
-        "WARNING: Filtered raw sequneces less than 0.1% of the total raw sequences. Possible contamination."
+        "WARNING: Filtered raw sequences less than 0.1% of the total raw sequences. Possible contamination."
     end
 
     common_keys.each do |seqtag|
@@ -401,7 +410,11 @@ begin
       when 4
         joined_sh = shp.join2(model: :indiv)
       end
-      return joined_sh
+      if joined_sh
+        return joined_sh
+      else
+        joined_sh = ViralSeq::SeqHash.new
+      end
     end
 
     if primer[:end_join]