loqusdb 2.7.16__py3-none-any.whl → 2.7.17__py3-none-any.whl

loqusdb/__init__.py

@@ -4,7 +4,7 @@ from pymongo import ASCENDING, IndexModel
 
 logger = logging.getLogger(__name__)
 
-__version__ = "2.7.16"
+__version__ = "2.7.17"
 
 INDEXES = {
     "variant": [

loqusdb/commands/load.py

@@ -88,6 +88,13 @@ def validate_profile_threshold(ctx, param, value):
     callback=validate_profile_threshold,
     help="profile hamming distance to store similar individuals (0-1)",
 )
+@click.option(
+    "--snv-gq-only",
+    is_flag=True,
+    default=False,
+    show_default=True,
+    help="Apply GQ threshold only to SNV variants",
+)
 @click.pass_context
 def load(
     ctx,
@@ -104,6 +111,7 @@ def load(
     hard_threshold,
     soft_threshold,
     qual_gq,
+    snv_gq_only,
 ):
     """Load the variants of a case
 
@@ -145,6 +153,7 @@ def load(
             skip_case_id=skip_case_id,
             case_id=case_id,
             gq_threshold=gq_threshold,
+            snv_gq_only=snv_gq_only,
             qual_gq=qual_gq,
             max_window=max_window,
             profile_file=variant_profile_path,

loqusdb/utils/load.py

@@ -31,6 +31,7 @@ def load_database(
     family_type="ped",
     skip_case_id=False,
     gq_threshold=None,
+    snv_gq_only=False,
     qual_gq=False,
     case_id=None,
     max_window=3000,
@@ -88,11 +89,13 @@ def load_database(
             adapter, profiles, hard_threshold=hard_threshold, soft_threshold=soft_threshold
         )
 
-    # If a gq threshold is used the variants needs to have GQ
+    # If a gq threshold is used the variants need to have GQ (only SNVs if snv_gq_only)
     for _vcf_file in vcf_files:
-        # Get a cyvcf2.VCF object
-        vcf = get_vcf(_vcf_file)
+        is_sv = _vcf_file == sv_file
+        if snv_gq_only and is_sv:
+            continue  # skip GQ check for SV VCF
 
+        vcf = get_vcf(_vcf_file)
         if gq_threshold and not vcf.contains("GQ") and not qual_gq:
             LOG.warning("Set gq-threshold to 0 or add info to vcf {0}".format(_vcf_file))
             raise SyntaxError("GQ is not defined in vcf header")
@@ -144,7 +147,7 @@ def load_database(
                 vcf_obj=vcf_obj,
                 case_obj=case_obj,
                 skip_case_id=skip_case_id,
-                gq_threshold=gq_threshold,
+                gq_threshold=gq_threshold if not snv_gq_only or variant_type == "snv" else None,
                 qual_gq=qual_gq,
                 max_window=max_window,
                 variant_type=variant_type,

{loqusdb-2.7.16.dist-info → loqusdb-2.7.17.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: loqusdb
-Version: 2.7.16
+Version: 2.7.17
 Summary: A simple observation count database
 License: MIT
 Author: Your Name

{loqusdb-2.7.16.dist-info → loqusdb-2.7.17.dist-info}/RECORD

@@ -1,4 +1,4 @@
-loqusdb/__init__.py,sha256=D1R2iXsppIzZJxfYNMSdRlrzsMKf-4z-tSMYAFWtEPw,1688
+loqusdb/__init__.py,sha256=OKW6E-XujZmZRrXRSOInDNnGqybJIkbeIBvnFcWqhuo,1688
 loqusdb/__main__.py,sha256=8FGKySAGaWSzAYMj6HRsxeyiME3V01Idt7HrmN7pSYY,397
 loqusdb/build_models/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 loqusdb/build_models/case.py,sha256=AByutEYK2N3kS9JFvyZfPKNZdCpZHCSD0nNHAgaU1Cs,4127
@@ -10,7 +10,7 @@ loqusdb/commands/cli.py,sha256=wJD5S1BoCxtRTAobd1QtmQpzlngJg-mt1nsyD92fDD4,3176
 loqusdb/commands/delete.py,sha256=R6ysHKSMw1mmL4ZbktoUIKzdzDLQ3314YPYhIy1myic,1979
 loqusdb/commands/export.py,sha256=HKoRzUo_BHNOdw_TcKUId9TTowi8VJVGqnuDlK-FqFE,3531
 loqusdb/commands/identity.py,sha256=KLA9c8e6cJFDxtqIa1G6zdHTHK1sz2b3v1Utdtik_4k,787
-loqusdb/commands/load.py,sha256=sDGCOGbukdfebLn1iBGyWoI6t3tijLzywP1tg3WheQI,4431
+loqusdb/commands/load.py,sha256=6rrt_XIdUHhNc-HisMRS7G-bWkgYhqggDjbEQubYmUQ,4635
 loqusdb/commands/load_profile.py,sha256=cflCbF9f77_HCH8xPnN8zSSocvIffRMnC2LPE0j7Xq8,3336
 loqusdb/commands/migrate.py,sha256=2C8YL-zVqnpnqg3JIyUr0rbVnb8-AGPVWNhicHnPKLo,667
 loqusdb/commands/restore.py,sha256=eqPX0yao0IAYS5SbjCdlsfSJRBbRByBLISUU2hTzqqs,1492
@@ -43,14 +43,14 @@ loqusdb/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 loqusdb/utils/annotate.py,sha256=cPNWlhsv6yoe3lxNfa9DytO5eACuM_mOJJw_mglVMN0,2646
 loqusdb/utils/case.py,sha256=aeTvyACJTDjzl-aOjAZaUzFMLisgFKMfcoXSvNAZz4s,2168
 loqusdb/utils/delete.py,sha256=-ddBM_QXKzlUN6egEJggKzXX1P-WEdi92HgaD1DJRtg,4843
-loqusdb/utils/load.py,sha256=ZdsnmHTVGOjwG24Fm5gcUu7tcLwYSqvh3W9j6ixqyF0,8239
+loqusdb/utils/load.py,sha256=1UVB5jUy1DESZdWnRrMcyfFiUgqmj9l-Bm0DZOHetfg,8428
 loqusdb/utils/migrate.py,sha256=9Q6kdIi9TpFVzDYptlEE8RqPPS5wyzfM3F8egzmmBBk,1113
 loqusdb/utils/profiling.py,sha256=3OizF7CpYvSl9kyl2g4KGJxbIRUqWfmfLxn3843XYDk,9164
 loqusdb/utils/update.py,sha256=1edJG-u24FgOSxyXAQEiyTG4IyK-Uo3lSIl5qyzcXsI,4433
 loqusdb/utils/variant.py,sha256=U6nMZRUf5NDDQ74nG0HBCLMnFQVgFAT6eHll_F2uiwc,2087
 loqusdb/utils/vcf.py,sha256=ybmrTBEPYa0FbUXo8ttlwATk13RnKjX9eIDbRDwCiVE,5175
-loqusdb-2.7.16.dist-info/LICENSE,sha256=urpFcJXw3elN9kV2fFutc-lXegjuu2lqP_GSy8_CAbs,1054
-loqusdb-2.7.16.dist-info/METADATA,sha256=Aayv73fzmF_pCfXHNoRZEzen9CbZqZccc5pfwzMXaCE,5321
-loqusdb-2.7.16.dist-info/WHEEL,sha256=fGIA9gx4Qxk2KDKeNJCbOEwSrmLtjWCwzBz351GyrPQ,88
-loqusdb-2.7.16.dist-info/entry_points.txt,sha256=wFoWzEFjsSgXkj9FMQA8C9ihZoJ9R1XvbGuX9hEEI6E,52
-loqusdb-2.7.16.dist-info/RECORD,,
+loqusdb-2.7.17.dist-info/LICENSE,sha256=urpFcJXw3elN9kV2fFutc-lXegjuu2lqP_GSy8_CAbs,1054
+loqusdb-2.7.17.dist-info/METADATA,sha256=WH2QGz7PhKif97YZxHnOOxmiaN96ZAd_zGCfUnm2GGU,5321
+loqusdb-2.7.17.dist-info/WHEEL,sha256=b4K_helf-jlQoXBBETfwnf4B04YC67LOev0jo4fX5m8,88
+loqusdb-2.7.17.dist-info/entry_points.txt,sha256=wFoWzEFjsSgXkj9FMQA8C9ihZoJ9R1XvbGuX9hEEI6E,52
+loqusdb-2.7.17.dist-info/RECORD,,

{loqusdb-2.7.16.dist-info → loqusdb-2.7.17.dist-info}/WHEEL

@@ -1,4 +1,4 @@
 Wheel-Version: 1.0
-Generator: poetry-core 2.1.2
+Generator: poetry-core 2.1.3
 Root-Is-Purelib: true
 Tag: py3-none-any