aryarm / aryarm/as_analysis

support for single-end RNA-seq and ATAC-seq in the WASP subworkflow

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

Description

It should be possible to support single-ended data, since WASP itself supports it.
- [x] create a config parameter to allow switching to single-ended support within the WASP subworkflow
- [x] parse the samples file differently depending on the config parameter
- [x] [read the FASTQs into the same dictionary as usual](https://github.com/aryarm/as_analysis/blob/faac4573c25e61b156c5fa7640b62cd7b16db3e5/Snakefile#L16)
- [x] also [here](https://github.com/aryarm/as_analysis/blob/faac4573c25e61b156c5fa7640b62cd7b16db3e5/Snakefile#L84)
- [x] execute the aligners with only one fastq file
- [x] [STAR1](https://github.com/aryarm/as_analysis/blob/994d3f21201ddea5313c429ea44d46c968f1d905/Snakefiles/Snakefile-WASP#L48)
- [x] [STAR2](https://github.com/aryarm/as_analysis/blob/994d3f21201ddea5313c429ea44d46c968f1d905/Snakefiles/Snakefile-WASP#L135)
- [x] [BWA1](https://github.com/aryarm/as_analysis/blob/994d3f21201ddea5313c429ea44d46c968f1d905/Snakefiles/Snakefile-WASP#L67)
- [x] [BWA2](https://github.com/aryarm/as_analysis/blob/994d3f21201ddea5313c429ea44d46c968f1d905/Snakefiles/Snakefile-WASP#L154)
- [x] remove [the `--is_paired_end` parameter](https://github.com/aryarm/as_analysis/blob/faac4573c25e61b156c5fa7640b62cd7b16db3e5/Snakefiles/Snakefile-WASP#L108) from `find_intersecting_snps.py` and have it output only one fastq file
- [x] use the `rmdup.py` script instead of the [`rmdup_pe.py` script](https://github.com/aryarm/as_analysis/blob/994d3f21201ddea5313c429ea44d46c968f1d905/Snakefiles/snp2h5_rules.smk#L15) within the `rmdup` rule of the WASP subworkflow
- [ ] add a config option for [the SNP filtering criteria in `prepare_counts-rna.r` and `prepare-counts.r`](https://github.com/aryarm/as_analysis/blob/faac4573c25e61b156c5fa7640b62cd7b16db3e5/scripts/prepare_counts-rna.r#L57) (see #77)
- [ ] make it default to 10, if not provided

Contributor guide

No contributing guide indexed for this repository

Research direction

Start by reading the SNP filtering logic at the referenced line in scripts/prepare_counts-rna.r and scripts/prepare-counts.r, then review issue #77 for the intended configuration option. Done means both scripts accept the configurable SNP filtering criterion and default it to 10 when no value is provided.

Written by the indexing model from the issue text.

Assessment

Tech stack
r
Domain
bioinformatics
Issue type
Feature
Difficulty
2/5
Estimated time
1-3 hours
Activity status
Stale
Clarity
Mostly clear
Newbie friendliness
45/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.