Build a RefSeq Homology Alignment and FHIR Molecular Evidence Indexer

Align homologous transcript sequences retrieved from NCBI Datasets using Biopython pairwise routines. Generate FHIR R5 Molecular Observation entities capturing sequence divergence metrics.

Comparative Genomics Researcher

Ingest FASTA files for homologous eukaryotic gene packages fetched via NCBI Datasets taxonomy flags.

Execute pairwise sequence alignments across homologs using Biopython alignment classes.

Calculate sequence identity percentages, insertion count, and deletion gaps across CDS regions.

Export divergence summaries as FHIR Observation resources linked to ResearchStudy identifiers.

Web workspace with sequence comparison matrix viewer and FHIR bundle inspection drawer.

Download ACRV1 gene sequences across human and murine taxa via NCBI Datasets command mock directory. Starter examples below are synthetic. They are not downloaded records.

NCBI Datasets gene download Leverages NCBI Datasets taxon flag '--taxon human' and multi-accession download conventions for homologous gene extraction.

dspy, biopython, qdrant. Optimizes prompt pipelines that convert complex alignment matrices into plain-language molecular divergence summaries. Fast remote instruction-tuned language model with 16k context window for multi-sequence summaries.

Submit runnable source, pinned dependencies, setup commands, a working interface, example outputs, and a recorded end-to-end demo. Include tests and an action log.

Calculates exact pairwise identity percentage over Biopython sequence alignments.

Structures divergence scores into compliant FHIR Molecular Observation code-value pairs.

Detects frameshift indels and marks affected sequence records with distinct warning flags.

```json [ { "input": { "query_gene": "ACRV1", "query_seq": "ATGAAACTCCTTCTCCTGGCTTCAGTGTTGTTTTCATTGGTGTTCACCCAAGAATCT", "target_seq": "ATGAAACTCCTTCTCCTGGCTTCAGTGTTGTTTTCATTGGTGTTCACCCAAGAATCT" }, "expected_output": { "alignment_score": 57, "percent_identity": 100, "gap_count": 0, "frameshift_detected": false, "fhir_observation_code": "sequence-homology" } }, { "input": { "query_gene": "ACRV1", "query_seq": "ATGAAACTCCTTCTCCTGGCTTCAGTGTTGTTTTCATTGGTGTTCACCCAAGAATCT", "target_seq": "ATGAAACTCCTTCTCCTGGCTTC---TTGTTTTCATTGGTGTTCACCCAAGAATCT" }, "expected_output": { "alignment_score": 51, "percent_identity": 94.74, "gap_count": 3, "frameshift_detected": true, "fhir_observation_code": "sequence-homology" } } ] ```

Operates strictly on de-identified public reference genomes; not used for diagnostic variant interpretation. Use mock external writes. Require human review before any real deployment.