Skip to content

Problem with cell spolit code #3

@Yassy-Ubuntu

Description

@Yassy-Ubuntu

Hi,

I was executing your code with raw data (Thp1_macro_DLO_Hi-C1_repaired_R1.fastq.gz and Thp1_macro_DLO_Hi-C1_repaired_R2.fastq.gz) from your paper (Decoding the spatial chromatin organization and dynamic epigenetic landscapes of macrophage cells during differentiation and immune activation) and I get this error, could you please help me to fix it?

Error executing process > 'split_cells (1)'

Caused by:
Missing output file(s) Thp1_macro_DLO_Hi-C1_repaired_R.cell.*.pairs expected by process split_cells (1)

Command executed:

spcell Thp1_macro_DLO_Hi-C1_repaired_R.pairs /home/msbb/Yasmine/sciDLO/test_data/barcodes.txt -o Thp1_macro_DLO_Hi-C1_repaired_R.cell -t 20

Command exit status:
0

Command output:
2023-12-19 17:18:05,545 INFO [spcell] pairs_file: Thp1_macro_DLO_Hi-C1_repaired_R.pairs barcodes_file: /home/msbb/Yasmine/sciDLO/test_data/barcodes.txt
dist_thresh: 1, max_diff_b1b2: 2, max_diff_r1r2: 4
threads: 20 wait_timeout: 2000
2023-12-19 17:18:07,732 INFO [spcell] End split cell.
2023-12-19 17:18:07,732 INFO [spcell] Count result:
valid: 0 0.00%
b1b2_not_match: 417219 99.77%
r1r2_not_match 0 0.00%
barcode_not_found 949 0.23%
total readsYas418168

Barcode counts:

Work dir:
/home/msbb/Yasmine/test_DLO/work/bd/bce4786dbf8330553c6d82f887f025

Tip: view the complete command output by changing to the process work dir and entering the command cat .command.out

P.S.
The test run with your test data worked great, so I am wondering what could be the problem.

Thank you!!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions