ViralSeq
A Ruby Gem containing bioinformatics tools for processing viral NGS data.
Specifically for Primer-ID sequencing and HIV drug resistance analysis.
Installation
$ gem install viral_seq
Usage
Load all ViralSeq classes by requiring 'viral_seq.rb'
#!/usr/bin/env ruby
require 'viral_seq'
Use executable locator
to get the coordinates of the sequences on HIV/SIV reference genome from a FASTA file through a terminal
$ locator -i sequence.fasta -o sequence.fasta.csv
Use executable tcs
pipeline to process Primer ID MiSeq sequencing data. Parameter json file can be generated using tcs_json_generator
or at https://tcs-dr-dept-tcs.cloudapps.unc.edu/generator.php
$ tcs params.json
Use executable tcs_json_generator
to generate params .json file for the tcs
pipeline.
$ tcs_json_generator
Some Examples
Load nucleotide sequences from a FASTA format sequence file
my_seqhash = ViralSeq::SeqHash.fa('my_seq_file.fasta')
Make an alignment (using MUSCLE)
aligned_seqhash = my_seqhash.align
Filter nucleotide sequences with the reference coordinates (HIV Protease)
qc_seqhash = aligned_seqhash.hiv_seq_qc(2253, 2549, false, :HXB2)
Further filter out sequences with Apobec3g/f hypermutations
qc_seqhash = qc_seqhash.a3g
Calculate nucleotide diveristy π
qc_seqhash.pi
Calculate cut-off for minority variants based on Poisson model
cut_off = qc_seqhash.pm
Examine for drug resistance mutations for HIV PR region
qc_seqhash.sdrm_hiv_pr(cut_off)
Updates
Version 1.0.9-07182020:
Change ViralSeq::SeqHash#stop_codon and ViralSeq::SeqHash#a3g_hypermut return value to hash object.
TCS pipeline updated to version 2.0.1. Add optional
export_raw: TRUE/FALSE
in json params. Ifexport_raw
isTRUE
, raw sequence reads (have to pass quality filters) will be exported, along with TCS reads.
Version 1.0.8-02282020:
TCS pipeline (version 2.0.0) added as executable. tcs - main TCS pipeline script. tcs_json_generator - step-by-step script to generate json file for tcs pipeline.
Methods added: ViralSeq::SeqHash#trim
Bug fix for several methods.
Version 1.0.7-01282020:
- Several methods added, including ViralSeq::SeqHash#error_table ViralSeq::SeqHash#random_select
- Improved performance for several functions.
Version 1.0.6-07232019:
- Several methods added to ViralSeq::SeqHash, including ViralSeq::SeqHash#size ViralSeq::SeqHash#+ ViralSeq::SeqHash#write_nt_fa ViralSeq::SeqHash#mutation
- Update documentations and rspec samples.
Version 1.0.5-07112019:
- Update ViralSeq::SeqHash#sequence_locator.
Program will try to determine the direction (
+
or-
of the query sequence) - update executable
locator
to have a column ofdirection
in output .csv file
Version 1.0.4-07102019:
- Use home directory (Dir.home) instead of the directory of the script file for temp MUSCLE file.
- Fix bugs in bin
locator
Version 1.0.3-07102019:
- Bug fix.
Version 1.0.2-07102019:
- Fixed a gem loading issue.
Version 1.0.1-07102019:
- Add keyword argument :model to ViralSeq::SeqHashPair#join2.
- Add method ViralSeq::SeqHash#sequence_locator (also: #loc), a function to locate sequences on HIV/SIV reference genomes, as HIV Sequence Locator from LANL.
- Add executable 'locator'. An HIV/SIV sequence locator tool similar to LANL Sequence Locator.
- update documentations
Version 1.0.0-07092019:
- Rewrote the whole ViralSeq gem, grouping methods into modules and classes under main Module::ViralSeq
Development
Bug reports and pull requests are welcome on GitHub at https://github.com/ViralSeq/viral_seq. This project is intended to be a safe, welcoming space for collaboration, and contributors are expected to adhere to the Contributor Covenant code of conduct.
License
The gem is available as open source under the terms of the MIT License.
Code of Conduct
Everyone interacting in the viral_seq project’s codebases, issue trackers, chat rooms and mailing lists is expected to follow the code of conduct.