All notable changes to Quasitools will be documented in this file.
Commands:
- complexity bam: measures the complexity of a quasispecies across a reference genome using aligned reads in BAM format
Options:
- complexity: now has the option to filter out low frequency haplotypes
complexity
has been replaced withcomplexity fasta
andcomplexity bam
- Fixed a complexity bug with haplotype counting and consensus sequence generation.
- Fixed a complexity bug when parsing haplotypes with indels.
- Integrated PyAAVF and changed output of mutation report in
hydra
command to produce AAVF file along withcall aavar
to produced AAVF file instead of HMCF file. - Made various changes to the documentation to improve clarity.
- Fixed quasitools distance sometimes crashing because of uneven coverage.
- Improved the execution time of quasitools distance.
- Fixed the Hamming distance calculation using integers instead of floats, depending on the Python version.
- Use MkDocs for documentation
- Make VCF dependency optional in aavar and add optional VCF dependency in codonvar
- Fix hydra consensus output: exclude header line if there are no sequence lines
- Remove a deprecated regular expression sequence in reference_parser.py
2018-09-25
Commands:
- complexity: measures the complexity of a quasispecies
- Improved read parsing, pileup generation, and related interal data structures.
2018-07-25
- Fix bug in quality_control where iteratively trimmed reads are not used when printing output file.
2018-07-24
- Updated hydra command and quality command to fix boolean flag for enabling median or mean score cutoff value.
2018-06-29
- commands:
- quality: performs quality control on FASTQ reads
- modules:
- quality_control.py: contains QualityControl class
- Updated hydra command to use new quality control class, masking and/or trimming reads based on command-line options that user has specified.
- Fixed bug in cmd_consensus.py which occured when there are multiple records in the reference file and the user is writing to a file.
2018-04-13
- Updated hydra command to accept an id to be used as the sequence identifier in the consensus report, and as the RG-id in the bam file
2018-04-10
- commands:
- distance: returns the evolutionary distances between viral quasispecies as a distance matrix
- modules:
- pileup.py: contains Pileup and Pileup_List classes
- distance.py: contains Distance_Matrix class
2018-03-16
- Fixed FASTA identifier in consensus output for the consensus and hydra commands
- User can pass in a default id as an option
- RG tag is used if present in the supplied .bam file
2017-11-27
- Updated samtools dependency to v1.3
2017-11-22
- fixed call codonvar command to properly output to either file or stdout depending on flag
2017-11-03
- commands:
- hydra: which identifies HIV Drug Resistance in a NGS sample
- call: new base command for nt, aa, and codon variant calling
- ntvar: formally just callntvar
- aavar: formally aavariants
- codonvar: reports codon variants for use with dnds command
- dnds: calculates dnds values per coding region using codonvar output
- Renamed aa_census command to aa_coverage
- Hardended command options and parameters
2017-09-06
This is the initial release of Quasitools