跳到正文
cccskills
免费在 GitHub 公开

bioconductor-rsamtools

This package provides an interface to the 'samtools', 'bcftools', and 'tabix' utilities for manipulating SAM (Sequence Alignment / Map), FASTA, binary variant call (BCF) and compressed indexed tab-delimited (tabix) files.

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

查看安装方法

包含的文件(1)

  • SKILL.md5.1 KB

SKILL.md(原文)

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

Rsamtools

Dependencies & Environment

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

  • Version: 2.28.0 · Bioconductor: 3.23 · R: ≥ 4.6
  • Depends: Seqinfo, GenomicRanges, Biostrings
  • Imports: BiocGenerics, S4Vectors, IRanges, XVector, bitops, BiocParallel
  • System requirements: GNU make
  • Install: BiocManager::install("Rsamtools")

When to Use

  • Low-level Alignment Import: Importing and parsing BAM files directly into R lists or DataFrame objects using scanBam.
  • Targeted Read Extraction: Extracting specific genomic coordinates and alignment fields (e.g., read sequence, strand, position) using ScanBamParam.
  • BAM File Management: Managing collections of BAM files and their metadata programmatically using BamViews.
  • Coverage Calculation: Calculating read coverage over specific genomic ranges using coverage on imported BAM data.

When NOT to Use

  • For high-level representation of gapped alignments, use GenomicAlignments (specifically readGAlignments) because Rsamtools provides lower-level list-based outputs.
  • For iterating through BAM files in parallel, use GenomicFiles because it implements higher-level strategies for file iteration.
  • For whole-genome manipulation and DNA sequence operations, use Biostrings or BSgenome instead.

Data Requirements

  • Input Format: BAM files (e.g., .bam) and their corresponding index files (.bai).
  • Genomic Coordinates: Specified as GRanges objects for targeted queries.

Key Parameters

  • file: Path to the BAM file to be parsed by scanBam.
  • param: A ScanBamParam object determining which genomic coordinates and components to input.
  • which: A GRanges object specifying the genomic ranges to extract.
  • what: A character vector (e.g., c("rname", "pos", "qwidth")) specifying which BAM record fields to retrieve.
  • flag: A filter created by scanBamFlag (e.g., isUnmappedQuery=FALSE) to restrict reads.
  • bamRanges: A GRanges object defining the genomic rows for a BamViews instance.

Best Practices

  • Always ensure a .bai index file is available locally when querying specific ranges with the which argument.
  • Restrict imported fields using the what argument in ScanBamParam to minimize memory consumption.
  • Process large BAM files in chunks (e.g., by chromosome) using summary functions to avoid memory exhaustion.
  • Use BamViews to marshal references to multiple BAM files without loading them immediately into memory.

Common Pitfalls

  • Querying by range without an index: Fails because providing a which argument to ScanBamParam requires a BAM index. Fix: Ensure the .bai file exists or create one using indexBam.
  • "EOF marker is absent" warning: Occurs with some BAM files. Fix: This message can safely be ignored as noted in the vignette.
  • Chromosome naming mismatches: Ensembl vs UCSC naming (e.g., "chr1" vs "1") causes empty queries. Fix: Use seqlevels<- to map between naming schemes before querying.

Alternatives

  • GenomicAlignments: Provides a more useful, higher-level representation of BAM files in R (e.g., GAlignments objects).
  • GenomicFiles: Useful for iterating through BAM and other files, including in parallel.
  • ShortRead: Better suited for I/O and quality assessment of ungapped short read alignments.

Citations

  • Morgan M, Pagès H, Obenchain V, Hayden N. Rsamtools: Binary alignment (BAM), FASTA, variant call (BCF), and tabix file import. Bioconductor.

References

<!-- biomate-cta -->

Run this on BioMate

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

▶ Open rsamtools 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 的全部技能

报告此技能的问题