loqusdb 2.6.10__py3-none-any.whl → 2.7.0__py3-none-any.whl

loqusdb/build_models/variant.py

@@ -139,7 +139,7 @@ def get_coords(variant):
     return coordinates
 
 
-def build_variant(variant, case_obj, case_id=None, gq_threshold=None, genome_build=None):
+def build_variant(variant, case_obj, case_id=None, gq_threshold=None, gq_qual=False, genome_build=None):
     """Return a Variant object
 
     Take a cyvcf2 formated variant line and return a models.Variant.
@@ -188,7 +188,13 @@ def build_variant(variant, case_obj, case_id=None, gq_threshold=None, genome_bui
             ind_id = ind_obj["ind_id"]
             # Get the index position for the individual in the VCF
             ind_pos = ind_obj["ind_index"]
-            gq = int(variant.gt_quals[ind_pos])
+
+            if gq_qual:
+                gq = int(variant.QUAL)
+
+            if not gq_qual:
+                gq = int(variant.gt_quals[ind_pos])
+
             if gq_threshold and gq < gq_threshold:
                 continue
 

loqusdb/commands/load.py

@@ -55,6 +55,7 @@ def validate_profile_threshold(ctx, param, value):
 )
 @click.option("--ensure-index", is_flag=True, help="Make sure that the indexes are in place")
 @click.option("--gq-threshold", default=20, show_default=True, help="Threshold to consider variant")
+@click.option("--qual-gq", is_flag=True, default=False, show_default=True, help="Use QUAL tag instead of GQ value for quality filter")
 @click.option(
     "--max-window",
     "-m",
@@ -96,6 +97,7 @@ def load(
     check_profile,
     hard_threshold,
     soft_threshold,
+    qual_gq
 ):
     """Load the variants of a case
 
