Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error with snATAC peakcall: ZeroDivisionError #150

Open
itsyogeshd opened this issue Dec 22, 2021 · 0 comments
Open

Error with snATAC peakcall: ZeroDivisionError #150

itsyogeshd opened this issue Dec 22, 2021 · 0 comments

Comments

@itsyogeshd
Copy link

itsyogeshd commented Dec 22, 2021

I am trying to run the pipeline on snRNA and snATAC data generated using the 10x multi-ome kit, the scRNA part completed successfully but the snATAC part exited with the following error.

  File "/home/ydhungan/miniconda3/envs/MAESTRO/bin/macs2", line 653, in <module>
    main()
  File "/home/ydhungan/miniconda3/envs/MAESTRO/bin/macs2", line 51, in main
    run( args )
  File "/home/ydhungan/miniconda3/envs/MAESTRO/lib/python3.8/site-packages/MACS2/callpeak_cmd.py", line 64, in run
    if options.PE_MODE: (treat, control) = load_frag_files_options (options)
  File "/home/ydhungan/miniconda3/envs/MAESTRO/lib/python3.8/site-packages/MACS2/callpeak_cmd.py", line 350, in load_frag_files_options
    treat = tp.build_petrack()
  File "MACS2/IO/Parser.pyx", line 614, in MACS2.IO.Parser.BEDPEParser.build_petrack
  File "MACS2/IO/Parser.pyx", line 642, in MACS2.IO.Parser.BEDPEParser.build_petrack
ZeroDivisionError: float division
[Tue Dec 21 16:45:24 2021]
Error in rule scatac_allpeakcall:
    jobid: 22
    output: Result/Analysis/data/sample1_all_peaks.narrowPeak, Result/Analysis/data/sample1_all_treat_pileup.bdg
    log: Result/Log/sample1_macs2_allpeak.log (check log file(s) for error message)
    shell:
        macs2 callpeak -f BEDPE -g hs --outdir Result/Analysis/sample1 -n sample1_all -B -q 0.05 --nomodel --extsize=50 --SPMR --keep-dup all -t Result/Mapping/data/fragments_corrected_dedup_count.tsv
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)
@itsyogeshd itsyogeshd changed the title Error with scATAC peakcall: ZeroDivisionError Error with snATAC peakcall: ZeroDivisionError Dec 22, 2021
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant