scout-browser 4.85__py3-none-any.whl → 4.86__py3-none-any.whl

scout/demo/drop/fraser_top_hits_clinical.tsv

@@ -0,0 +1,5 @@
+hgnc_id	geneID	hgncSymbol	gene_type	gene_name_orig	sampleID	seqnames	start	end	width	strand	type	pValue	psiValue	deltaPsi	counts	totalCounts	meanCounts	meanTotalCounts	nonsplitCounts	nonsplitProportion	nonsplitProportion_99quantile	annotatedJunction	pValueGene	padjustGene	PAIRED_END	DNA_ID	DROP_GROUP	SPLICE_COUNTS_DIR	HPO_TERMS	GENE_COUNTS_FILE	GENE_ANNOTATION	GENOME	isExternal	potentialImpact	causesFrameshift	UTR_overlap	blacklist
+2439	ENSG00000119535.18	CSF3R	protein_coding	CSF3R	ADM1059A2	chr1	36479517	Imp	1961	-	jaccard	1.6652e-06	0.49	-0.39	127	258	4237.77	4570.55	3	0.01	0.02	both	9.9912e-06	0.013423	True		outrider,fraser						False	annotatedIntron_reducedUsage	unlikely	5'-UTR	False
+4831	ENSG00000213934.9	HBG1	protein_coding	HBG1	ADM1059A2	chr11	5248488	5254291	5804	-	jaccard	1.918e-12	0.36	0.35	35	96	22.84	6902.91	0	0.0	0.0	end	5.7541e-12	1.2885e-08	True		outrider,fraser						False	exonSkipping	inconclusive	3'-UTR	False
+4832	ENSG00000196565.15	HBG2	protein_coding	HBG2	ADM1059A2	chr11	5248488	5254291	5804	-	jaccard	1.918e-12	0.36	0.35	35	96	22.84	6902.91	0	0.0	0.0	end	3.836e-12	1.2885e-08	True		outrider,fraser						False	exonSkipping	inconclusive	3'-UTR	False
+17284	ENSG00000213934.9	POT1	protein_coding	POT1	ADM1059A2	chr7	124532319	124532434	115	-	jaccard	1.918e-12	0.36	0.35	35	96	22.84	6902.91	0	0.0	0.0	end	5.7541e-12	1.2885e-08	True		outrider,fraser						False	exonSkipping	inconclusive	3'-UTR	False

scout/demo/drop/outrider_top_hits_clinical.tsv

@@ -0,0 +1,10 @@
+hgnc_id	seqnames	start	end	strand	geneID	hgncSymbol	gene_type	gene_name_orig	sampleID	pValue	padjust	zScore	l2fc	rawcounts	normcounts	meanCorrected	theta	aberrant	AberrantBySample	AberrantByGene	padj_rank	FDR_set	foldChange
+25415	chr4	88257620	88284769	-	ENSG00000163644.15	PPM1K	protein_coding	PPM1K	ADM1059A2	0.0016124374690447165	1.0	-5.92	-0.9	27	317.46	601.46	139.77	False	4.0	0.0	6110.0	transcriptome-wide	0.54
+10019	chr6	3063824	3115187	+	ENSG00000137275.16	RIPK1	protein_coding	RIPK1	ADM1059A2	0.0009997468998232232	1.0	6.33	0.53	104	1891.09	1308.4	547.35	False	4.0	0.0	6110.0	transcriptome-wide	1.44
+4827	chr11	5225464	5229395	-	ENSG00000244734.4	HBB	protein_coding	HBB	ADM1059A2	3.417496739472308e-25	4.169327890311004e-20	-12.61	-10.14	61	2625.42	3016434.64	9.21	True	4.0	1.0	1.0	transcriptome-wide	0.0
+4831	chr11	5248269	5249857	-	ENSG00000213934.9	HBG1	protein_coding	HBG1	ADM1059A3	0.0016205310428587193	1.0	2.23	2.86	73	151.01	19.19	0.91	False	4.0	0.0	6110.0	transcriptome-wide	7.26
+16860	chr12	108522214	108561400	-	ENSG00000075856.12	SART3	protein_coding	SART3	ADM1059A2	0.0014636643867152977	1.0	7.02	0.51	97	1412.71	989.2	1000.0	False	4.0	0.0	6110.0	transcriptome-wide	1.42
+4824	chr16	172876	173710	+	ENSG00000188536.13	HBA2	protein_coding	HBA2	ADM1059A2	2.15484140654196e-24	1.2729742965811128e-19	-12.53	-10.01	262	6781.29	7019330.9	8.97	True	4.0	1.0	2.5	transcriptome-wide	0.0
+4823	chr16	176680	177522	+	ENSG00000206172.8	HBA1	protein_coding	HBA1	ADM1059A2	3.1302782768232926e-24	1.2729742965811128e-19	-12.5	-9.97	246	6497.61	6540277.08	8.95	True	4.0	1.0	2.5	transcriptome-wide	0.0
+9543	chr17	4796144	4798502	+	ENSG00000142507.10	PSMB6	protein_coding	PSMB6	ADM1059A2	0.0010771639306525651	1.0	-10.54	-1.14	13	342.89	783.82	1000.0	False	4.0	0.0	6110.0	transcriptome-wide	0.45
+17284	chr7	124532319	124532434	-	ENSG00000213934.9	POT1	protein_coding	POT1	ADM1059A2	0.0016205310428587193	1.0	2.23	2.86	73	151.01	19.19	0.91	False	4.0	0.0	6110.0	transcriptome-wide	7.26

