Skip to content

Commit c52d7b6

Browse files
authored
Merge pull request #58 from biodiversitycellatlas/dev
Include alevin data in reporting dashboard
2 parents a54f509 + 40d046b commit c52d7b6

9 files changed

Lines changed: 442 additions & 211 deletions

File tree

bin/dashboard_mappingstats.py

Lines changed: 8 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -497,7 +497,7 @@ def parse_alevin_quant(quant_path: Path) -> Dict[str, Optional[int]]:
497497

498498
def summarize_cell_meta(cell_meta_path: Path) -> Dict[str, Optional[int]]:
499499
"""
500-
Summarize alevin-fry cell_meta.tsv-like file.
500+
Summarize alevin-fry cell_meta.tsv file.
501501
502502
Returns:
503503
- n_cells
@@ -551,9 +551,9 @@ def parse_alevinfry_sample(sample_root: Path) -> Dict[str, object]:
551551
Parse one alevin-fry sample directory under mapping_alevin/<sample_id>.
552552
553553
Uses:
554-
- <sample_id>_run/aux_info/meta_info.json total & mapped reads
555-
- <sample_id>_counts/quant.json N cells, Total Genes Detected
556-
- cell_meta.tsv per-cell summaries
554+
- <sample_id>_run/aux_info/meta_info.json : total & mapped reads
555+
- <sample_id>_counts/quant.json : N cells, Total Genes Detected
556+
- cell_meta.tsv : per-cell summaries
557557
"""
558558
sample_name = sample_root.name
559559
row: Dict[str, object] = {
@@ -573,12 +573,12 @@ def parse_alevinfry_sample(sample_root: Path) -> Dict[str, object]:
573573
run_dir = first_existing(d for d in sample_root.glob("*_run") if d.is_dir())
574574
counts_dir = first_existing(d for d in sample_root.glob("*_counts") if d.is_dir())
575575

576-
# ---- meta_info.json total & mapped reads ----
576+
# ---- meta_info.json : total & mapped reads ----
577577
meta_candidates: List[Path] = []
578578
if run_dir is not None:
579579
meta_candidates.extend(
580580
[
581-
run_dir / "aux_info" / "meta_info.json", # your real path
581+
run_dir / "aux_info" / "meta_info.json",
582582
run_dir / "meta_info.json",
583583
run_dir / "quant" / "meta_info.json",
584584
]
@@ -606,7 +606,7 @@ def parse_alevinfry_sample(sample_root: Path) -> Dict[str, object]:
606606
if uniq_frac is not None:
607607
row["% uniquely mapped reads"] = convert_to_pct(uniq_frac)
608608

609-
# ---- quant.json in *_counts N cells, Total Genes Detected ----
609+
# ---- quant.json in *_counts : N cells, Total Genes Detected ----
610610
quant_candidates: List[Path] = []
611611
if counts_dir is not None:
612612
quant_candidates.extend(
@@ -625,7 +625,7 @@ def parse_alevinfry_sample(sample_root: Path) -> Dict[str, object]:
625625
if quant_summary.get("num_genes") is not None:
626626
row["Total Genes Detected"] = quant_summary["num_genes"]
627627

628-
# ---- cell_meta.tsv per-cell summaries (mean / medians) ----
628+
# ---- cell_meta.tsv : per-cell summaries (mean / medians) ----
629629
cell_meta_candidates: List[Path] = []
630630
if counts_dir is not None:
631631
cell_meta_candidates.extend(

0 commit comments

Comments
 (0)