varvamp 1.2.1__py3-none-any.whl → 1.2.2__py3-none-any.whl

varvamp/__init__.py

@@ -1,3 +1,3 @@
 """Tool to design amplicons for highly variable virusgenomes"""
 _program = "varvamp"
-__version__ = "1.2.1"
+__version__ = "1.2.2"

varvamp/command.py

@@ -119,7 +119,7 @@ def get_args(sysargs):
         type=int,
         metavar="100",
         default=100,
-        help="min overlap of the amplicons"
+        help="min overlap of the amplicon inserts"
     )
     SINGLE_parser.add_argument(
         "-n",
@@ -505,8 +505,8 @@ def main():
     # write files that are shared in all modes
     reporting.write_regions_to_bed(primer_regions, args.name, data_dir)
     reporting.write_alignment(data_dir, alignment_cleaned)
-    reporting.write_fasta(data_dir, f"majority_consensus", f"{args.name}_consensus",majority_consensus)
-    reporting.write_fasta(results_dir, f"ambiguous_consensus", f"{args.name}_consensus", ambiguous_consensus)
+    reporting.write_fasta(data_dir, f"majority_consensus", f"{args.name}_majority_consensus",majority_consensus)
+    reporting.write_fasta(results_dir, f"ambiguous_consensus", f"{args.name}_ambiguous_consensus", ambiguous_consensus)
 
     # Functions called from here on return lists of amplicons that are refined step-wise into final schemes.
     # These lists that are passed between functions and later used for reporting consist of dictionary elemnts,

varvamp/scripts/logging.py

@@ -147,7 +147,7 @@ def raise_arg_errors(args, log_file):
                 log_file,
                 exit=True
             )
