Feature Request/Support: Support for mixed fragment_single and fragment_paired keys in AnnDataSet
Nobody has claimed this yet.
- Dominant language
- Python
- Stars
- 323
- Forks
- 43
- PR merge metrics
- No merged PRs in 30d
Description
I am attempting to combine Single-end (SE) and Paired-end (PE) ATAC-seq samples into a single AnnDataSet for downstream analysis using snapatac2. The Problem is SE samples store fragment information in adata.obsm['fragment_single'] while PE samples store fragment information in adata.obsm['fragment_paired']. When these are combined into an AnnDataSet, the disparate keys in .obsm lead to issues where downstream functions cannot locate the fragment data because they expect a consistent key across all underlying AnnData objects.
What is the best practice for joint analysis of mixed-end libraries in SnapATAC2? Is there a recommended way to "standardize" these keys before merging (e.g., treating PE as SE) without losing data integrity? Does AnnDataSet support mapping these two different keys to a universal fragment interface?
Thank you for your help!
Contributor guide
No contributing guide indexed for this repository
First steps
- Read the whole issue, then the project's contributing guide.
- Comment on the issue to say you are picking it up — it saves two people doing the same work.
- Fork the repository and make your change on a branch.
- Open a pull request that references the issue number.
Research direction
No files or tests are named. Start by tracing how AnnDataSet combines .obsm entries and how downstream functions locate fragment_single or fragment_paired. Done would require a documented, data-preserving approach for mixed single-end and paired-end libraries, with coverage for merging and downstream analysis.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- bioinformatics
- Issue type
- Feature
- Difficulty
- 5/5
- Estimated time
- Over a week
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 25/100