scout/load/setup.py

@@ -51,12 +51,12 @@ def setup_scout(
 
     WARNING: If the instance is populated all collections will be deleted
 
-    Build insert a institute and an admin user.
-    There are multiple sources of information that is used by scout and that needs to exist for
-    scout to work proper.
+    Build and insert an institute and an admin user.
 
+    Multiple sources of information that are used by scout need to exist for
+    scout to work properly.
     Genes:
-         Scout uses HGNC as the source for gene identifiers en ensembl as source for coordinates.
+         Scout uses HGNC as the source for gene identifiers and ENSEMBL as source for coordinates.
          Additional information of disease connections for genes if fetched from OMIM.
          Link between hpo terms and genes is fetched from HPO
          For more details check the documentation.

scout/models/case/case_loading_models.py

@@ -15,7 +15,7 @@ except ImportError:
 
 from pydantic import BaseModel, Field, field_validator, model_validator
 
-from scout.constants import ANALYSIS_TYPES, FILE_TYPE_MAP
+from scout.constants import ANALYSIS_TYPES, FILE_TYPE_MAP, OMICS_FILE_TYPE_MAP
 from scout.exceptions import PedigreeError
 from scout.utils.date import get_date
 
@@ -41,15 +41,17 @@ CASE_FILE_PATH_CHECKS = [
     "cnv_report",
     "coverage_qc_report",
     "delivery_report",
+    "exe_ver",
+    "fraser_tsv",
     "gene_fusion_report",
     "gene_fusion_report_research",
     "madeline_info",
     "multiqc",
     "multiqc_rna",
+    "outrider_tsv",
     "peddy_ped",
     "peddy_ped_check",
     "peddy_sex_check",
-    "exe_ver",
     "smn_tsv",
     "reference_info",
     "RNAfusion_inspector",
@@ -59,6 +61,7 @@ CASE_FILE_PATH_CHECKS = [
 ]
 
 VCF_FILE_PATH_CHECKS = FILE_TYPE_MAP.keys()
+OMICS_FILE_PATH_CHECKS = OMICS_FILE_TYPE_MAP.keys()
 
 GENOME_BUILDS = ["37", "38"]
 TRACKS = ["rare", "cancer"]
@@ -153,6 +156,25 @@ class Mitodel(BaseModel):
     ratioppk: Optional[float] = None
 
 
+class OmicsFiles(BaseModel):
+    """Represents multiple kinds of omics files, e.g. RNA expression outliers for aberrant splicing
+    and aberrant expression."""
+
+    fraser: Optional[str] = None
+    fraser_research: Optional[str] = None
+    outrider: Optional[str] = None
+    outrider_research: Optional[str] = None
+
+    @model_validator(mode="before")
+    def validate_file_path(cls, values: Dict) -> "OmicsFiles":
+        """Make sure that VCF file exists on disk."""
+        for item in OMICS_FILE_PATH_CHECKS:
+            item_path: str = values.get(item)
+            if item_path:
+                values[item] = _resource_abs_path(item_path)
+        return values
+
+
 class REViewer(BaseModel):
     alignment: Optional[str] = None
     alignment_index: Optional[str] = None
@@ -392,6 +414,7 @@ class CaseLoader(BaseModel):
     madeline_info: Optional[str] = Field(None, alias="madeline")
     multiqc: Optional[str] = None
     multiqc_rna: Optional[str] = None
+    omics_files: Optional[OmicsFiles] = None
     owner: Optional[str] = None
     peddy_ped: Optional[str] = None  # Soon to be deprecated
     peddy_ped_check: Optional[str] = Field(None, alias="peddy_check")  # Soon to be deprecated

scout/models/omics_variant.py

@@ -0,0 +1,227 @@
+""" OMICS variant
+
+    For potentially causative variants that are not yet in ClinVar
+    and have yet not been marked causative in any existing case.
+
+"""
+
+import logging
+from datetime import datetime
+from typing import List, Optional
+
+from pydantic import BaseModel, Field, field_validator, model_validator
+
+LOG = logging.getLogger(__name__)
+
+
+class OmicsVariantLoader(BaseModel):
+    """Omics variants loader
+    OmicsVariants are e.g. RNA expression outliers as identified by the DROP pipeline.
+
+    Variable names are as found in the original files, plus a set common to all mixed in after file parsing,
+    but before model validation by this class.
+
+    The serialisation names will be used when dumping the model for e.g. db storage.
+    """
+
+    case_id: str
+    institute: str
+    build: str = "38"
+    variant_type: str = "clinical"
+    category: str  # eg "outlier"
+    sub_category: str  # eg "splicing"
+    date: datetime = datetime.now()
+    display_name: str
+    omics_variant_id: str
+
+    # DROP Fraser and Outrider outlier TSVs
+
+    # sample id is mandatory: each row pertains to one outlier event in one individual as compared to others
+    # In the db object, this will be replaced with a "samples" array of individual dict.
+    sampleID: str
+
+    # outlier variants must identify the gene they pertain to, primarily with an hgnc_id
+    hgnc_ids: Optional[List[int]] = Field(alias="hgnc_id", serialization_alias="hgnc_ids")
+    geneID: Optional[str]
+
+    hgnc_symbols: Optional[List[str]] = Field(
+        alias="hgncSymbol", serialization_alias="hgnc_symbols"
+    )
+    gene_name_orig: Optional[str]
+
+    gene_type: Optional[str]
+
+    # coordinates if applicable
+    chromosome: Optional[str] = Field(alias="seqnames", serialization_alias="chromosome")
+    position: Optional[int] = Field(alias="start", serialization_alias="position")
+    end: Optional[int]
+    width: Optional[int] = None
+    strand: Optional[str] = None
+
+    p_value: Optional[float] = Field(alias="pValue", serialization_alias="p_value", default=None)
+
+    # Fraser specific
+    type: Optional[str] = None
+    psi_value: Optional[float] = Field(
+        alias="psiValue", serialization_alias="psi_value", default=None
+    )
+    delta_psi: Optional[float] = Field(
+        alias="deltaPsi", serialization_alias="delta_psi", default=None
+    )
+    counts: Optional[int] = None
+    total_counts: Optional[int] = Field(
+        alias="totalCounts", serialization_alias="total_counts", default=None
+    )
+    mean_counts: Optional[float] = Field(
+        alias="meanCounts", serialization_alias="mean_counts", default=None
+    )
+    mean_total_counts: Optional[float] = Field(
+        alias="meanTotalCounts", serialization_alias="mean_total_counts", default=None
+    )
+    nonsplit_counts: Optional[int] = Field(
+        alias="nonsplitCounts", serialization_alias="nonsplit_counts", default=None
+    )
+    nonsplit_proportion: Optional[float] = Field(
+        alias="nonsplitProportion", serialization_alias="nonsplit_proportion", default=None
+    )
+    nonsplit_proportion_99quantile: Optional[float] = Field(
+        alias="nonsplitProportion_99quantile",
+        serialization_alias="nonsplit_proportion_99quantile",
+        default=None,
+    )
+    annotated_junction: Optional[str] = Field(
+        alias="annotatedJunction", serialization_alias="annotated_junction", default=None
+    )
+    p_value_gene: Optional[float] = Field(
+        alias="pValueGene", serialization_alias="p_value_gene", default=None
+    )
+    p_adjust_gene: Optional[float] = Field(
+        alias="padjustGene", serialization_alias="p_adjust_gene", default=None
+    )
+    paired_end: Optional[str] = Field(
+        alias="PAIRED_END", serialization_alias="paired_end", default=None
+    )
+    is_external: Optional[bool] = Field(
+        alias="isExternal", serialization_alias="is_external", default=None
+    )
+    potential_impact: Optional[str] = Field(
+        alias="potentialImpact", serialization_alias="potential_impact", default=None
+    )
+    causes_frameshift: Optional[str] = Field(
+        alias="causesFrameshift", serialization_alias="causes_frameshift", default=None
+    )
+    utr_overlap: Optional[str] = Field(
+        alias="UTR_overlap", serialization_alias="utr_overlap", default=None
+    )
+
+    # Outrider specific
+    padjust: Optional[float] = None
+    zscore: Optional[float] = Field(alias="zScore", serialization_alias="zscore", default=None)
+    l2fc: Optional[float] = None
+    rawcounts: Optional[int] = None
+    normcounts: Optional[float] = None
+    meanCorrected: Optional[float] = None
+    theta: Optional[float] = None
+    aberrant: Optional[bool] = None
+    aberrant_by_sample: Optional[float] = Field(
+        alias="aberrantBySample", serialization_alias="aberrant_by_sample", default=None
+    )
+    aberrant_by_gene: Optional[float] = Field(
+        alias="aberrantByGene", serialization_alias="aberrant_by_gene", default=None
+    )
+    padj_rank: Optional[float] = None
+    fdr_set: Optional[str] = Field(alias="FDR_set", serialization_alias="fdr_set", default=None)
+    fold_change: Optional[float] = Field(
+        alias="foldChange", serialization_alias="fold_change", default=None
+    )
+
+    @field_validator("chromosome")
+    def strip_chr(cls, chrom: str) -> str:
+        """We store chromosome names without a chr prefix internally."""
+        return chrom.lstrip("chr")
+
+    @model_validator(mode="before")
+    def ensure_end(cls, values):
+        """End is not always set, but sometimes width is.
+        Sometimes Imp is given as end. Worst case we default to width 1."""
+        end_guess = int(values.get("start")) + int(values.get("width", 1))
+        if "end" not in values:
+            values["end"] = end_guess
+
+        if isinstance(values["end"], str):
+            if values["end"].isdigit():
+                values["end"] = int(values["end"])
+            if values["end"] == "Imp":
+                # imprecise?
+                values["end"] = end_guess
+
+        return values
+
+    @model_validator(mode="before")
+    def genes_become_lists(cls, values):
+        """HGNC ids and gene symbols are found one on each line in DROP tsvs.
+        Convert to a list with a single member in omics_variants for storage."""
+
+        if "hgnc_id" in values:
+            values["hgnc_id"] = [int(values.get("hgnc_id"))]
+
+        if "hgncSymbol" in values:
+            values["hgncSymbol"] = [str(values.get("hgncSymbol"))]
+
+        return values
+
+    @model_validator(mode="before")
+    def set_display_name(cls, values) -> "OmicsVariantLoader":
+        """Set a free text qualification, depending on the kind of variant."""
+
+        values["display_name"] = "_".join(
+            [
+                values.get("hgncSymbol"),
+                values.get("category"),
+                values.get("sub_category"),
+                get_qualification(values=values),
+                values.get("seqnames"),  # chrom, unserialised
+                str(values.get("start")),
+                str(values.get("end")),
+                values.get("variant_type"),
+            ]
+        )
+        return values
+
+    @model_validator(mode="before")
+    def set_omics_variant_id(cls, values) -> "OmicsVariantLoader":
+        """Set OMICS variant id based on the kind of variant."""
+
+        values["omics_variant_id"] = "_".join(
+            [
+                values.get("seqnames"),  # chrom, unserialised
+                str(values.get("start")),
+                str(values.get("end")),
+                values.get("build"),
+                values.get("hgncSymbol"),
+                values.get("sub_category"),
+                get_qualification(values=values),
+                values.get("variant_type"),
+            ]
+        )
+        return values
+
+    @model_validator(mode="before")
+    def set_sample_display_name(cls, values) -> "OmicsVariantLoader":
+        """Set a display name."""
+        values["display_name"] = values.get(
+            "display_name", values.get("sample_name", values.get("individual_id"))
+        )
+        return values
+
+
+def get_qualification(values: dict) -> str:
+    """Get qualification string for ID and display name.
+    This string further qualifies the kind of omics event,
+    e.g. for an expression outlier it could be 'up' or 'down'."""
+    qualification = "affected"
+    if values.get("sub_category") == "expression":
+        qualification = "up" if float(values.get("zScore", 0)) > 0 else "down"
+    if values.get("sub_category") == "splicing":
+        qualification = values.get("potentialImpact")
+    return qualification

scout/parse/omics_variant/__init__.py

@@ -0,0 +1,11 @@
+from typing import Dict, Iterable, List
+
+from .drop import parse_omics_tsv
+
+OMICS_CATEGORY_PARSER = {"tsv": parse_omics_tsv}
+
+
+def parse_omics_file(omics_lines: Iterable[str], omics_file_type: dict) -> List[Dict[str, str]]:
+    """Call appropriate parser for omics variants file, depending on the file format anticipated."""
+    parser = OMICS_CATEGORY_PARSER[omics_file_type.get("format")]
+    return parser(omics_lines)

scout/parse/omics_variant/drop.py

@@ -0,0 +1,19 @@
+from typing import Dict, Iterable, List
+
+
+def parse_omics_tsv(lines: Iterable[str]) -> List[Dict[str, str]]:
+    """Parse a DROP Outrider or Fraser TSV file."""
+    omics_infos = []
+    header = []
+
+    for i, line in enumerate(lines):
+        line = line.rstrip()
+        if i == 0:
+            # Header line
+            header = line.split("\t")
+            continue
+
+        info = dict(zip(header, line.split("\t")))
+        omics_infos.append(info)
+
+    return omics_infos

scout/parse/variant/callers.py

@@ -22,6 +22,7 @@ def parse_callers(variant, category="snv"):
     """
     relevant_callers = CALLERS[category]
     callers = {caller["id"]: None for caller in relevant_callers}
+    callers_keys = set(callers.keys())
 
     other_info = variant.INFO.get("FOUND_IN")
     svdb_origin = variant.INFO.get("svdb_origin")
@@ -30,10 +31,12 @@ def parse_callers(variant, category="snv"):
     if other_info:
         for info in other_info.split(","):
             called_by = info.split("|")[0]
-            callers[called_by] = "Pass"
+            if called_by in callers_keys:
+                callers[called_by] = "Pass"
     elif svdb_origin:
         for called_by in svdb_origin.split("|"):
-            callers[called_by] = "Pass"
+            if called_by in callers_keys:
+                callers[called_by] = "Pass"
     elif raw_info:
         info = raw_info.split("-")
         for call in info:
@@ -47,7 +50,7 @@ def parse_callers(variant, category="snv"):
                 for caller in callers:
                     if caller in call:
                         callers[caller] = "Filtered"
-            elif call in set(callers.keys()):
+            elif call in callers_keys:
                 callers[call] = "Pass"
 
     if raw_info or svdb_origin or other_info:

scout/parse/variant/frequency.py

@@ -11,8 +11,14 @@ EXAC_KEYS = ["EXACAF"]
 EXAC_MAX_KEYS = ["ExAC_MAX_AF", "EXAC_MAX_AF"]
 
 # gnomAD has both SNV and SV
-GNOMAD_INFO_KEYS = ["GNOMADAF", "GNOMAD_AF", "gnomADg_AF", "gnomad_svAF"]
-GNOMAD_INFO_MAX_KEYS = ["gnomADg_AF_POPMAX", "GNOMADAF_popmax", "GNOMADAF_POPMAX", "GNOMADAF_MAX"]
+GNOMAD_INFO_KEYS = ["GNOMADAF", "GNOMAD_AF", "gnomADg_AF", "gnomad_svAF", "gnomad_af"]
+GNOMAD_INFO_MAX_KEYS = [
+    "gnomADg_AF_POPMAX",
+    "GNOMADAF_popmax",
+    "GNOMADAF_POPMAX",
+    "GNOMADAF_MAX",
+    "gnomad_popmax_af",
+]
 
 # SV
 CLINGEN_BENIGN_KEYS = [
@@ -69,6 +75,7 @@ def parse_frequencies(variant, transcripts):
     # These are SV-specific frequencies
     update_frequency_from_vcf(frequencies, variant, ["left_1000GAF"], "thousand_g_left")
     update_frequency_from_vcf(frequencies, variant, ["right_1000GAF"], "thousand_g_right")
+    update_frequency_from_vcf(frequencies, variant, ["colorsdb_af"], "colorsdb_af")
 
     # Search transcripts CSQ if not found in VCF INFO
     if not frequencies:
@@ -117,6 +124,7 @@ def parse_sv_frequencies(variant: cyvcf2.Variant) -> Dict:
     update_sv_frequency_from_vcf(sv_frequencies, variant, SWEGEN_KEYS, "swegen")
     update_sv_frequency_from_vcf(sv_frequencies, variant, DECIPHER_KEYS, "decipher")
     update_sv_frequency_from_vcf(sv_frequencies, variant, CG_KEYS, "clingen_mip")
+    update_sv_frequency_from_vcf(sv_frequencies, variant, ["colorsdb_af"], "colorsdb_af")
 
     return sv_frequencies
 

scout/server/app.py

@@ -7,7 +7,7 @@ from typing import Dict, Union
 from urllib.parse import parse_qsl, unquote, urlsplit
 
 import coloredlogs
-from flask import Flask, current_app, redirect, request, url_for
+from flask import Flask, redirect, request, url_for
 from flask_cors import CORS
 from flask_login import current_user
 from markdown import markdown as python_markdown
@@ -25,6 +25,7 @@ from .blueprints import (
     institutes,
     login,
     managed_variants,
+    omics_variants,
     panels,
     phenomodels,
     phenotypes,
@@ -98,6 +99,7 @@ def configure_extensions(app):
     extensions.store.init_app(app)
     extensions.login_manager.init_app(app)
     extensions.mail.init_app(app)
+    extensions.clinvar_api.init_app(app)
 
     if app.config.get("SQLALCHEMY_DATABASE_URI"):
         extensions.chanjo_report.init_app(app)
@@ -177,6 +179,7 @@ def register_blueprints(app):
     app.register_blueprint(diagnoses.omim_bp)
     app.register_blueprint(institutes.overview)
     app.register_blueprint(managed_variants.managed_variants_bp)
+    app.register_blueprint(omics_variants.omics_variants_bp)
 
 
 def register_filters(app):

scout/server/blueprints/alignviewers/controllers.py

@@ -1,7 +1,7 @@
 # -*- coding: utf-8 -*-
 import logging
 import os.path
-from typing import Dict
+from typing import Dict, Optional
 
 from flask import flash, session
 from flask_login import current_user
@@ -12,7 +12,7 @@ from scout.server.utils import case_append_alignments, find_index
 from scout.utils.ensembl_rest_clients import EnsemblRestApiClient
 
 LOG = logging.getLogger(__name__)
-CUSTOM_TRACK_NAMES = ["Genes", "ClinVar", "ClinVar CNVs"]
+DEFAULT_TRACK_NAMES = ["Genes", "ClinVar", "ClinVar CNVs"]
 
 
 def check_session_tracks(resource):
@@ -51,33 +51,39 @@ def set_session_tracks(display_obj: dict):
     session["igv_tracks"] = session_tracks
 
 
-def make_igv_tracks(case_obj, variant_id, chrom=None, start=None, stop=None):
+def make_igv_tracks(
+    case_obj: dict,
+    variant_id: str,
+    chrom: Optional[str] = None,
+    start: Optional[int] = None,
+    stop: Optional[int] = None,
+) -> dict:
     """Create a dictionary containing the required tracks for displaying IGV tracks for case or a group of cases
 
     Args:
-        institute_id(str): institute _id
+        institute_id: institute _id
         case_obj(scout.models.Case)
-        variant_id(str): _id of a variant
-        chrom(str/None): requested chromosome [1-22], X, Y, [M-MT]
-        start(int/None): start of the genomic interval to be displayed
-        stop(int/None): stop of the genomic interval to be displayed
+        variant_id: _id of a variant
+        chrom: requested chromosome [1-22], X, Y, [M-MT]
+        start: start of the genomic interval to be displayed
+        stop: stop of the genomic interval to be displayed
 
     Returns:
-        display_obj(dict): A display object containing case name, list of genes, lucus and tracks
+        display_obj: A display object containing case name, list of genes, locus and tracks
     """
     display_obj = {}
     variant_obj = store.variant(document_id=variant_id)
 
+    chromosome = "All"
     if variant_obj:
         # Set display locus
         start = start or variant_obj["position"]
         stop = stop or variant_obj["end"]
+        chrom = chrom or variant_obj.get("chromosome")
 
-        chromosome = chrom or variant_obj.get("chromosome")
-        chromosome = chromosome.replace("MT", "M")
+    if all([start, stop, chrom]):
+        chromosome = chrom.replace("MT", "M")
         display_obj["locus"] = "chr{0}:{1}-{2}".format(chromosome, start, stop)
-    else:
-        chromosome = "All"
 
     # Set genome build for displaying alignments:
     if "38" in str(case_obj.get("genome_build", "37")) or chromosome == "M":
@@ -115,20 +121,27 @@ def make_igv_tracks(case_obj, variant_id, chrom=None, start=None, stop=None):
     return display_obj
 
 
-def make_sashimi_tracks(case_obj, variant_id=None):
+def make_sashimi_tracks(
+    case_obj: dict, variant_id: Optional[str] = None, omics_variant_id: Optional[str] = None
+):
     """Create a dictionary containing the required tracks for a splice junction plot
+    If either a regular variant_id or an omics variant id is passed, set display to a particular locus.
+    Otherwise defaults to whole genome "All" view.
 
-    Args:
-        case_obj(scout.models.Case)
-        variant_id(str) _id of a variant
     Returns:
-        display_obj(dict): A display object containing case name, list of genes, lucus and tracks
+        display_obj(dict): A display object containing case name, list of genes, locus and tracks
     """
     build = "38"  # This feature is only available for RNA tracks in build 38
 
     locus = "All"
+    variant_obj = None
+
     if variant_id:
         variant_obj = store.variant(document_id=variant_id)
+    if omics_variant_id:
+        variant_obj = store.omics_variant(variant_id=omics_variant_id)
+
+    if variant_obj:
         locus = make_locus_from_variant(variant_obj, case_obj, build)
 
     display_obj = {"locus": locus, "tracks": []}
@@ -234,9 +247,9 @@ def set_common_tracks(display_obj, build):
     # Set up IGV tracks that are common for all cases:
     display_obj["reference_track"] = HUMAN_REFERENCE[build]  # Human reference is always present
 
-    # if user settings for igv tracks exist -> use these settings, otherwise display all tracks
+    # if user settings for igv tracks exist -> use these settings, otherwise display default tracks ---> Genes, ClinVar and ClinVar CNVs
     custom_tracks_names = (
-        user_obj.get("igv_tracks") if "igv_tracks" in user_obj else CUSTOM_TRACK_NAMES
+        user_obj.get("igv_tracks") if "igv_tracks" in user_obj else DEFAULT_TRACK_NAMES
     )
 
     display_obj["custom_tracks"] = []
@@ -300,7 +313,6 @@ def set_config_custom_tracks(display_obj: dict, build: str):
     """Set up custom public or private tracks stored in a cloud bucket or locally. These tracks were those specified in the Scout config file.
     Respect user's preferences."""
     user_obj = store.user(email=current_user.email)
-    custom_tracks_names = user_obj.get("igv_tracks")
 
     config_custom_tracks = []
 
@@ -308,8 +320,7 @@ def set_config_custom_tracks(display_obj: dict, build: str):
         build_tracks = config_igv_tracks.tracks.get(build, [])
         for track in build_tracks:
             # Do not display track if user doesn't want to see it
-            if custom_tracks_names and track["name"] not in custom_tracks_names:
-                continue
-            config_custom_tracks.append(track)
+            if "igv_tracks" not in user_obj or track["name"] in user_obj.get("igv_tracks"):
+                config_custom_tracks.append(track)
     if config_custom_tracks:
         display_obj["config_custom_tracks"] = config_custom_tracks

scout/server/blueprints/alignviewers/templates/alignviewers/igv_sashimi_viewer.html

@@ -79,6 +79,7 @@
                                     type: 'wig',
                                     format: "bigwig",
                                     url: "{{ url_for('alignviewers.remote_static', file=track.coverage_wig) }}",
+                                    height: 500,
                                   },
                                   {
                                     type: 'spliceJunctions',
@@ -88,24 +89,27 @@
                                     labelUniqueReadCount: true,
                                     url: "{{ url_for('alignviewers.remote_static', file=track.splicej_bed) }}",
                                     indexURL: "{{ url_for('alignviewers.remote_static', file=track.splicej_bed_index) }}",
-                                    minUniquelyMappedReads: 1
+                                    minUniquelyMappedReads: 1,
+                                    height: 500,
                                   },
                                 ]
                               }, // end of sashimi track with data