-        if args.overlap < 50:
+        if args.overlap < 10:
             raise_error(
                 "small overlaps might hinder downstream analyses. Consider increasing.",
                 log_file

varvamp/scripts/reporting.py

@@ -53,7 +53,7 @@ def write_regions_to_bed(primer_regions, scheme_name, path, mode=None):
     with open(outfile, 'w') as o:
         for counter, region in enumerate(primer_regions):
             print(
-                f"{scheme_name}_consensus",
+                f"{scheme_name}_ambiguous_consensus",
                 region[0],
                 region[1],
                 "REGION_"+str(counter),
@@ -68,9 +68,7 @@ def write_primers_to_bed(outfile, scheme_name, primer_name, primer_properties, n
     """
     with open(outfile, 'a') as o:
         # write header for primer bed
-        if os.path.getsize(outfile) == 0 and sequence is not None:
-            print("#chrom\tchromStart\tchromEnd\tprimer-name\tpool\tstrand\tprimer-sequence", file=o)
-        data = [f"{scheme_name}_consensus",
+        data = [f"{scheme_name}_ambiguous_consensus",
             primer_properties[1],  # start
             primer_properties[2],  # stop
             primer_name,
@@ -150,7 +148,7 @@ def write_qpcr_to_files(path, final_schemes, ambiguous_consensus, scheme_name, l
             amp_name = f"{scheme_name}_{n}"
             # write bed amplicon file
             print(
-                f"{scheme_name}_consensus",
+                f"{scheme_name}_ambiguous_consensus",
                 amp["LEFT"][1],
                 amp["RIGHT"][2],
                 amp_name,
@@ -329,7 +327,7 @@ def write_scheme_to_files(path, amplicon_scheme, ambiguous_consensus, scheme_nam
         # write amplicon bed with amplicons sorted by start position
         for record in sorted(amplicon_bed_records, key=lambda x: x[0]):
             print(
-                f"{scheme_name}_consensus",
+                f"{scheme_name}_ambiguous_consensus",
                 *record,
                 ".",
                 sep="\t",

varvamp/scripts/scheme.py

@@ -102,14 +102,14 @@ def create_amplicon_graph(amplicons, min_overlap):
 
     # add the maximum len of a primer to ensure that possible amplicon starts
     # before the min overlap
-    min_overlap = min_overlap + config.PRIMER_SIZES[2]
+    min_overlap = min_overlap + config.PRIMER_SIZES[1]
 
     for current_amplicon in amplicons:
         # remember all vertices
         amplicon_id = current_amplicon["id"]
         nodes.append(amplicon_id)
         start = current_amplicon["LEFT"][1] + current_amplicon["length"]/2
-        stop = current_amplicon["RIGHT"][2] - min_overlap
+        stop = current_amplicon["RIGHT"][1] - min_overlap
         for next_amplicon in amplicons:
             # check if the next amplicon lies within the start/stop range of
             # the current amplicon and if its non-overlapping part is large

{varvamp-1.2.1.dist-info → varvamp-1.2.2.dist-info}/METADATA

@@ -1,6 +1,6 @@
-Metadata-Version: 2.1
+Metadata-Version: 2.2
 Name: varvamp
-Version: 1.2.1
+Version: 1.2.2
 Summary: Variable VirusAMPlicons (varVAMP) is a tool to design primers for highly diverse viruses
 Home-page: https://github.com/jonas-fuchs/varVAMP
 Author: Dr. Jonas Fuchs
@@ -9,12 +9,21 @@ Classifier: Programming Language :: Python :: 3.9
 Classifier: License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)
 Requires-Python: >=3.9
 Description-Content-Type: text/markdown
-Requires-Dist: biopython >=1.79
-Requires-Dist: matplotlib >=3.5.1
-Requires-Dist: primer3-py >=1.1.0
-Requires-Dist: pandas >=1.4.4
-Requires-Dist: numpy >=1.23.3
-Requires-Dist: seqfold >=0.7.15
+Requires-Dist: biopython>=1.79
+Requires-Dist: matplotlib>=3.5.1
+Requires-Dist: primer3-py>=1.1.0
+Requires-Dist: pandas>=1.4.4
+Requires-Dist: numpy>=1.23.3
+Requires-Dist: seqfold>=0.7.15
+Dynamic: author
+Dynamic: author-email
+Dynamic: classifier
+Dynamic: description
+Dynamic: description-content-type
+Dynamic: home-page
+Dynamic: requires-dist
+Dynamic: requires-python
+Dynamic: summary
 
 
 **var**iable **V**irus**AMP**licons (varVAMP) is a tool to design primers for highly diverse viruses. The input is an alignment of your viral (full-genome) sequences.

{varvamp-1.2.1.dist-info → varvamp-1.2.2.dist-info}/RECORD

@@ -1,21 +1,21 @@
-varvamp/__init__.py,sha256=dcU1t4AmPwg3LoG01RN2Q4VvzPKoZ00gRfg-qoofv28,107
+varvamp/__init__.py,sha256=wCHzWXwkVwFmeC5JBU1NcQH5A-wYZkWuit357pqziQA,107
 varvamp/__main__.py,sha256=9R3mbX2_Q5LByPx8WLoTbvZ-G2dbRSMpDlgze0Wm5Fc,98
-varvamp/command.py,sha256=Oh6WdYgbPKRChrxCoOOEFJHjINF3BIXahmSlYILDfPU,19723
+varvamp/command.py,sha256=9SOR3G9Ew0bJMQt_osWI9QW2Hh671kjez02SwhcnOAY,19749
 varvamp/scripts/__init__.py,sha256=DtRsgfnA60BvccKMuD-Ueo1UpxaeANUfIWxGi1xPcV0,46
 varvamp/scripts/alignment.py,sha256=w2I7BnaUt_rh1EaDQaU9Q_dmL_q019NtnaC8YMccOgM,7317
 varvamp/scripts/blast.py,sha256=yF1-FjL7-P8sxSO6x4Q15gzj_I3LbE8ZfP37w6K3Akk,8434
 varvamp/scripts/consensus.py,sha256=eU5eKU9iXyWeln2BW_LMPI8JiEq0_NXBI0jy4kMOvQg,3375
 varvamp/scripts/default_config.py,sha256=8yjZ1ugZ5pus1daWRXeQMo-lVIikJXV4_q7FPEkTS64,3766
 varvamp/scripts/get_config.py,sha256=Kvg3YhttUbDeRxjKhfxLjI5JSKpPoiLSEYqcE6rT4Xk,2940
-varvamp/scripts/logging.py,sha256=bk7SA7DtNtKtAJXRGLAQyA0stpLqBmRUydBG0dpoZQ0,20820
+varvamp/scripts/logging.py,sha256=Vg_iojqrTejjKYOwba5VGm0gAi_-cQ5h7WNx0bMnmBc,20820
 varvamp/scripts/param_estimation.py,sha256=jpfmbjCFp23V07af3y1uO2dQZJiWLlyOBa11aMphkko,3694
 varvamp/scripts/primers.py,sha256=zLNs3lWiHS0dBmRIMpi8v0XxWi26JrS43YJW7xhEneM,13422
 varvamp/scripts/qpcr.py,sha256=qO6JTl36WMUZGoz6NHIjRrtMOF87yQedmj-sE85Yay8,14453
 varvamp/scripts/regions.py,sha256=VgYAlZmkC3rT946yhAdxATi-YT4oGcEwAPa0BkkfWIg,3239
-varvamp/scripts/reporting.py,sha256=_o-A-X1OPEB6ym0w-q3Tb0nehVJwuCbsoyspZjzmvs4,22921
-varvamp/scripts/scheme.py,sha256=lB-HT-rFSjN0o1KemzojyMT3jBw-LSrDGakshBXCAR4,16198
-varvamp-1.2.1.dist-info/METADATA,sha256=tqQlI16IR0-DMDNYV1e6Ka9SzbzKbtqF1zV4gwBiVRM,5382
-varvamp-1.2.1.dist-info/WHEEL,sha256=cVxcB9AmuTcXqmwrtPhNK88dr7IR_b6qagTj0UvIEbY,91
-varvamp-1.2.1.dist-info/entry_points.txt,sha256=puzW-basyBexZT4JVRUfUEqobvFmEyfqRQaqFjp7rB0,49
-varvamp-1.2.1.dist-info/top_level.txt,sha256=11oVwE3SBUB9aTmvpvEDru95Tc5GZqQikzzFjw2eVGc,8
-varvamp-1.2.1.dist-info/RECORD,,
+varvamp/scripts/reporting.py,sha256=KGy3Nicb8q3dgyzkNqYQa6BYlUsWzERv_5O2TDkpY5s,22792
+varvamp/scripts/scheme.py,sha256=8pg8n5WMVUx23Pb7fNjUlLIpV_TqszVIAvvlcgqg8mo,16198
+varvamp-1.2.2.dist-info/METADATA,sha256=dgoF0hWirC5BQAQwMUOaMiLV6K2C16yUTotKOAacajo,5573
+varvamp-1.2.2.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+varvamp-1.2.2.dist-info/entry_points.txt,sha256=puzW-basyBexZT4JVRUfUEqobvFmEyfqRQaqFjp7rB0,49
+varvamp-1.2.2.dist-info/top_level.txt,sha256=11oVwE3SBUB9aTmvpvEDru95Tc5GZqQikzzFjw2eVGc,8
+varvamp-1.2.2.dist-info/RECORD,,

{varvamp-1.2.1.dist-info → varvamp-1.2.2.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (74.1.2)
+Generator: setuptools (75.8.0)
 Root-Is-Purelib: true
 Tag: py3-none-any