tooluniverse-sequence-retrieval
Retrieve DNA/RNA/protein sequences from NCBI and ENA with disambiguation. Quality hierarchy: RefSeq (NM_/NP_) > RefSeq predicted (XM_/XP_) > GenBank submissions. Use for fetching specific sequences by accession, gene-symbol-to-sequence lookup, transcript-isoform retrieval, and curated-vs-raw-submission preference.
Security Assessment
About tooluniverse-sequence-retrieval
This skill guides an agent through retrieving DNA, RNA, and protein sequences from the NCBI and ENA biological databases using the ToolUniverse toolkit. It solves a recurring problem in bioinformatics work: sequence identifiers are ambiguous, databases overlap, and naive lookups return low-quality or wrong-organism records. The skill encodes a disambiguation-first workflow and an explicit quality hierarchy (curated RefSeq NM_/NP_ accessions rank above predicted RefSeq XM_/XP_, which rank above raw GenBank submissions), so the agent fetches the best available sequence rather than the first one it finds.
The workflow runs through defined phases: optional clarification, gene and organism disambiguation via an accession-prefix decision tree, silent search-and-retrieve against NCBI Nucleotide (with ENA as a fallback for non-RefSeq accessions), and a structured Sequence Profile Report. It documents specific ToolUniverse tools (NCBI_search_nucleotide, NCBI_fetch_accessions, NCBI_get_sequence, ena_get_entry, ena_get_sequence_fasta), fallback chains, MANE Select transcript preference for human isoforms, curation-level tiers, error handling for common failures like ENA 404s on RefSeq accessions, and cross-database reconciliation (e.g. GenBank U00096 mapping to RefSeq NC_000913 for E. coli). Worked examples cover fetching the E. coli K-12 genome and the human BRCA1 gene.
It targets bioinformaticians, computational biologists, and researchers who need reproducible, correctly-disambiguated sequence retrieval for downstream BLAST, alignment, or annotation work. It is a read-only data-retrieval helper — it queries public reference databases and formats results; it does not modify data or run analyses itself.
FAQ
Which databases and identifiers does it work with?
It retrieves from NCBI (Nucleotide) and ENA/EMBL. RefSeq accessions (NC_/NM_/NR_/NP_/XM_) are NCBI-only, while GenBank-format accessions are mirrored in ENA/EMBL. The skill warns never to send a RefSeq accession to ENA tools because they return 404.
How does it decide which sequence is best?
It applies a quality hierarchy: curated RefSeq (NM_/NP_) is preferred over predicted RefSeq (XM_/XP_), which is preferred over submitted GenBank records. It also checks version numbers and prefers the MANE Select transcript for human canonical isoforms.
What does it require to run?
It depends on the ToolUniverse Python package (from tooluniverse import ToolUniverse; tu.load_tools()) and network access to the NCBI and ENA services. It uses their documented tool functions rather than raw HTTP.
What does the output look like?
A Sequence Profile Report covering a search summary, the primary sequence's metadata and curation level, a truncated FASTA preview, annotation counts, ranked alternative sequences, cross-database references, and download-format options for FASTA or GenBank.
What are its limitations?
It is read-only retrieval and disambiguation, not analysis. Very large sequences are noted and offered as a download link rather than inlined, and results depend on the availability of the upstream NCBI/ENA services.
Install tooluniverse-sequence-retrieval
Quick Setup:
- Copy the skill folder to
.claude/skills/ - Claude will automatically detect and use the skill
Repository
mims-harvard/tooluniverse