-                              { // genes track
-                                name: geneTrack.name,
-                                type: geneTrack.type,
-                                format: geneTrack.format,
-                                sourceType: geneTrack.sourceType,
-                                url: geneTrack.url,
-                                indexURL: geneTrack.indexURL,
-                                displayMode: geneTrack.displayMode,
-                                visibilityWindow: 300000000,
-                                height: 100,
-                                searchable: true,
-                                order: {{counter.loop}},
-                                infoURL: "https://www.ncbi.nlm.nih.gov/gene/?term=$$"
-                              }, // end of genes track
+                              {% if custom_tracks|selectattr("name","equalto", "Genes")|list|length > 0 %}
+                                { // genes track
+                                  name: geneTrack.name,
+                                  type: geneTrack.type,
+                                  format: geneTrack.format,
+                                  sourceType: geneTrack.sourceType,
+                                  url: geneTrack.url,
+                                  indexURL: geneTrack.indexURL,
+                                  displayMode: geneTrack.displayMode,
+                                  visibilityWindow: 300000000,
+                                  height: 100,
+                                  searchable: true,
+                                  order: {{counter.loop}},
+                                  infoURL: "https://www.ncbi.nlm.nih.gov/gene/?term=$$"
+                                }, // end of genes track
+                              {% endif %}
                             {% endif %}
                         {% endfor %}
                     ] // end of tracks