跳到正文
cccskills
免费在 GitHub 公开

bioconductor-biostrings

Memory efficient string containers, string matching algorithms, and other utilities, for fast manipulation of large biological sequences or sets of sequences.

中文简介正在准备中。当前显示原文说明。

实测评分 未实测

查看安装方法

包含的文件(1)

  • SKILL.md5.1 KB

SKILL.md(原文)

安装前可以查看交给 AI 代理的指令内容。

Biostrings

Dependencies & Environment

Package-intrinsic requirements from the Bioconductor landing page — reproduce in any R environment.

  • Version: 2.80.1 · Bioconductor: 3.23 · R: ≥ 4.6
  • Depends: BiocGenerics, S4Vectors, IRanges, XVector, Seqinfo
  • Imports: crayon
  • Install: BiocManager::install("Biostrings")

When to Use

  • Representing and manipulating biological sequences using specialized classes like BString, DNAString, and RNAString.
  • Matching sets of probes against each other using PDict and vcountPDict.
  • Computing base content (e.g., GC content) of sequences using alphabetFrequency.
  • Handling and masking multiple sequence alignments (DNA, RNA, or amino acids) using DNAMultipleAlignment and rowmask/colmask.
  • Calculating consensus matrices and strings from masked alignments using consensusMatrix and consensusString.

When NOT to Use

  • For standard character string manipulation not involving biological sequences; use base R character vectors or stringr instead because they have less overhead.
  • For ultra-fast, genome-scale alignment of millions of short sequencing reads; use external command-line aligners (e.g., Minimap2, Bowtie2) instead.

Data Requirements

  • Input sequences in FASTA, Clustal, Phylip, or Stolkholm formats for multiple alignments.
  • Character vectors or probe sequence data packages (e.g., hgu95av2probe) for sequence manipulation.

Key Parameters

  • baseOnly: Logical to compute frequencies of only the standard bases in alphabetFrequency.
  • collapse: Logical to sum frequencies across all sequences in alphabetFrequency.
  • filepath: Path to the alignment file in readDNAMultipleAlignment.
  • format: Format of the alignment file (e.g., "clustal", "phylip").
  • invert: Logical to specify rows/columns to keep rather than hide when setting masks.
  • append: Specifies how new mask ranges interact with existing masks ("union", "intersect", "replace").
  • min.fraction: Fraction of gaps required to mask a column in maskGaps.
  • min.block.width: Minimum width of gap blocks to mask in maskGaps.

Best Practices

  • Use DNAStringSet instead of lists of DNAString objects to leverage vectorized operations.
  • Apply rowmask and colmask to MultipleAlignment objects to non-destructively hide uninformative regions or gaps before clustering.
  • Use maskGaps to automatically mask columns based on gap fractions.
  • Cast MultipleAlignment objects to DNAStringSet before exporting to FASTA with writeXStringSet.

Common Pitfalls

  • Incompatible sequence alphabets: Attempting to compare a DNAString with an RNAString directly. Fix: Convert sequences using RNAString(d) before comparison.
  • Clustering unmasked alignments: Generating phylogenetic trees with long uninformative regions, leading to inflated distances. Fix: Use maskGaps to remove gap-heavy columns before calculating stringDist.
  • Losing masks on export: Exporting to formats that don't support masks. Fix: Use write.phylip to preserve column masking, or cast to matrix to drop masked regions entirely.

Alternatives

  • seqinr: For biological sequence retrieval and analysis, but uses standard R lists/characters which are less memory-efficient.
  • ape: For phylogenetics and sequence manipulation, but lacks the specialized S4 infrastructure and non-destructive masking of Biostrings.
  • DECIPHER: For large-scale sequence alignment, database management, and oligonucleotide design, built on top of Biostrings.

Citations

  • Pagès, H. et al. (2024). Biostrings: Efficient manipulation of biological strings. R package version 2.74.1.

References

<!-- biomate-cta -->

Run this on BioMate

This skill is the knowledge layer — when, why, and how to use biostrings. To run this analysis on your own data with managed compute, automated QC, and reproducible outputs, use BioMate — free to start.

▶ Open biostrings on BioMate →

评价

还没有评价。欢迎分享你的使用感受。

评论

无需登录即可发表。欢迎分享使用感受、提问,或说明在哪些环境下无法运行。

还没有评论。来写第一条吧。

提交后会立即公开。请勿填写邮箱、电话等个人信息。URL 最多 2 个。

0 / 1000 字

此表单受 reCAPTCHA 保护,适用 Google 的隐私权政策和服务条款。

同一仓库的技能

The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题