@@ -137,6 +139,7 @@ def load(
             skip_case_id=skip_case_id,
             case_id=case_id,
             gq_threshold=gq_threshold,
+            qual_gq=qual_gq,
             max_window=max_window,
             profile_file=variant_profile_path,
             hard_threshold=hard_threshold,

loqusdb/utils/load.py

@@ -32,6 +32,7 @@ def load_database(
     family_type="ped",
     skip_case_id=False,
     gq_threshold=None,
+    qual_gq=False,
     case_id=None,
     max_window=3000,
     profile_file=None,
@@ -49,6 +50,7 @@ def load_database(
           family_type(str): Format of family file
           skip_case_id(bool): If no case information should be added to variants
           gq_threshold(int): If only quality variants should be considered
+          qual_gq(bool): Use QUAL field instead of GQ format tag to gate quality
           case_id(str): If different case id than the one in family file should be used
           max_window(int): Specify the max size for sv windows
           check_profile(bool): Does profile check if True
@@ -92,7 +94,7 @@ def load_database(
         # Get a cyvcf2.VCF object
         vcf = get_vcf(_vcf_file)
 
-        if gq_threshold and not vcf.contains("GQ"):
+        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,6 +146,7 @@ def load_database(
                 case_obj=case_obj,
                 skip_case_id=skip_case_id,
                 gq_threshold=gq_threshold,
+                qual_gq=qual_gq,
                 max_window=max_window,
                 variant_type=variant_type,
                 genome_build=genome_build,
@@ -190,6 +193,7 @@ def load_variants(
     case_obj,
     skip_case_id=False,
     gq_threshold=None,
+    qual_gq=False,
     max_window=3000,
     variant_type="snv",
     genome_build=None,
@@ -222,7 +226,7 @@ def load_variants(
     with click.progressbar(vcf_obj, label="Inserting variants", length=nr_variants) as bar:
 
         variants = (
-            build_variant(variant, case_obj, case_id, gq_threshold, genome_build=genome_build)
+            build_variant(variant, case_obj, case_id, gq_threshold, qual_gq, genome_build=genome_build)
             for variant in bar
         )
 

{loqusdb-2.6.10.dist-info → loqusdb-2.7.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: loqusdb
-Version: 2.6.10
+Version: 2.7.0
 Summary: Store observations of vcf variants in a mongodb
 Home-page: https://github.com/moonso/loqusdb
 Author: Måns Magnusson
@@ -29,7 +29,6 @@ Requires-Dist: numpy (==1.21.4)
 Requires-Dist: coloredlogs (==14.0)
 Requires-Dist: pyyaml (==5.4.0)
 Requires-Dist: vcftoolbox (==1.5)
-Requires-Dist: pip (==21.3.1)
 Requires-Dist: setuptools (==59.2.0)
 Requires-Dist: mongo-adapter (>=0.3.3)
 Requires-Dist: ped-parser

{loqusdb-2.6.10.dist-info → loqusdb-2.7.0.dist-info}/RECORD

@@ -4,14 +4,14 @@ loqusdb/log.py,sha256=CDcrCjzs9ef-d5Wg8Q_41bCOZRM5j8PyP06kNcynTj0,1691
 loqusdb/build_models/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 loqusdb/build_models/case.py,sha256=P3sfQkI_fH8u5iqecYWhV866lcHz4upWkepaea5MMIw,4255
 loqusdb/build_models/profile_variant.py,sha256=TbSxfVjESstS_FgbkOW4NQwMQVeTyhn9oc9yPZmDhzI,1021
-loqusdb/build_models/variant.py,sha256=9OhR1mbolMPHyChyDWLFfS3eaiz_N6BfJOWBUWi7JKQ,6558
+loqusdb/build_models/variant.py,sha256=JE3o_htRQjpwwUPxmP5HJ_0Ax8BwKp610uBIo3HTaVA,6671
 loqusdb/commands/__init__.py,sha256=BXAN3UADgqPrkGczzjlLO9GyyQ96dnLnP7n92JlYHgo,603
 loqusdb/commands/annotate.py,sha256=748kImopE5WbaO1nuv3WUgIqezWFSsi7SBeWhOz26-s,1384
 loqusdb/commands/cli.py,sha256=wJD5S1BoCxtRTAobd1QtmQpzlngJg-mt1nsyD92fDD4,3176
 loqusdb/commands/delete.py,sha256=R6ysHKSMw1mmL4ZbktoUIKzdzDLQ3314YPYhIy1myic,1979
 loqusdb/commands/export.py,sha256=0V3S3QU9LKlR13w3KCGfqFliTYiDRCgNwusA27AEvmE,3254
 loqusdb/commands/identity.py,sha256=hzbnvniKgSNEwSeYHsxdNvVlqu_vXeOSLlNFnDXTQjA,779
-loqusdb/commands/load.py,sha256=4ypfiar2Si76R3r5P_ncO7vig4GoQyEnfRm3l79WvGU,4231
+loqusdb/commands/load.py,sha256=cVDdY7meBfcv8nMEGsjAX6aE-SKDOceGqM2vAvXPhko,4407
 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
@@ -40,15 +40,15 @@ 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=Xj1QpgMO_UCs0gEZC_dZz_p0kJPbQyZi9vlTNIOX-mc,8075
+loqusdb/utils/load.py,sha256=rujUk89lL5jO35QfO8NKGJkjEJmN6PQxohKkxtpnw6Y,8252
 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=Lq5x9egVB-3rExBiRceF67TaL4Hp2gGoWMSRBEcnm4Q,2088
 loqusdb/utils/vcf.py,sha256=ybmrTBEPYa0FbUXo8ttlwATk13RnKjX9eIDbRDwCiVE,5175
-loqusdb-2.6.10.dist-info/LICENSE,sha256=urpFcJXw3elN9kV2fFutc-lXegjuu2lqP_GSy8_CAbs,1054
-loqusdb-2.6.10.dist-info/METADATA,sha256=ntsTkcLmU07bQFrYSk958_59He_pwSTxiRAwwsKnkEg,6057
-loqusdb-2.6.10.dist-info/WHEEL,sha256=2wepM1nk4DS4eFpYrW1TTqPcoGNfHhhO_i5m4cOimbo,92
-loqusdb-2.6.10.dist-info/entry_points.txt,sha256=39QklW01vy9ilBLcRPgCP18kN6oKXXLoOK50gZE7Jbs,59
-loqusdb-2.6.10.dist-info/top_level.txt,sha256=lRdRO6hqPhJEjFhfNsbCgVWMztvkYsjiOGK9DAL0UAI,8
-loqusdb-2.6.10.dist-info/RECORD,,
+loqusdb-2.7.0.dist-info/LICENSE,sha256=urpFcJXw3elN9kV2fFutc-lXegjuu2lqP_GSy8_CAbs,1054
+loqusdb-2.7.0.dist-info/METADATA,sha256=Sg2aINOps7cTErVgDl6xszkTG7GWTPrGdwFUQ9onR50,6026
+loqusdb-2.7.0.dist-info/WHEEL,sha256=pkctZYzUS4AYVn6dJ-7367OJZivF2e8RA9b_ZBjif18,92
+loqusdb-2.7.0.dist-info/entry_points.txt,sha256=39QklW01vy9ilBLcRPgCP18kN6oKXXLoOK50gZE7Jbs,59
+loqusdb-2.7.0.dist-info/top_level.txt,sha256=lRdRO6hqPhJEjFhfNsbCgVWMztvkYsjiOGK9DAL0UAI,8
+loqusdb-2.7.0.dist-info/RECORD,,

{loqusdb-2.6.10.dist-info → loqusdb-2.7.0.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.38.4)
+Generator: bdist_wheel (0.40.0)
 Root-Is-Purelib: true
 Tag: py3-none-any