aryarm / aryarm/as_analysis

allow BAM files as input to the WASP subworkflow

Open
#71 1 comment 0 reactions 0 assignees View on GitHub
enhancement gcp
Dominant language
R
Stars
10
Forks
9
PR merge metrics
No merged PRs in 30d

Description

The WASP subworkflow requires RNA-seq FASTQ files as input. But what if the alignment step has already been performed, and the user only has BAM files, instead?

We could automatically skip the first mapping step and go right to `find_intersecting_SNPs`. To ensure the second mapping step is performed identically to the first, we can add a config option allowing the user to provide the argument for the second mapping step or we could try to take it from the `@PG` tags in the BAM header.

Contributor guide

No contributing guide indexed for this repository

Research direction

Start at the WASP subworkflow and its find_intersecting_SNPs entry point; inspect how FASTQ inputs and the first mapping step are currently wired. Decide whether the second mapping argument should come from a config option or BAM @PG tags, then verify that BAM input skips the initial mapping while preserving identical second-step arguments.

Written by the indexing model from the issue text.

Assessment

Tech stack
r
Domain
bioinformatics
Issue type
Feature
Difficulty
5/5
Estimated time
Over a week
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
25/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.