GangCaoLab / sciDLO

5 stars 1 forks source link

Problem with cell spolit code #3

Open Yassy-Ubuntu opened 11 months ago

Yassy-Ubuntu commented 11 months ago

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!!

Nanguage commented 11 months ago

It looks like the barcode and the sequence don't match.