Project for centrally supported wrapper of Mathijs Sanders' hairpin flagging algorithms, known as "Mathijs' scripts"
CLI for statistical detection and flagging of variants caused by hairpin/cruciform artifacts in LCM sequencing. Previously, this task was being performed with diffuse versions of a small pipeline built by Mathijs and added to by others over many years. This repository brings together the components of that pipeline, strips out the extraneous functionality, and packages them into a single command for ease of use. The components are inherited and remain in their original form - these are runScriptImitateANNOVAR.sh and additionalBamStatistics.jar. The parent repositories can be found here and here respectively. There is an associated paper here; information on the calculated statistics can be found in the 'SNV filtering' section.
The various components that have been removed or pared down either produced statistics that are not utilised for hairpin detection, or were program functionality which would be better placed elsewhere, such as prefiltering of input VCFs according to CPLM.
Java >= 8
Python >= 3.10
samtools == 1.14
pysam == 0.19.1
vcfpy == 0.13.4
clone repository and cd into bin/ and run the following to install into a virtual environment:
python3.10 -m venv pyenv
source pyenv/bin/activate
pip install -r requirements.txt
deactivate
Usage:
hairpin \
[ Mandatory ]
-v input VCF \
-b BAM file corresponding to VCF \
-g path to reference genome fasta \
[ Optional ]
-o output directory (defaults to current working dir) \
-m set java heap memory (default 10G) \
-h display usage \
You must have the .bai, .bas, and .met.gz files associated with the .bam file in the same directory as the .bam file specified with -b.
The tool will output the input VCF file updated with annotations indicating potential hairpin artificats.
"Script has failed" If the core component of the hairpin process fails it will usually output this cryptic error. Unfortunately, this is caused by external code and we cannot make this error more informative at this time. It is most likely caused by running the script without enough memory available.
Copyright (c) 2023 Genome Research Ltd.
Author: CASM/Cancer IT <cgphelp@sanger.ac.uk>
This file is part of hairpin-wrapper.
hairpin-wrapper is free software: you can redistribute it and/or modify it under
the terms of the GNU Affero General Public License as published by the Free
Software Foundation; either version 3 of the License, or (at your option) any
later version.
This program is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
FOR A PARTICULAR PURPOSE. See the GNU Affero General Public License for more
details.
You should have received a copy of the GNU Affero General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>.