nf-core / nf-core/atacseq

Latest version still generated MACS2, not MACS3

Open
#389 1 comment 0 reactions 0 assignees View on GitHub

Nobody has claimed this yet.

bug
Dominant language
Nextflow
Stars
229
Forks
141
Avg merge
13h 41m
Merged PRs (30d)
1

Description

Description of the bug

Hi,

I want to ask why I use latesest version, but my results is still macs2, not macs3.

Here is my command line:

nextflow run nf-core/atacseq -resume -c nextflow.config --input /project/sli68423_1316/users/yang/workspace/nanome/related_proje
ct/aml_nanopore/atac/00samplesheet_atac.csv --outdir results_mll_atac_hg38 -profile singularity --genome GRCh38 --read_length 50 -
-narrow_peak

Here is output directory:

ls bwa/merged_library/
ataqv                           D3_REP1.mLb.clN.sorted.bam      D6_REP1.mLb.clN.sorted.bam      deeptools
bigwig                          D3_REP1.mLb.clN.sorted.bam.bai  D6_REP1.mLb.clN.sorted.bam.bai  macs2
D1_REP1.mLb.clN.sorted.bam      D4_REP1.mLb.clN.sorted.bam      D7_REP1.mLb.clN.sorted.bam      picard_metrics
D1_REP1.mLb.clN.sorted.bam.bai  D4_REP1.mLb.clN.sorted.bam.bai  D7_REP1.mLb.clN.sorted.bam.bai  samtools_stats
D2_REP1.mLb.clN.sorted.bam      D5_REP1.mLb.clN.sorted.bam      D8_REP1.mLb.clN.sorted.bam
D2_REP1.mLb.clN.sorted.bam.bai  D5_REP1.mLb.clN.sorted.bam.bai  D8_REP1.mLb.clN.sorted.bam.bai
Command used and terminal output

nextflow run nf-core/atacseq -resume -c nextflow.config --input /project/sli68423_1316/users/yang/workspace/nanome/related_proje
ct/aml_nanopore/atac/00samplesheet_atac.csv --outdir results_mll_atac_hg38 -profile singularity --genome GRCh38 --read_length 50 -
-narrow_peak

Relevant files

No response

System information

No response

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 by reproducing the reported command with the supplied configuration, profile, genome, and samplesheet, then inspect the bwa/merged_library/macs2 output. Trace which peak-calling component the workflow invokes and compare that behavior with the requested latest version; done means identifying the cause and confirming the expected MACS version is used.

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
25/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.