Skip to content

Problem with cell spolit code #3

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

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions