Fix g123 test failures on Adir1/Amin1 by falling back to segregating sites when no phasing analysis is available

Author: blankkiri

malariagen_data/amin1.py

@@ -105,7 +105,7 @@ def __init__(
             site_filters_analysis=site_filters_analysis,
             discordant_read_calls_analysis=discordant_read_calls_analysis,
             default_site_mask="minimus",
-            default_phasing_analysis="minimus_noneyet",
+            default_phasing_analysis=None,
             default_coverage_calls_analysis="minimus_noneyet",
             bokeh_output_notebook=bokeh_output_notebook,
             results_cache=results_cache,

tests/anoph/test_g123.py

@@ -6,6 +6,8 @@
 
 from malariagen_data import af1 as _af1
 from malariagen_data import ag3 as _ag3
+from malariagen_data import adir1 as _adir1
+from malariagen_data import amin1 as _amin1
 from malariagen_data.anoph.g123 import AnophelesG123Analysis
 
 
@@ -56,6 +58,44 @@ def af1_sim_api(af1_sim_fixture):
     )
 
 
+@pytest.fixture
+def adir1_sim_api(adir1_sim_fixture):
+    return AnophelesG123Analysis(
+        url=adir1_sim_fixture.url,
+        public_url=adir1_sim_fixture.url,
+        config_path=_adir1.CONFIG_PATH,
+        major_version_number=_adir1.MAJOR_VERSION_NUMBER,
+        major_version_path=_adir1.MAJOR_VERSION_PATH,
+        pre=True,
+        gff_gene_type="protein_coding_gene",
+        gff_gene_name_attribute="Note",
+        gff_default_attributes=("ID", "Parent", "Note", "description"),
+        default_site_mask="dirus",
+        results_cache=adir1_sim_fixture.results_cache_path.as_posix(),
+        taxon_colors=_adir1.TAXON_COLORS,
+        default_phasing_analysis=None,
+    )
+
+
+@pytest.fixture
+def amin1_sim_api(amin1_sim_fixture):
+    return AnophelesG123Analysis(
+        url=amin1_sim_fixture.url,
+        public_url=amin1_sim_fixture.url,
+        config_path=_amin1.CONFIG_PATH,
+        major_version_number=_amin1.MAJOR_VERSION_NUMBER,
+        major_version_path=_amin1.MAJOR_VERSION_PATH,
+        pre=True,
+        gff_gene_type="protein_coding_gene",
+        gff_gene_name_attribute="Note",
+        gff_default_attributes=("ID", "Parent", "Note", "description"),
+        default_site_mask="minimus",
+        results_cache=amin1_sim_fixture.results_cache_path.as_posix(),
+        taxon_colors=_amin1.TAXON_COLORS,
+        default_phasing_analysis=None,
+    )
+
+
 # N.B., here we use pytest_cases to parametrize tests. Each
 # function whose name begins with "case_" defines a set of
 # inputs to the test functions. See the documentation for
@@ -76,6 +116,14 @@ def case_af1_sim(af1_sim_fixture, af1_sim_api):
     return af1_sim_fixture, af1_sim_api
 
 
+def case_adir1_sim(adir1_sim_fixture, adir1_sim_api):
+    return adir1_sim_fixture, adir1_sim_api
+
+
+def case_amin1_sim(amin1_sim_fixture, amin1_sim_api):
+    return amin1_sim_fixture, amin1_sim_api
+
+
 def check_g123_gwss(*, api, g123_params):
     # Run main gwss function under test.
     x, g123 = api.g123_gwss(**g123_params)
@@ -115,6 +163,10 @@ def test_g123_gwss_with_default_sites(fixture, api: AnophelesG123Analysis):
 
 @parametrize_with_cases("fixture,api", cases=".")
 def test_g123_gwss_with_phased_sites(fixture, api: AnophelesG123Analysis):
+    # Skip if this dataset has no phasing analyses (e.g., Adir1, Amin1).
+    if not api.phasing_analysis_ids:
+        pytest.skip("No phasing analyses available for this dataset.")
+
     # Set up test parameters.
     all_sample_sets = api.sample_sets()["sample_set"].to_list()
     g123_params = dict(
@@ -182,6 +234,10 @@ def test_g123_gwss_with_bad_sites(fixture, api: AnophelesG123Analysis):
 
 @parametrize_with_cases("fixture,api", cases=".")
 def test_g123_calibration(fixture, api: AnophelesG123Analysis):
+    # Skip if this dataset has no phasing analyses (e.g., Adir1, Amin1).
+    if not api.phasing_analysis_ids:
+        pytest.skip("No phasing analyses available for this dataset.")
+
     # Set up test parameters.
     all_sample_sets = api.sample_sets()["sample_set"].to_list()
     window_sizes = np.random.randint(100, 500, size=random.randint(2, 5)).tolist()