We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 2bd58d2 commit 6dc2781Copy full SHA for 6dc2781
src/dcd_mapping/align.py
@@ -363,6 +363,11 @@ def align(
363
msg = f"BLAT result {target_label} matches multiple target gene names in scoreset {scoreset_metadata.urn}"
364
target_gene = scoreset_metadata.target_genes[target_label]
365
alignment_results[target_label] = _get_best_match(blat_result, target_gene)
366
+ # confirm that there is an alignment result for each target gene
367
+ for target_gene in scoreset_metadata.target_genes:
368
+ if target_gene not in alignment_results:
369
+ msg = f"No BLAT result found for target gene {target_gene} in scoreset {scoreset_metadata.urn}"
370
+ raise AlignmentError(msg)
371
return alignment_results
372
373
0 commit comments