nf-core / nf-core/atacseq

limiting to NFR in the dev version

Open
#424 0 comments 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

Dominant language
Nextflow
Stars
230
Forks
142
Avg merge
13h 41m
Merged PRs (30d)
1

Description

How could I run the analysis with NFR only? My pipeline_info/paramsXXX.json has the right range of minFragmentLength(0) and maxFragmentLength (120). But the input bam file for macs3 callpeak is not limited to NFR. Are there any parameter(s) that I missed? Here is my command:
nextflow run dev_version/atacseq -profile singularity --input ./fastq.csv --min_trimmed_reads 0 --aligner bowtie2 --keep_dups TRUE --narrow_peak TRUE --genome hg38 --read_length 150 --outdir ./outDir_

Thanks,

Contributor guide

Open the contributing guide

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

Research direction

Start with dev_version/atacseq and pipeline_info/paramsXXX.json, then trace how minFragmentLength and maxFragmentLength reach the macs3 callpeak input. Reproduce the command with the supplied hg38 and Bowtie2 options and verify that the input BAM contains only the requested NFR range.

Written by the indexing model from the issue text.

Assessment

Domain
bioinformatics
Issue type
Bug
Difficulty
4/5
Estimated time
3-5 days
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
30/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.