本文へ移動
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(原文)

インストールする前に、エージェントに与えられる指示の中身を確認できます。

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 →

レビュー

まだレビューはありません。使ってみた感想をお寄せください。

同じリポジトリのスキル

概要と使いどころ

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-kb8042026年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-kb8042026年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-kb8042026年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-kb8042026年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-kb8042026年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-kb8042026年6月21日 更新

bioMate-AI のスキルをすべて見る

このスキルの問題を報告する