resolve merge conflicts

README.md

@@ -1,9 +1,5 @@
 VILOCA: VIral LOcal haplotype reconstruction and mutation CAlling for short and long read data
 ===============
-[![Build Status](https://travis-ci.org/cbg-ethz/shorah.svg?branch=master)](https://travis-ci.org/cbg-ethz/shorah)
-[![Bioconda package](https://img.shields.io/conda/dn/bioconda/shorah.svg?label=Bioconda)](https://bioconda.github.io/recipes/shorah/README.html)
-[![Docker container](https://quay.io/repository/biocontainers/shorah/status)](https://quay.io/repository/biocontainers/shorah)
-
 
 VILOCA is an open source project for the analysis of next generation sequencing
 data. It is designed to analyse genetically heterogeneous samples. Its tools
@@ -17,22 +13,26 @@ genetic variants present in a mixed sample.
 For installation miniconda is recommended: https://docs.conda.io/en/latest/miniconda.html.
 We recommend to install VILOCA in a clean conda environment:
 ```
-conda create --name env_viloca libshorah
+conda create --name env_viloca --channel conda-forge --channel bioconda libshorah
 conda activate env_viloca
-pip install git+https://github.com/LaraFuhrmann/VILOCA@master
+pip install git+https://github.com/cbg-ethz/VILOCA@master
 ```
 
 ### Example
-To test your installation, we recommend running the program on `tests/data_1`.
+To test your installation run VILOCA `tests/data_1`:
+```
+viloca run -b test_aln.cram -f test_ref.fasta -z scheme.insert.bed --mode use_quality_scores
+```
+
 
 If the sequencing amplicon strategy is known, we recommend using the amplicon-mode of the program, which takes as input the `<smth>.insert.bed` - file:
-`shorah shotgun -b test_aln.cram -f test_ref.fasta -z scheme.insert.bed --mode use_quality_scores`
+`viloca run -b test_aln.cram -f test_ref.fasta -z scheme.insert.bed --mode use_quality_scores`
 
 If the sequencing quality scores are not trustable, the sequencing error parameters can also be learned:
-`shorah shotgun -b test_aln.cram -f test_ref.fasta -z scheme.insert.bed --mode learn_error_params`.
+`viloca run -b test_aln.cram -f test_ref.fasta -z scheme.insert.bed --mode learn_error_params`.
 
 If there is no information on the sequencing amplicon strategy available, run:
-`shorah shotgun -b test_aln.cram -f test_ref.fasta --mode use_quality_scores`
+`viloca run -b test_aln.cram -f test_ref.fasta --mode use_quality_scores`
 
 ### Parameters
 There are several parameters available:  

pyproject.toml

@@ -1,16 +1,16 @@
 [tool.poetry]
-name = "ShoRAH"
+name = "VILOCA"
 version = "0.1.0"
 description = "SHOrt Reads Assembly into Haplotypes"
 license = "GPL-3.0-only"
-authors = ["Benjamin Langer <[email protected]>"]
+authors = ["Benjamin Langer <[email protected]>, Lara Fuhrmann <[email protected]>"]
 build = "build.py"
 packages = [
-    { include = "shorah" }
+    { include = "viloca" }
 ]
 
 [tool.poetry.scripts]
-shorah = 'shorah.cli:main'
+viloca = 'viloca.cli:main'
 
 [tool.poetry.dependencies]
 python = ">=3.9.9,<3.11"

tests/data_1/shotgun_test.sh

@@ -1,4 +1,4 @@
 #!/bin/bash
 
-shorah shotgun -a 0.1 -w 201 -x 100000 -p 0.9 -c 0 \
+viloca run -a 0.1 -w 201 -x 100000 -p 0.9 -c 0 \
 -r HXB2:2469-3713 -R 42 -f test_ref.fasta -b test_aln.cram --out_format csv "$@"

tests/data_5/shotgun_prepare.sh

@@ -1,3 +1,3 @@
 #!/bin/bash
 
-shorah shotgun -a 0.1 -w 42 -x 100000 -p 0.9 -c 0 -r REF:43-273 -R 42 -b test_aln.cram -f ref.fasta
+viloca run -a 0.1 -w 42 -x 100000 -p 0.9 -c 0 -r REF:43-273 -R 42 -b test_aln.cram -f ref.fasta

tests/test_b2w.py

@@ -2,7 +2,7 @@
 import filecmp
 import os
 import glob
-from shorah import b2w, tiling
+from viloca import b2w, tiling
 import math
 import libshorah
 

tests/test_b2w_mapping.py

@@ -1,7 +1,7 @@
 from array import array
 import pytest
 from cigar import Cigar
-from shorah import b2w
+from viloca import b2w
 import hashlib
 
 class MockAlignedSegment:

tests/test_envp_post.py

@@ -1,5 +1,5 @@
 from unittest.mock import patch, mock_open
-from shorah import envp_post
+from viloca import envp_post
 
 DEFAULT_MOCK_DATA = "default mock data"
 

tests/test_pooled_post.py

@@ -1,5 +1,5 @@
 from unittest.mock import patch, mock_open
-from shorah import pooled_post
+from viloca import pooled_post
 import numpy as np
 
 DEFAULT_MOCK_DATA = "default mock data"
@@ -48,4 +48,4 @@ def open_side_effect(name):
 #                                      open("debug/w-HXB2-2938-3138.dbg"),
 #                                      open("support/w-HXB2-2938-3138.reads-support.fas"),
 #                                      open("corrected/w-HXB2-2938-3138.reads-cor.fas"),
-#                                      "shorah") # TODO
+#                                      "shorah") # TODO

tests/test_pooled_pre.py

@@ -1,6 +1,6 @@
 import pysam
 import os
-from shorah import pooled_pre
+from viloca import pooled_pre
 
 def test__annotate_alignment_file():
     out = "out.bam"
@@ -48,4 +48,4 @@ def test_pre_process_pooled():
     os.remove(out + ".bai")
 
     assert a[0] == a[1] != 0
-    assert a[2] == 0
+    assert a[2] == 0

tests/test_shorah_snv.py

@@ -1,5 +1,5 @@
 import pytest
-from shorah.shorah_snv import _compare_ref_to_read, SNP_id, SNV
+from viloca.shorah_snv import _compare_ref_to_read, SNP_id, SNV
 
 
 @pytest.mark.parametrize("ref, seq, spec", [
@@ -39,4 +39,4 @@ def test_compare_ref_to_read(ref, seq, spec):
 
     assert snp == spec
 
-    assert tot_snv == len(snp)
+    assert tot_snv == len(snp)

tests/test_tiling.py

@@ -1,4 +1,4 @@
-from shorah import tiling
+from viloca import tiling
 import pytest
 
 def test_equispaced():

shorah/__main__.py → viloca/__main__.py

@@ -8,7 +8,7 @@
 - https://docs.python.org/2/using/cmdline.html#cmdoption-m
 - https://docs.python.org/3/using/cmdline.html#cmdoption-m
 """
-from shorah.cli import main
+from viloca.cli import main
 
 if __name__ == "__main__":
     main()